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FOREWORD 

This  report  provides  results  of  a  detailed  study  to  estimate  lifetime 
costs  for  flexible  pavements  as  a  function  of  legal  axle  limits  using 
an  improved  version  of  FHWA  Computer  Program  VESYS  IIM. 

VESYS  IIM  was  modified  to  include  capabilities  for:  1)  seasonal 
characterizations  of  pavement  materials,  2)  a  discretized  representation 
of  axle  load  distribution,  and  3)  low-temperature  cracking  predictions. 

A  combined  survey  of  the  literature  and  laboratory  testing  program  was 
conducted  to  define  the  variations  in  material  permanent  deformation 
parameters  as  important  material  characteristics  vary  seasonally  with 
the  environment.  The  resulting  data  and  other  information  and  experience 
were  applied  to  arrive  at  input  data  that  would  yield  realistic  performance 
predictions. 

A  factorial  of  64  solutions  was  obtained  using  the  improved  version  of 
VESYS  IIM  and  the  input  data  developed  to  study  the  effects  of  truck 
traffic  consistent  with  four  levels  of  legal  axle  limits  (18,  20,  22,  and 
24  kips),  two  levels  of  traffic,  two  levels  of  pavement  section  thickness 
and  four  environmental  zones  (Wet-Freeze,  Dry-Freeze,  Wet-No  Freeze  and 
Dry-No  Freeze).  Where  failures  were  predicted,  an  overlay  was  applied 
and  a  new  solution  obtained  until  a  pavement  life  of  at  least  20  years  was 
attained.  The  initial  and  overlay  costs  were  estimated  and  these  costs 
for  20  years  of  pavement  service  were  related  to  the  legal  axle  limits. 

Copies  of  the  report,  intended  primarily  for  research  and  development 
audiences,  are  being  distributed  by  transmittal  memorandum  to  such 
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COMMENTS  BY  FHWA 


The  improved  VESYS  structural  subsystem  generated  in  this  research  study, 
also  referred  to  as  VESYS  A,  was  verified  with  comparisons  of  predicted 
and  actual  distresses  in  twelve  flexible  pavement  test  sections  throughout 
the  country.   Four  of  these  sections  were  subjected  to  an  indepth  analysis 
of  the  effect  of  increasing  axle  weights  representing  the  following  four 
different  environmental  zones: 

FREEZE  WET  Illinois  IH-80,  Section  38 

FREEZE  DRY  Utah  IH-80 

NO  FREEZE  WET  Florida  1-10 

NO  FREEZE  DRY  Texas  IH-20 

The  lifetimes  and  costs  generated  for  these  environmental  zones  reflect 
the  behavior  of  the  four  specific  test  sections  only.   These  specific 
conditions  should  not  be  associated  with  conditions  existing  in  other 
environmental  zones . 

The  study  demonstrates  the  usefulness  of  predictive  type  models  for  future 
planning  of  roadway  systems.  The  criteria  for  overlay,  costs,  etc.,  are 
based  on  the  best  data  available.   However,  additional  information  is 
required  to  provide  a  more  precise  estimate  of  pavement  performance. 

Research  and  development  of  VESYS  A  is  continuing  to  expand  its  present 
analysis  capability  from  three  to  seven  layers.   This  will  lead  to  develop- 
ment of  the  new  generation  of  the  VESYS  system,  VESYS  III. 
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CHAPTER  I 
INTRODUCTION 
Background 


Recent  Federal  Legislation  (Title  23,  Sec.  127  as  amended   1975) 
has  permitted  States  to  raise  allowable  axle  loads  on  Interstate  highways 
from  18,000/32,000  lb  (8,160/14,510  kg)   to  20,000/34,000  lb  (9,070/15,420 
kg)  on  single/dual  axles,  respectively.   There  is  considerable  concern 
for  pavement  damage  in  those  States  where  this  option  has  been  exercised; 
many  States  already  had  legal  loads  of  20,000  lb  (9,080  kg)  or  more 
under  the  "grandfather"  clause. 

A  considerable  effort  has  been  expended  already  to  predict  service 
life  reductions  as  a  function  of  increases  in  axle  load  limits  based 
on  data  from  the  AASHO  Road  Test,  Loadometer  studies  and  traffic  count 
surveys.   This  methodology  has  developed  a  predicted  64  percent  increase 
in  the  equivalent  effect  on  a  pavement  for  a  fully  loaded  truck  with  an 
axle  weighing  2,000  lbs.   (908  kg)  more  than  the  present  limits. 
However,  these  predictions  are  extrapolated  past  the  range  in  numbers 
of  axle  loads  for  which  the  data  were  developed  and  are  possibly 
suspect  as  a  result. 

There  are  at  present  several  FHWA  research  projects  underway 
aimed  at  producing  the  needed  information  using  a  range  of  techniques. 
FHWA  Research  Project  "Effects  on  Highways  of  Increased  Legal 
Vehicle  Weights  Using  VESYS  IIM"  is  one  of  these  and  has  as  its 
objectives:   (1)  modification  and  improvement  of  Computer  Program 
VESYS  IIM  so  that  it  can  be  used  in  a  realistic  manner  to  establish 
relationships  between  increased  legal  axle  loadings  and  flexible 
pavement  life  cycle  costs,  and  (2)  application  of  VESYS  IIM  to  pro- 
duce these  relationships. 

VESYS  IIM  as  modified  (modified  program  will  subsequently  be 
called  VESYS  A  for  simplicity)  is  a  sophisticated  computer  code  that 
accepts  some  23  control  variables  and  44  independent  variables  de- 
scribing a  flexible  pavement  structure,  the  traffic  loadings  it  en- 
dures and  indirectly  (through  input  of  pavement  temperatures  and 
seasonal  materials  characterizations)   the  environment  in  which  it 
exists.   Using  this  input,  it  then  predicts  fatigue  cracking,  rut 
depth,  slope  variance,  Present  Serviceability  Index  and  expected  life 
as  functions  of  time  in  terms  of  truck  traffic  and  axle  load  distribu- 
tion.  VESYS  IIM  in  the  format  extant  prior  to  this  project  is  describ- 
ed in  varying  degrees  of  detail  in  a  number  of  references  including 


1      2 
References  1  and  2  .   The  modifications  and  improvements  in  VESYS 

resulting  from  this  project  will  be  described  subsequently. 

Program  Modifications  to  VESYS  IIM 

The  modifications  to  improve  VESYS  IIM  toward  more  accurate 
simulation  of  real  pavement  structures  included  (1)  discretized 
representation  of  axle  load  distribution,  (2)  seasonal  characteriza- 
tions of  stiffness  and  permanent  deformation  properties  of  materials, 
and  (3)  low-temperature  cracking  predictions. 

VESYS  IIM  introduced  the  entire  spectrum  of  axle  loads  experi- 
enced by  the  pavement  as  a  mean  axle  load  and  a  variance.   This  was 
considered  to  be  inadequate  for  discriminating  among  the  differences 
due  to  variations  in  the  small  percent  of  higher  axle  loads,  so  a 
capability  was  introduced  for  considering  the  effects  due  to  dis- 
crete segments  of  the  axle  load  distribution  separately  and  com- 
bining them.   This  allows  the  theoretical  model  to  more  accurately 
assess  the  responses  of  the  pavement  structure  to  the  various  mag- 
nitudes of  axle  loading. 

All  analytical  procedures  based  on  elastic  layer  theory  require 
definition  of  the  stiffnesses  of  the  various  layers  as  input  data. 
Most  operate  only  with  single  values  for  each  layer,  but  VESYS  IIM 
has  the  additional  capability  of  considering  the  effects  of  tem- 
perature variations  on  surface  layer  stiffness  and  of  including 
these  variations  in  stiffness  in  the  analysis  as  damage  is  accu- 
mulated in  time.   As  the  stiffnesses  of  base  and  subgrade  materials 
also  undergo  seasonal  changes  in  moisture  content  and  density  due 
to  variations  in  rainfall,  frost  penetration,  spring  breakup  and 
other  such  causes,  a  more  accurate  simulation  results  if  they  also 
change  rationally  on  a  seasonal  basis. 

VESYS  IIM  included  a  capability  for  rationally  predicting  rut 
depts  in  time,  partially  based  on  permanent  deformation  characteriza- 
tions of  the  three  layer  materials  in  the  form  of  two  new  parameters 
ALPHA  (N)  and  GNU  (N)  for  N  =  1,  2,  and  3  (See  Reference  1  for 


Rauhut,  J.  B. ,  J.  C.  O'Quin  and  W.  R.  Hudson,  "Sensitivity 
Analysis  of  FHWA  Structural  Model  VESYS  IIM,  Vol.  1,  Preparatory 
and  Related  Studies,"  Report  No.  FHWA-RD-76-23,  March  1976. 

2 
Kenis,  W.  J.,  "Predictive  Design  Procedures  -  A  Design  Method 

for  Flexible  Pavements  Using  the  VESYS  Structural  Subsystem,"  Pro- 
ceedings, Fourth  International  Conference,  Structural  Design  of 
Asphalt  Pavements,  Volume  I,  August  1977. 


detailed  discussion).  As  the  magnitudes  of  these  parameters  are 
heavily  dependent  on  such  seasonal  variables  as  moisture  content 
and  temperature  and  are  themselves  very  important  to  the  predic- 
tions of  rut  depth,  slope  variance,  present  serviceability  index 
and  expected  life,  arrays  of  values  instead  of  single  values  are 
required  to  rationally  characterize  these  materials  across  the 
seasons  of  the  year.  VESYS  IIM  was  modified  in  this  project  to 
allow  this  capability. 

Cracking  of  the  asphalt  concrete  surface  layer  due  to  traffic- 
induced  fatigue  is  a  major  form  of  flexible  pavement  distress  and 
this  phenomena  is  modeled  in  VESYS  IIM.   However,  low-temperature 
transverse  cracking  is  another  important  form  of  cracking  distress 
that  sometimes  predominates  over  fatigue  cracking  in  northern 
states.   It  was  desirable  to  broaden  VESYS  IIM  to  predict  low- 
temperature  cracking  and  a  statistical  model  was  selected  and  added 
to  VESYS  IIM  during  this  project. 

These  three  additions  in  VESYS  IIM  capabilities  are  believed 
to  have  considerably  improved  its  ability  to  simulate  the  per- 
formance of  real  flexible  pavements. 

Development  of  Input  Data  for  Modified 
VESYS  IIM  Analyses 

The  improvements  to  Program  VESYS  IIM  resulting  in  Program 
VESYS  A  discussed  above  may  not  be  employed  effectively  unless  the 
magnitudes  of  the  input  parameters  are  known.   There  is  a  very 
considerable  body  of  information  available  for  changes  in  base  and 
subgrade  material  stiffness  with  changes  in  moisture  content  and 
dry  density,  but  there  is  little  available  information  on  the 
permanent  deformation  parameters  ALPHA  (N)  and  GNU  (N)  due  to 
their  recent  definition  and  application  to  rutting  predictions. 

Review  of  available  information  on  permanent  deformations  in 
well  compacted  base  materials  indicated  that  a  reasonable  level  of 
information  was  available  in  the  literature  (References  1,  3^,  and  42) 


Proceedings  from  a  Conference  on  Utilization  of  Graded  Aggre- 
gate Base  Materials  in  Flexible  Pavements,  held  March  25-26,  1974 
in  Oakbrook,  Illinois,  published  by  the  National  Crushed  Stone  Asso- 
ciation, National  Sand  and  Gravel  Association  and  the  National  Slag 
Association. 

Kalcheff ,  I.  V,  and  R.  G.  Hicks,  "A  Test  Procedure  for  Deter- 
mining the  Resilient  Properties  of  Granular  Materials,"  Journal  of 
Testing  and  Evaluation,  Vol.  I,  No.  6,  ASTM,  1973. 


Consequently,  it  was  decided  that  additional  testing  of  base 
materials  to  establish  seasonal  trend  information  could  be  omitted 
in  favor  of  testing  of  typical  subgrade  and  asphalt  concrete 
specimens. 

Testing  was  conducted  on  a  typical  silty  clay  soil  in  a 
small  factorial  of  varying  moisture  content  and  dry  density  to  pro- 
vide some  insight  as  to  variation  of  ALPHA  (3)  and 'GNU  (3)  with 
these  parameters.   Similarly,  asphalt  concrete  specimens  were 
tested  with  variations  in  temperature  and  stress  level  to  gain 
insight  as  to  variation  of  ALPHA  (1)  and  GNU  (1)  with  those 
parameters.   The  results  of  this  testing,  research  of  the  litera- 
ture and  previous  experience  were  used  to  assign  rational  values 
of  ALPHA  (N)  and  GNU  (N)  for  the  materials  in  the  three  layers 
and  for  the  four  seasons. 

The  testing  performed,  its  results  and  results  of  the  litera- 
ture research  are  discussed  in  detail  in  the  Appendices.   The 
resulting  values  selected  appear  and  are  also  discussed  in  Chap- 
ters IV  and  V. 

Much  of  the  input  data  resulted  from  selection  of  the  four 
interstate  highways  to  represent  the  four  environments  wet-freeze, 
dry-freeze,  wet-no  freeze  and  dry-no  freeze  and  was  furnished  by 
personnel  of  the  various  State  Departments  of  Transportation  and 
other  sources.   Other  data  such  as  pavement  temperatures  and  axle 
load  distributions  for  various  legal  axle  limits  required  some 
research  to  develop  and  their  development  is  discussed  in  more 
detail  subsequently. 

Input  data  was  developed  for  four  test  sections  each  from  the 
AASHO  Road  Test  and  the  Brampton  Test  Road  as  well  as  the  four 
interstate  highways.   As  previous  comparisons  had  been  made  between 
VESYS  IIM  predictions  and  measured  performance  for  the  same  AASHO 
Road  Test  and  Brampton  Test  Road  sections  (Reference  1) ,  it  was 
possible  to  observe  the  effects  of  using  the  improved  program.   It 
was  also  possible  and  necessary  to  compare  solutions  for  the  four 
interstate  highway  sections  to  performance  observations  made  avail- 
able by  state  personnel.   Where  predictions  appeared  errant,  the 
input  values  were  studied  to  isolate  inadequacies  and  to  correct 
the  input  values  where  more  rational  values  could  be  developed  by 
further  study.   This  iterative  approach  of  comparing  predictions 
and  observed  performance,  studying  the  input  parameters  apparently 
causing  differences  to  see  if  the  input  values  were  adequate  or 
could  be  rationally  improved,  new  predictions  and  comparisons,  etc., 
was  used  to  effect  a  "calibration"  of  the  procedure  so  that  rational 
predictions  could  be  expected  for  the  test  sections  to  be  studied. 


Relating  Increased  Legal  Axle  Loads  to 
Flexible  Pavement  Life  Cycle  Costs 

As  previously  stated,  the  goals  of  this  research  project  were 
to  effect  certain  improvements  to  VESYS  IIM  so  that  it  could  be  used 
to  relate  increased  legal  axle  loads  to  flexible  pavement  life 
cycle  costs  and  to  apply  these  capabilities  to  that  task.   The 
application  of  the  improved  program  was  made  in  a  factorial  of  64 
solutions  including  the  four  interstate  highways,  two  levels  of 
pavement  section  thickness,  two  levels  of  traffic  and  the  four 
levels  of  legal  axle  limits.   Where  either  predicted  cracking  or 
rutting  caused  failure  before  the  twenty-year  design  life  was 
reached,  overlays  were  added  and  new  solutions  obtained.   This 
effort  is  discussed  in  Chapter  VI. 

Results  in  the  form  of  performance  and  cost  data  to  provide 
satisfactory  service  across  the  design  life  of  the  pavement  have 
been  accumulated  and  are  presented  in  Chapter  VII. 


CHAPTER  II 
PROGRAM  MODIFICATIONS  TO  VESYS  IIM 


A  complete  description  of  the  VESYS  IIM  program  can  be  found 
in  Reference  1;  this  should  be  read  prior  to  reading  the  following 
description  of  changes. 

The  modifications  performed  on  the  program  under  this  project 
fall  into  three  major  categories:   (1)  the  seasonal  modification 
of  materials  properties;  (2)  the  breakdown  of  the  load  distribution 
by  tire  radius  and  corresponding  tire  pressure;  and  (3)  the  addition 
of  low  temperature  cracking  as  a  distress  mechanism.   The  first  two 
involved  additional  static  solutions  and  thus  required  many  changes 
in  dimensions  and  indexing,  whereas  the  third  was  a  relatively  sim- 
ple addition  to  the  repeated  load  portion  of  the  program.   These 
three  changes  will  be  discussed  separately  below. 

Seasonal  Characterization  of  Pavement  Materials 

VESYS  A  considers  a  three  layer  pavement  system:   an  asphalt 
layer,  a  base,  and  a  subgrade.   For  each  layer  the  permanent  deforma- 
tion parameters  are  input  for  each  season,  instead  of  remaining  con- 
stant throughout  the  year;  in  addition,  the  creep  compliance  func- 
tions for  the  base  and  subgrade  layers  are  modified  by  separate  sea- 
sonal multipliers  for  each  layer  and  season.   These  seasonal  multi- 
pliers permit  modeling  of  the  effects  of  environmental  variables  on 
the  stiffness  of  the  layers  (e.g.,  spring  breakup  of  frozen  soil). 

The  implementation  of  these  additional  parameters  requires  that 
a  static  solution  be  performed  for  each  season,  whereas  before  only 
one  static  solution  was  required  for  the  entire  simulation.   A  Type  I 
"static  solution"  for  the  VESYS  system  means  the  solution  for  the 
set  of  desired  responses  to  a  loading  of  1  psi  (  6.9  kN/m  )  over  a 
circle  of  specified  radius  (permanent  deformation  at  the  surface  and 
radial  strain  at  the  bottom  of  the  surface  layer)  as  a  function  of  load 
duration,  these  responses  expressed  as  sets  of  coefficients  of  a 
Dirichlet  series.   Because  of  this  increase  in  the  number  of  static 
solutions,  many  changes  in  indexing  and  in  array  dimensions  were  re- 
quired; these  changes  will  be  detailed  later. 

Since  there  was  already  provision  for  up  to  24  seasons,  i.e., 
periods  in  the  year  for  which  a  temperature  could  be  input,  the 
number  of  possible  seasons,  or  static  solutions,  was  maintained  at 
24  for  future  flexibility.   This  many  seasons  would  be  too  expensive 
to  run  in  a  large  factorial  for  a  limited  research  effort  and  is  not 
necessary.   Use  of  input  values  averaged  over  three-month  periods 
representing  the  four  seasons  of  the  year  may  be  expected  to  provide 
characterizations  whose  accuracy  is  consistent  with  the  "state-of-the- 


art"  in  materials  characterization  and  model  simulations.   The  decision 
was  also  made  at  this  time  based  on  previous  experience  with  the  program 
to  replace  the  capability  of  executing  a  Type  3  run  (repeated  load 
solution  only)  based  on  card  input  of  results  from  a  previous  Type  1 
run  by  a  similar  capability  to  obtain  these  results  from  a  computer  file. 
It  was  decided  that  too  many  chances  for  error  existed  in  the  selection 
of  input  values;  and  also,  with  the  addition  of  multiple  static  solu- 
tions, that  the  volume  of  input  required  was  too  great  to  be  done  con- 
veniently with  cards.   This  decision  made  possible  the  removal  of  some 
of  the  input  variables  and  keywords,  and  thus  the  program  has  now  after 
modification  approximately  the  same  number  of  input  variables  (physical 
and  control)  as  before. 

A  list  of  those  key  words  added  and  those  removed  will  appear 
at  the  end  of  this  chapter,  as  well  as  a  list  of  variables  added 
and  those  whose  dimensions  were  changed. 

Changes  to  the  main  program  involved  additions  to  the  input 
guide,  removal  of  the  read  statements  corresponding  to  the  key 
words  removed,  modifications  to  the  read  statements  which  now  read 
data  for  each  season  (ALPH,  GGU,  and  SMLT)  and  corresponding  changes 
to  the  printout  and  data  validity  checks,  and  insertion  in  this 
routine  of  the  call  to  the  master  repeated  load  subroutine  for  the 
second  part  of  a  TYPE  1  run  (this  had  been  in  subroutine  STATIC, 
although  the  call  to  the  same  routine  in  the  case  of  a  TYPE  3  run 
was  in  the  main  program) . 

Subroutine  STATIC  is  the  main  subroutine  for  the  static  load 
portion  of  the  program.   Changes  to  this  routine  involve,  beyond 
the  changes  in  dimension  and  common  blocks,  the  addition  of  a  new 
loop  around  most  of  the  routine,  the  index  on  which  serves  to 
identify  the  season  for  which  a  solution  is  being  made.   All 
references  to  the  values  of  input  permanent  deformation  parameters 
ALPHA  and  GNU  must  be  made  by  season  as  well  as  by  layer,  and  the 
results  of  the  solutions  (the  coefficients  for  the  series  represent- 
ing the  variation  of  response  with  load  duration,  and  coefficients 
of  variation  of  these  series)  are  also  stored  by  season;  the  two 
quantities  ®SYS   and  ysys,  curve  fitting  parameters  representing 
the  permanent  deformation  performance  of  the  entire  pavement  sys- 
tem, also  are  stored  by  season.   In  performing  the  solutions  for 
each  season,  the  coefficients  of  the  series  representing  creep  com- 
pliance for  the  base  and  subgrade  are  multiplied  by  a  seasonal  mul- 
tiplier for  the  given  layer  (multiplying  each  coefficient  by  a  con- 
stant has  the  effect  of  multiplying  the  value  of  the  evaluated 
series  by  the  same  constant) .   This  modified  creep  compliance  func- 
tion is  used  then  in  the  calculations  for  the  desired  responses 
for  that  season  (radial  strain  at  the  bottom  and  deformation  at 
the  top  of  the  asphalt  surface  layer) ;  the  coefficients  are  then 
returned  to  their  original  values. 


Another  change,  not  related  to  seasonal  variation  of  input 
parameters,  is  the  revision  of  the  limits  on  the  loop  performed  to 
evaluate  the  system  permanent  deformation  parameters;  this  is  the 
"unload"  portion  of  the  static  solution  for  vertical  deformation, 
in  which  the  layer  moduli  or  creep  compliances  are  modified  accord- 
ing to  (1)  the  permanent  deformation  parameters  for  the  particular 
layer  and  season  and  (2)  the  number  of  loads  for  which  the  simula- 
tion is  desired.   It  is  felt  that  since  the  actual  number  of  load- 
ings considered  is  large,  the  values  of  repeated  axle  loads  used 
to  obtain  aSyS  and  ySvs  as  a  straight  line  fit  to  the  logarithms 
of  predicted  surface  deformation  as  a  function  of  the  logarithms  of 
the  number  of  axle  loads  should  be  correspondingly  large.   Conse- 
quently, the  values  for  LY1  and  LY2  used  throughout  this  project 
were  6  and  8,  corresponding  to  values  for  numbers  of  axle  loadings 
of  1()5,  106,  and  10?,  respectively.   These  values  were  suggested  by 
ARE  in  a  previous  study  (Ref.  1)  but  were  replaced  at  FHWA  by 
values  for  LY1  and  LY2  of  1  and  5 ,  corresponding  to  a  range  in 
loadings  from  1  to  10^,  due  to  numerical  difficulties  encountered 
on  the  computer  in  use  there. 

In  the  repeated  (or  random)  load  portion  of  VESYS  A,  subrou-  - 
tine  RANDOM  serves  as  a  calling  routine  for  the.  rest  of  the  section, 
as  well  as  performing  certain  calculations  needed  for  both  the  rut 
depth  and  the  fatigue  cracking  calculations.   In  this  routine, 
unnecessary  common  blocks  were  removed,  the  variable  BYPASS  was 
added  to  the  calling  sequence,  and  an  error  printout  previously  in 
subroutine  ERROR  was  added  here  instead  so  that  in  case  an  error 
occurs,  the  user  can  tell  for  which  season  (which  static  solution) 
the  problem  exists.   (The  errors  considered  are:   (1)  a  negative 
sum  of  Dirichlet  series  coefficients  for  a  given  response  and  sea- 
son; and  (2)  a  negative  value  in  the  array  DELTA,  which  contains 
values  used  in  evaluating  these  series.)   Calls  to  subroutines 
CRACKS  and  ROUGH  were  removed  from  subroutine  SERVC  and  added  to 
RANDOM  to  clarify  the  sequence  of  operations. 

Subroutines  ROUGH  and  CRACKS  have  several  features  in  common; 
a  short  discussion  here  will  save  duplication  later.   ROUGH  deals 
with  surface  deformation;  CRACKS  deals  with  radial  strain  at  the 
bottom  of  the  asphalt;  they  each  evaluate  the  physical  result 
(deformation  or  strain)  at  each  season  using  the  appropriate 
Dirichlet  series,  the  value  of  input  load  duration,  and  the  time- 
temperature  shift  calculated  in  RANDOM  for  the  temperature  cor- 
responding to  that  season.   The  desired  result  and  its  variance  are 
calculated  in  subroutine  GRESP  (Generalized  RESPonse)  and  are 
printed.   In  CRACKS  this  result  (axial  strain)  is  used  to  calcu- 
late the  damage  per  cycle  of  loading,  while  in  ROUGH  the  surface 
deflection  is  used  more  directly  in  the  next  portion  of  the  rou- 
tine. 


In  each  routine  a  loop  is  then  performed  over  the  number  of 
seasons  required  to  attain  the  last  time  input  in  the  TRANDOM  array. 
After  the  accumulated  damage  is  calculated  and  stored  for  each  of 
these  seasons,  a  second  loop  is  performed  over  the  number  of  TRANDOM 
values  for  which  output  is  desired.   At  each  value,  the  damage  is 
linearly  interpolated  between  the  two  closest  "season"  times,  where 
the  time  associated  with  each  season  is  the  end  of  that  season. 

Subroutine  ROUGH,  which  calculates  the  accumulated  permanent 
deformation  (rutting)  and  associated  statistical  quantities  leading 
to  slope  variance,  was  not  greatly  changed.   The  common  blocks  were 
modified  as  previously  discussed;  the  variables  ALPHA  and  GNU  (pre- 
viously the  single  system  permanent  deformation  parameters)  were 
set  successively  to  values  for  the.  current  season  from  the  arrays 
ALF  and  GMU  in  the  process  of  generating  the  rutting. 

The  rutting  calculation  was  not  changed  except  for  the  indexing 
of  variables  which  previously  had  only  one  value.   The  interpolation 
routine  was  changed  so  that  it  was  no  longer  possible  to  reference 
the  rutting  for  a  season  beyond  those  for  which  it  had  been  calcu- 
lated.  Under  certain  conditions  (not  unusual),  this  could  and  did 
happen  previously;  the  results  were  not  affected  but  the  reference 
made  debugging  elsewhere  more  difficult  (a  common  debugging  tool  is 
to  preset  the  computer  memory  to  a  value  which  will  cause  an  error 
exit  if  the  value  is  used;  this  helps  to  isolate  undefined  variables 
and  the  use  of  undefined  locations  in  arrays). 

The  only  remaining  change  in  ROUGH  was  to  define  GVNRM,  a  mul- 
tiplier of  the  variance  of  rut  depth  in  the  calculation  of  slope 
variance,  to  be  the  RMS  average  for  all  seasons  of  the  values  of 
GCVAR  (K, 4)  ,  the  coefficients  of  variation  of  the  system's  vertical 
deformation.   This  choice  was  made  rather  than  a  straight  mean 
because  it  is  the  square  of  GVNRM  that  is  actually  used  in  the  cal- 
culation of  slope  variance. 

Changes  to  subroutine  CRACKS  were  very  similar  to  those 
applied  to  ROUGH.   There  is  no  reference  to  seasonally  indexed 
quantities  except  for  the  strain  and  its  variance,  however,  so  the 
changes  are  fewer  in  number.   The  interpolation  for  the  damage  and 
its  variance  at  time  TRANDOM  is  performed  in  the  same  way  as  for 
subroutine  ROUGH, 

Subroutine  GRESP  was  modified  to  include  the  seasonal  index  K 
on  the  Dirichlet  series  Gs  viz.  G(K,  I);  G  is  either  GNORM  (in  the 
call  from  ROUGH)  or  GRAB  (in  the  call  from  CRACKS).  GRESP  already 
returned  results  for  each  season  due  to  the  inclusion  of  the  time- 
temperature  shift.  The  internal  dimension  statement  for  G  was 
changed  from  G(l) ,  simply  indicating  that  G  is  a  vector,  to 


G(25,  20).   Another  formal  parameter,  S,  previously  undlmensioned, 
was  changed  to  S(25),  and  all  references  to  S  were  changed  to 
S(K).   (S  corresponds  to  the  variables  GCVAR  (K,4)  or  GCVAR  (K,7), 
representing  the  coefficient  of  variation  of  the  Dirichlet  series 
coefficients. ) 

The  only  other  change  was  the  removal  of  calls  to  subroutines 
ROUGH  and  CRACKS  from  SERVC,  corresponding  to  their  insertion  in 
subroutine  RANDOM. 

Discretized  Representation  of  Axle  Load  Distributions 

VESYS  IIM  uses  just  three  parameters  to  describe  distribution 
of  wheel  load  (taken  as  one  half  of  axle  load)  in  the  flow  of 
traffic:   RADIUS,  the  radius  of  a  circle  of  area  equal  to  that  of 
the  tire  footprint;  AMPLITUD,  strictly  the  multiplier  of  the  value 
of  LOADING  used  in  the  static  solutions,  the  product  then  represent- 
ing tire  pressure  (in  practice,  since  a  value  of  1.0  is  normally 
used  for  LOADING,  and  is  the  default  value,  AMPLITUD  may  be  more 
simply  thought  of  as  the  actual  tire  pressure);  and  VCAMP,  the 
variance  of  AMPLITUD.   AMPLITUD  (for  LOADING  =  1.0,  as  will  be 
assumed  throughout  the  remainder  of  this  report)  and  RADIUS  together 
define  a  specific  wheel  load,  given  by 

wheel  load  =  IT  (RADIUS)2  •  AMPLITUD 

VCAMP  then  must  represent  the  entire  variation  of  wheel  (or  axle) 
load,  whereas  in  fact  the  loaded  area  also  varies  with  wheel  load. 
Because  the  histogram  of  wheel  load  is  often  both  skew  and  double 
humped,  the  representation  of  this  histogram  simply  by  a  mean  and 
a  variance  is  poor. 

Initially  it  was  planned  simply  to  input  a  series  of  AMPLITUD 
values  and  to  represent  the  variation  in  wheel  load  this  way,  but 
these  values  become  unrealistic  when  tied  to  a  single  radius,  so 
it  was  decided  to  allow  static  solutions  to  be  made  for  more  than 
one  radius  value.   An  axle  load  distribution  can  be  broken  into 
several  parts  and  the  mean  load  determined  for  each  part,  which, 
with  the  assumption  of  a  reasonable  tire  pressure,  leads  to  a  value 
of  RADIUS  for  that  portion  of  the  load  distribution.   This  pro- 
cedure not  only  permits  a  better,  more  reasonable  representation  of 
load  distribution,  but  it  also  permits  isolation  of  the  effect  of 
different  portions  of  that  load  distribution  on  the  several  distress 
manifestations.   Of  course,  there  is  considerable  extra  computa- 
tional cost  involved  in  performing  multi-radius  solutions  over 
those  for  a  single  radius  since  the  cost  of  a  simulation  run  is 
roughly  proportional  to  the  total  number  of  static  solutions  (num- 
ber of  seasons  times  number  of  radii) . 
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As  might  be  expected,  the  implementation  of  multiple-radius 
simulations  involved  changes  similar  to  those  previously  mentioned 
concerning  the  seasonal  variation  of  materials  properties:   changes/ 
additions  to  input  variables,  common  blocks,  indexing,  and  output. 
In  both  cases,  the  major  change  is  the  increase  in  the  number  of 
static  solutions  and  associated  storage  and  indexing  (bookkeeping). 
The  changes  required  in  each  major  routine  or  group  of  routines  are 
outlined  below. 

In  VESYS,  the  additions  consisted  of  several  new  input  param- 
eters, new  variables  in  common  blocks,  and  removal  or  changes  in 
the  printed  output.   The  variables  added  are  listed  in  Table  1  and 
include:   (1)  NRADIUS  (RADIUS  and  AMPLTTUD  are  now  vectors  if 
NRADIUS  >  1);  (2)  LOADDIST,  which  gives  the  fraction  (not  percent) 
of  the  total  axles  per  day  to  be  associated  with  a  particular  wheel 
load  given  by  RADIUS  and  AMPLITUD  and  (3)  TRAFMULT,  which  permits  a 
seasonal  variation  in  the  traffic.   A  new  common  block,  /DIST/,  was 
added  to  hold  these  new  vectors  and  the  value  of  NRADIUS,  which  is 
limited  to  1^  NRADIUS  -12  (default  is  1).   Changes  in  the  lists  of 
default  and  limiting  values  as  well  as  in  the  list  of  acceptable  key 
words  were  made  in  BLOCKDATA.   The  new  arrays  were  also  preset  to 
appropriate  values. 

In  subroutine  STATIC,  changes  beyond  dimensions  and  common 
blocks  involved  primarily  indexing  and  access  to  some  of  the  new 
variables.   Since  in  the  previous  set  of  modifications  a  loop  had 
been  added  around  the  static  solution  calculation,  the  variation  of 
radius  was  handled  in  the  following  manner:   The  total  number 
(variable  name  NTXNR)  of  static  solutions  required  (equal  to 
NRADIUS  radii  times  NTEMPS  seasons,  using  keywords  rather  than 
actual  variable  names)  was  calculated  and  used  as  the  limit  of  the 
loop  on  static  solutions  in  place  of  the  number  of  seasons  alone. 
For  each  season,  NRADIUS  static  solutions  were  performed;  the  vari- 
able AAA,  previously  used  for  the  single  input  radius,  was  set  in 
turn  equal  to  each  value  in  the  array  RADIUS. 

The  indexing  was  done  as  follows:   ITM  is  the  sequential  index 

on  static  solutions  and  now  runs  from  1  to  NTXNR.  ITEMP  and  IRAD, 

the  indices  for  season  and  radius  corresponding  to  a  given  static 
solution,  are  calculated  by 

ITEMP  =  (ITM-1)/NRAD  +1 

and 

IRAD  =  M0D((ITM-1),  NRAD)  +1 
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where  NRAD  is  the  variable  containing  the  number  of  radii,  MOD  is 
the  modulo  arithmetic  function,  and  use  is  made  of  the  fact  that 
results  in  integer  arithmetic  in  FORTRAN  are  always  truncated. 
Assuming,  for  example,  a  simple  case,  NRAD  =  3,  NTEMP  (the  vari- 
able containing  the  number  of  seasons)  =  4,  we  derive  the  follow- 
ing table  of  indices: 


TABLE  1. 


INDEXING  EXAMPLE 


STATIC  SOLUTION 

SEASON  INDEX 

RADIUS  INDEX 

INDEX  (ITM) 

(ITEMP) 

(IRAD) 

1 

1 

1 

2 

1 

•   2 

3 

1 

3 

4 

2 

1 

5 

2 

2 

6 

2 

3 

7 

3 

1 

8 

3 

2 

9 

3 

3 

10 

4 

1 

11 

4 

2 

12 

4 

3 

The  second  equation  above  is  read  as:   IRAD  is  equal  to  the 
quantity  (ITM-1)  modulo  NRAD  plus  1,  and  can  also  be  expressed, 
using  the  truncation  of  integer  arithmetic,  as  (letting  IM  repre- 
sent (ITM-1)) 

IRAD  =  IM  -  (IM/NRAD)*NRAD  +  1 

All  the  output  for  the  static  solutions  was  already  stored  by 
index  ITM  so  the  only  remaining  change  was  to  the  indices  on  sea- 
sonally varying  input  parameters  (SMLT,  GGU,  and  ALPH)  from  the 
index  ITM  to  the  new  season  index  ITEMP. 

In  the  repeated  load  portion  of  the  program,  changes  in  sub- 
routine RANDOM  were  minor;  several  unused  common  blocks  were 
deleted  (/ SPACE/, /DERV/ ,  AND  /NERV/) ,  common  block  /DIST/  was 
added,  and  in  the  check  portion  mentioned  earlier  (in  the  first 
modification  set)  the  season  and  radius  indices  are  derived  from 
the  index  on  static  solutions  and  printed  if  any  solution  coeffi- 
cient set  fails  the  tests. 
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To  discuss  changes  in  CRACKS  and  ROUGH,  it  is  necessary  first 
to  describe  the  modifications  to  subroutine  GRESP.   GRESP  as 
originally  written  performed  the  following  calculations  for  each 
season:   (1)  evaluate  the  time -temperature  shift  to  derive  an 
effective  load  duration  from  the  input  duration;  (2)  evaluate 
the  Dirichlet  series  for  the  required  normalized  response  at  that 
effective  load  duration;  (3)  evaluate  several  related  stochastic 
quantities;  (A)  evaluate  the  true  (unnormalized)  response  and  its 
variance  using  AMPLITUD  (tire  pressure)  and  its  variance,  the 
variance  of  load  duration,  and  the  coefficient  of  variation  of  the 
normalized  response.   Another  term  in  the  variance  of  the  true 
response  involves  the  covariance  of  load  duration  and  amplitude, 
but  this  covariance  had  been  forced  equal  to  0  within  the  program 
as  received,  and  has  not  been  changed.   With  the  addition  of  vari- 
able radii  and  corresponding  static  solutions,  and  of  variable 
amplitude,  it  appeared  that  the  two  portions  of  GRESP  could  be 
split  into  two  routines;  one  (GRESP1)  would  evaluate  the  normalized 
strain  and  its  variance  for  all  the  static  solutions;  the  other 
(GRESP2)  would  be  called  in  a  loop  in  the  calling  program  which 
isolates  the  correct  static  solution,  season,  and  radius  by 
separate  indices,  and  provides  GRESP2  the  proper  value  of  ampli- 
tude and  of  static  solution  coefficient  of  variation.   By  splitting 
the  calculations  in  this  fashion,  certain  economies  of  calculation 
and  of  bookkeeping  can  be  made,  e.g.,  it  is  easier  to  avoid  redundant 
or  repetitive  calculation  of  the  same  quantities.   The  rewrite  of 
GRESP  also  permitted  a  change  in  the  interests  of  computation  effi- 
ciency from  many  divisions  by  the  same  quantity  to  many  multipli- 
cations by  the  inverse  of  that  quantity.   (The  variable  Y  in  GRESP 
was  replaced  by  the  variable  Q  =  l./Y  in  GRESP1) . 

Another  lesser  change  was  made  in  the  code  which  was  incor- 
porated in  GRESP2.   A  number  of  variables  used  in  the  calculation 
of  the  response  and  its  variance  were  then  normalized  by  the  final 
values  of  response  and  variance  and  returned  to  the  calling  pro- 
gram; however,  they  were  never  used,  and  only  served  needlessly  to 
lengthen  the  calling  sequence.   They  were  therefore  removed  from 
the  calling  sequence  of  GRESP2  and  left  unnormalized. 

Changes  in  subroutine  ROUGH  included:   (1)  an  expansion  of  the 
variable  list  in  the  common  block/RESPl/,  which  previously  had  had 
most  of  its  contents  contained  in  one  dummy  variable  name  (this 
permitted  access  to  the  variable  AMP,  which  was  set  to  values  of 
amplitude  AMPL  corresponding  to  the  current  value  of  the  radius 
index);  (2)  the  addition  of  block/DIST/;  (3)  the  call  to  GRESPl; 
(4)  the  loops  around  the  call  to  GRESP2  and  (5)  the  definition  of 
season  radius,  and  static  solution  indices  where  appropriate.   Fol- 
lowing these  calculations  a  table  is  printed  listing  for  each  static 
solution,  the  temperature,  radius,  total  deflection  (elastic  plus 
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permanent),  variance  of  deflection,  and  ag  „  and  ysys.   The  total 
number  of  seasons  over  which  damage  is  accumulated  is  determined 
by  the  last  value  of  TRANDOM  input  (program  variable  name  TIME); 
this  is  limited  by  the  dimensions  on  variables  to  480.   A  loop  for 
the  damage  estimation  at  each  of  these  seasons  is  then  initiated, 
and  within  this  is  a  loop  over  the  number  of  radii;  thus  the  static 
solution  for  each  radius  is  exercised  in  order  during  each  season. 

The  traffic  appropriate  at  each  season  is  determined  from  the 
array  ADT(I),  the  average  daily  traffic  over  the  period  TIME(I-l) 
to  TIME(I)  (input  under  the  LAMBDA  key  word).   This  is  done  by 
determining  which  is  the  first  value  of  TIME  greater  than  the  mid- 
time  of  the  current  season,  and  using  the  corresponding  value  of 
ADT .   For  each  season,  the  season  index  (for  position  within  the 
year,  as  opposed  to  position  with  the  entire  simulation  period)  is 
defined  by  the  modulo  function,  and  the  traffic  is  expressed  in 
total  axles  per  second  by  multiplying  the  ADT  value  by  the  inverse 
of  number  of  seconds  per  day  and  by  the  seasonal  traffic  multiplier, 

Within  each  season  in  the  simulation  period,  a  loop  over  the 
number  of  radii  is  executed.   The  season  index  has  been  calculated 
and  the  radius  index  is  equal  to  the  inner  loop  index,  hence  the 
static  solution  index  can  be  immediately  obtained  and  used  to 
retrieve  the  appropriate  values  of  ALPHA,  GNU,  response  (vertical 
deformation)  and  variance  of  response.   The  number  of  loads  per 
second  calculated  previously  is  now  multiplied  by  the  fraction  of 
traffic  assigned  to  this  radius-amplitude  (or  wheel  load)  combina- 
tion. 

The  damage  (rut  depth)  to  date  is  converted  to  the  time,  in 
seconds,  required  to  generate  this  damage  under  the  current  condi- 
tions of  environment  and  traffic;  this  time  is  added  to  the  number 
of  seconds  per  season  and  the  sum  is  then  used,  in  conjunction 
with  the  current  traffic,  to  generate  the  equivalent  total  number 
of  axles  under  current  conditions  of  traffic,  environment  and  wheel 
load  to  the  end  of  the  season.   This  number  of  vehicles  is  then 
entered  in  the  equation  for  permanent  deformation  along  with  the 
current  ALPHA  and  GNU  to  obtain,  after  all  the  solutions  for  dif- 
ferent radius  have  been  exercised,  the  rut  depth  and  its  variance 
at  the  end  of  the  current  season.   The  process  is  repeated  to  the 
end  of  the  simulation;  the  result  is  a  table  of  rut  depth  and  vari- 
ance of  rut  depth  for  the  end  of  each  season  in  the  simulation 
period.   Values  of  rut  depth  and  its  variance  are  then  interpolated 
for  the  specific  time  values  in  array  TIME. 

As  before,  GVNRM,  a  multiplicative  factor  in  the  conversion 
of  variance  of  rut  depth  to  slope  variance,  is  defined  to  be  the 
unweighted  rms  average  of  the  values  of  the  coefficients  of 
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variation  of  the  Dirichlet  series  for  vertical  deformation  for  all 
of  the  static  solutions.   Some  thought  might  be  given  towards 
weighting  the  values  for  a  given  radius  by  the  traffic  fraction 
corresponding  to  that  radius,  or  perhaps  by  the  average  (over  sea- 
son) of  the  effect  (vertical  deformation)  for  that  wheel  load. 
This  was  not  done  here;  it  was  observed  that  the  values  did  not 
vary  much  with  radius. 

Subroutine  CRACKS  was  modified  in  ways  similar  to  those  used 
for  ROUGH.   However,  there  exist  some  features  unique  to  each 
routine,  and  these  will  be  identified  and  discussed.   CRACKS  pro- 
vides an  algorithm  for  estimation  of  the  two  parameters  K-i  ,  K2  of 
the  usual  failure  equation 

N.  ..  =  K.  eK2 
fail    1 

as  functions  of  temperature  and  of  the  values  for  K  and  K2  at  70  °F 
(21  °C)  in  the  case  that  values  are  not  available  for  each  of  the 
seasonal  temperatures.   This  algorithm  and  related  statements  were 
moved  to  the  beginning  of  the  routine,  out  of  the  original  loop 
over  the  season  index.   When  values  either  input  or  calculated  for 
Ki  and  K2  are  available  for  each  season,  GRESP1  is  called;  then  a 
loop  over  season  index  and  an  inner  loop  over  radius  index  are 
initiated,  the  appropriate  amplitude  is  recalled  from  the  array  of 
amplitude  values,  and  GRESP2  is  called,  returning  the  radial  strain 
under  the  asphalt  layer  and  its  variance.   These  values,  along  with 
K,  and  K„  and  related  stochastic  parameters  for  the  current  season, 
are  entered  into  a  series  of  rather  long  and  complex  expressions 
to  obtain  Nfa:Q  and  its  variance  for  each  season  and  wheel  load 
(indexed  by  the  static  solution  index). 

The  values  of  strain,  variance  of  strain,  K^,  K2 ,  ^fa±±   and 
its  variance  are  printed  for  each  static  solution,  along  with  its 
identifying  temperature  and  radius.   Then,  as  in  ROUGH,  a  loop  is 
initiated  over  the  total  number  of  seasons  required  to  reach  the 
last  value  input  under  the  TRANDOM  keyword,  the  season  index  is 
generated  for  each  season,  and  within  that  season  a  loop  over 
radius  is  performed  with  its  accompanying  calculation  of  the  static 
solution  index.   The  appropriate  average  daily  traffic  is  selected 
just  as  before,  but  then  is  converted  to  total  axle  (loadings)  in 
the  season,  rather  than  to  axles  per  second.   This  is  converted  to 
axles  per  season  for  each  radius,  (Nj^ k )  and  Miner's  Law  is  then  in- 
voked to  obtain  the  damage  index  for  each  season. 


radius     N 


DI (season,  )  =  DI  (season,   )  +     £ 


i,  k 


N 

i=l      fail 


i,k 
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The  variance  of  damage  index  is  calculated  similarly. 

Since  cracking  damage  is  a  strictly  accumulative  process  (at 
least  as  modeled  here)  it  is  possible  to  separately  accumulate  the 
damage  and  its  variance  attributable  to  each  radius  (wheel  load) . 
This  has  been  done,  and  the  variables  DSUM  and  VSUMS  indexed  by  the 
radius  index,  contain  this  information  at  the  end  of  the  simula- 
tion.  These  are  passed  to  a  new  routine,  CRPRNT,  where  they  are 
normalized  by  the  total  damage  and  total  variance  for  all  radii 
and  are.  printed  in  a  table  with  radius,  amplitude,  resulting  wheel 
load  (the  only  place  wheel  load  is  printed)  and  fraction  of  traffic 
at  that  wheel  load.   This  allows  the  effect  of  higher  axle  loadings 
on  the  pavement  performance  in  to  be  quickly  ascertained  in  terms  of 
the  ratio  of  fractional  damage  to  fractional  number  of  axle  loads. 

Low  Temperature  Cracking  Predictions 

The  available  models  for  prediction  of  low-temperature  crack- 
ing include  Program  COLD  (Ref.  5  ),  the  Shahin-McCullough  Model 
(Refs.  62  and  7-*)  ,  and  the  Hajek-Haas  Model  (Refs.  84  and  9^)  .   Study 
of  these  models  indicated  the  following: 

1.  Program  COLD  is  aimed  at  predicting  if  cracking  will  occur, 
not  when  it  will  occur  and  at  what  crack  frequency. 

2.  The  Shahin-McCullough  Model  does  offer  cracking  frequencies 
and  output  oh  a  probabilistic  basis,  but  requires  quite  a 
bit  of  sophisticated  input  and  does  not  appear  to  offer 
much  more  reliable  predictions  than  the  much  simpler 
Hajek-Haas  Model. 
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3.   The  Hajek-Haas  Model  is  a  statistical  development  based 
on  measured  data  that  is  only  strictly  applicable  to  the 
areas  of  Canada  in  which  the  crack  observations  were 
taken.   However,  the  parameters  it  considers  are  those 
known  to  be  important  to  the  formation  of  transverse 
cracks  and  it  should  have  reasonable  applicability  to 
areas  of  the  northern  United  States  where  transverse 
cracking  is  most  prevalent. 

Mr.  E.  C.  Novak  reported  on  a  detailed  study  in  Michigan 
(Ref .  10  )  and  found  the  Hajek-Haas  equation  to  be  unreliable  for 
Michigan  pavements.   Unfortunately,  the  equation  published  in 
Reference  9  was  incorrect,  so  his  study  was  duplicated  during  this 
project  using  the  correct  equation  and  the  results  appear  in 
Appendix  D. 

It  appears  that  the  Hajek-Haas  Model  offers  almost  the 
reliability  of  the  more  rational  Shahin-McCullough  Model  and  has 
the  important  added  advantage  of  relative  simplicity  in  providing 
input  values  for  its  use.   Both  procedures  are  limited  by  the 
high  sensitivity  of  low-temperature  cracking  (and  indeed  any  frac- 
ture problem  involving  tensile  stresses)  to  many  more  parameters 
than  can  be  reasonably  or  accurately  defined  in  the  analysis,  and 
neither  will  consistently  give  reliable  predictions.   It  may  be 
noted  from  field  data  such  as  that  in  Table  34  of  Appendix  D  that 
there  are  wide  variations  in  actual  cracking  observed  for  apparently 
similar  pavement  sections,  so  it  is  not  surprising  that  predictive 
models  should  be  unable  to  attain  a  high  confidence  level  when  so 
much  variability  exists  in  nature. 

In  view  of  the  results  noted  in  Appendix  D,  the  considera- 
tions discussed  above  and  the  results  of  a  factorial  study  of  com- 
puted results  for  a  broad  range  of  input  values,  the  Hajek-Haas 
Model  was  selected  for  use  and  has  been  incorporated  into  VESYS  A. 

The  Hajek-Haas  Model  consists  of  a  single  rather  complex 
regression  equation  which  yields  Cracking  index  I,  the  number  of 
full  plus  half  the  number  of  half  transverse  cracks  per  500  feet 
(152  m)  of  roadway,  from  the  values  of:   m,  the  winter  design  tempera- 
ture, °C;  s,  the  stiffness  of  the  asphalt  cement  at  temperature  m, 
kg/cm  ,  as  defined  by  McLeod's  method;  a,  the  age  of  the  pavement, 
years;  t,  the  thickness  of  the  asphalt,  inches;  and  d,  a  dimensionless 
code  for  the  type  of  subgrade  (5-sand,  3-loam,  2-clay) .   The  correct 
equation  is: 


1 
Novak,  E.  C,  Jr.,  "Evaluation  of  a  Model  for  Predicting 

Transverse  Cracking  for  Flexible  Pavements,"  Research  Division 

Michigan  Highway  Commission,  February  1976. 
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101  =  2.497  x  1030  x  (O.ls)  (6-7966  "  -8740t  +  1"3388a) 

x  (7.054  x  10-V  x  (3.193  x  uf13)-0^  x  d  ("06026s) 

or,  in  a  more  usable  form  (taking  the  logarithms  of  both  sides), 
I  =  30.3974  +  (6.7966  -  .8704t  +  1.3388a)  (log1Qs  -  1) 
-  2.1516d  +  1.24958m  +  (,06026s)  (log  d) 

Subroutine  LOTMPC  was  written  to  evaluate  this  equation  from  the 
asphalt  thickness  THICK1,  the  times  TRANDOM,  and  the  values  of  stiff- 
ness, winter  design  temperature,  and  subgrade  code  input  under  the  new 
keywords  STIFF,  WINTTEMP,  and  SUBTYPE.   These  keywords  were  added  to 
the  appropriate  data  statement  and  suitable  default,  maximum  and  mini- 
mum values  were  provided  (See  Table  3  at  the  end  of  this  chapter). 
These  new  parameters  are  made  available  to  LOTMPCR  and  the  new  common 
block  /LOTCR/,  which  also  contains  the  computed  values  of  cracking  in- 
dex for  all  the  time  values  input  under  TRANDOM.   These  time  values 
are  converted  to  seconds  for  use  in  the  cracking  and  roughness  calcula- 
tions, but  a  loop  in  CRACKS  converts  them  to  years  before  LOTMPC  is 
called;  they  are  then  converted  to  seconds  on  return  from  LOTMPC. 

(A  note  of  warning:   The  variable  YEAR  which  is  used  for  the 
conversion  is  perhaps  incorrectly  named;  it  actually  contains  the 
constant  required  to  convert  TRANDOM  time  values  to  seconds  from 
whatever  units  they  are  input  in,  specified  by  the  TUNITS  keyword. 
If  low  temperature  cracking  is  to  be  calculated,  and/or  if  the 
column  headings  YEARS  in  several  of  the  repeated  load  subroutines 
are  to  be  accurate,  the  times  input  under  the  TRANDOM  keyword  must 
be  in  years,  and  TUNITS  must  be  6.) 

Cracking  index  values  returned  from  LOTMPC  are  converted  to 
units  of  yd2  cracked  per  1000  yd2(0.8  m2  per  800  m2)as  follows: 

Let  w  be  the  width  of  the  roadway  (ft), 

and  n  be  the  number  of  full  cracks  +  half  the  number  of 

half  cracks  in  a  500  foot  (152  m)  length  of  roadway 

(cracks  assumed  to  be  transverse) 

then  —5 —  is  the  area  in  square  yards  of  500  feet  (152  m)  of 
nw  is  the  total  crack  length  in  feet. 

nw     9nw    9n    18n 


thus 


500w/9  ""  500w   500   1000 
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is  the  total  crack  length  (ft)  per  square  yard  of  pavement,  aiul 
we  have  then 

18  n  ft  (1  yd) 

1000  sq.  yd.  x  (100°  sq'  yds'}  X  TTftT  =  6n 

yards  of  crack  length  per  1000  square  yards  of  pavement. 

If  we  now  assume  a  single  transverse  crack  6n  yards  long,  it 
will  pass  over  6n  different  "square  yards";  we  can  also  imagine 
transverse  cracks,  well  spaced,  whose  total  length  is  6n  yards,  and 
the  total  number  of  square  yards  transverse  will  remain  6n .   This 
will  remain  true  until  the  average  crack  spacing  becomes  small 
enough  that  the  random  placement  of  cracks  leaves  more  than  one 
crack  in  a  "square  yard."   The  smallest  average  crack  spacing  ob- 
served in  this  investigation  was  about  5  yards  (4.5  m) ;  therefore 
the  derivation  and  result  are  assumed  valid  here,  to  the  required 
degree  of  accuracy. 

Finally,  the  values  of  cracking  index,  average  crack  spacing, 
and  square  yards  cracked  per  1000  square  yards  for  low  temperature 
cracking  are  printed  with  the  results  for  fatigue  cracking;  an  addi- 
tional column  contains  the  total  square  yards  cracked  per  1000 
square  yards,  approximated  here  by  the  simple  sum  of  the  results 
from  the  two  distress  mechanisms.   This  tends  to  maximize  the 
effect  of  low  temperature  cracking,  since  the  two  mechanisms  yield 
cracks  in  perpendicular  directions,  and  a  single  square  yard  could 
easily  have  both  types  of  cracks.   This  total  is  limited  to  a  maxi- 
mum value  of  1000  before  being  printed. 

It  should  be  noted  that  in  the  derivation  of  the  regression 
equation  the  winter  design  temperature  was  defined  as  lowest  tem- 
perature at  or  below  which  only  1  percent  of  the  hourly  air  tem- 
peratures in  January  occur  for  the  severest  winter  during  a  10 
year  period.   Hajek  and  Haas  have  provided  a  relation  between  the 
design  temperature  defined  this  way  and  freezing  index,  and  also 
provide  freezing  index  maps  for  the  United  States  and  Canada. 
These  aids,  which  will  be  discussed  further  in  Chapter  IV  make  this 
method  useful  even  where  (as  is  probably  true  in  most  places)  such 
detailed  weather  records  are  not  available. 

Subroutine  to  Save  Static  Solutions 

In  the  discussion  of  program  modifications  above,  reference  was 
made  to  an  alternative  capability  to  that  for  isolated  Type  3  runs 
(those  not  following  Type  1  runs).   This  was  done  for  reasons  men- 
tioned there,  but  in  addition  the  isolated  Type  3  runs  (involving  use 
of  punched  computer  files)  are  not  pracitical  at  a  remote  job  entry 
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site  with  no  card  punch  facility.   However,  it  is  useful  and  con- 
venient to  be  able  to  rerun  simulations  in  which  only  parameters 
affecting  repeated  load  solutions  are  changed;  therefore,  a  new 
means  of  saving  the  results  of  the  relatively  more  time-consuming 
static  solution  was  designed  and  implemented  in  VESYS  A. 

The  keywords  SAVETAPE,  READTAPE  and  NOTAPE  control  the  genera- 
tion and  use  of  a  file  (TAPE  1;  or  on  logical  unit  1)  holding  the 
static  solution  data.   If  SAVETAPE  is  read  at  the  beginning  of  a 
type  1  run,  then  at  the  end  of  the  computations  for  static  solu- 
tions the  relevant  data  (detailed  in  Appendix  E)  are  written  to 
logical  unit  1,  followed  by  an  end  of  file.   If  this  is  the  first 
such  set  of  solutions  in  a  job,  unit  1  is  rewound  before  the  file 
is  written.   Physically  unit  1  can  be  either  a  tape  file  or  a  disk 
file,  but  one  must  be  careful  in  the  latter  case  to  declare  that 
file  to  be  a  permanent  file  (in  CDC  terminology)  or  otherwise  to 
ensure  that  it  will  remain  accessible  at  a  later  time.   In  the 
case  of  magnetic  tape  files,  the  initial  rewind  may  be  undesirable, 
as  positioning  after  other  files  already  on  the  tape  might  be 
wanted.   In  such  a  situation  the  statements 

IF(IENT.EQ.O)  REWIND  1 
IENT  =  1 

should  be  removed  from  subroutine  DOTAPE.   After  writing  to 
logical  unit  1,  the  program  then  writes  the  same  data  on  the 
output  file  on  a  special  page. 

The  keyword  READTAPE,  when  read  for  a  Type  3  run,  instructs 
the  program  to  read  data  from  logical  unit  1  into  special  arrays 
in  subroutine  DOTAPE,  after  all  other  card  input  data  for  the  run 
has  been  read.   Validity  checks  are  then  made  between  the  card 
input  data  and  the  data  read  from  the  external  file  as  follows: 

1.  Is  the  number  of  seasons  (keyword  NTEMPS)  the  same? 

2.  Is  the  number  of  radii  (NRADIUS)  the  same? 

3.  Is  the  number  of  Dirichlet  coefficients  (NDELTAS)  the 

same? 

4.  Is  the  thickness  of  asphalt  (THICK1)  the  same,  to  ±  0.2 

inches  (0.5  cm)? 

5.  Is  the  thickness  of  the  base  (THICK2)  the  same,  to  ±  0.2 

inches  (0.5  cm)  ? 

If  all  these  questions  are  answered  affirmatively,  the  data  is 
then  stored  in  the  normal  manner  for  use  in  the  repeated  load 
solution  just  as  if  it  had  been  calculated  on  the  same  computer  run. 
If  one  or  more  of  these  tests  are  failed,  then  the  parameters 
tested,  both  those  input  from  cards  and  those  input  from  the  data 
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file,  are  printed  and  the  program  terminates  its  execution.   In 
all  cases,  the  contents  of  the  title  array  from  the  data  file  are 
printed  along  with  the  above  parameters  and  the  values  of  tempera- 
ture and  radius. 

The  keyword  NOTAPE  is  used  in  the  second  and  following  data 
decks  (ended  by  ENDOFRUN,  not  by  ENDOFJOB,  which  terminates  execu- 
tion) of  a  multiple  Type  3  computer  run,  where  all  of  the  decks 
are  for  runs  making  use  of  the  same  set  of  static  solutions.   This 
causes  the  program  to  bypass  the  reading  of  the  file,  since  the 
data  are  already  in  the  memory.   It  is  required,  not  optional,  in 
such  cases.   Where  more  than  one  file  of  data  exists  on  logical 
unit  1,  which  is  never  rewound  in  the  read  (Type  3)  mode,  such 
data  could  be  read  and  used  either  by  omitting  the  NOTAPE  directive 
(in  the  case  of  one  simulation  run  per  set  of  static  solutions)  or 
by  entering  the  READTAPE  key  word  again  in  the  first  deck  making  use 
of  the  new  set  of  static  solutions. 

Appendix  E  describes  the  structure  of  these  static-solution 
files. 

Verification  of  Program  Changes 

The  many  changes  to  VESYS  IIM  to  produce  VESYS  A  required  that 
a  program  verification  analysis  be  performed:   the  changes  of 
indexing,  dimensions,  common  blocks  and  input,  even  more  than  those 
involved  in  direct  computation,  were  so  numerous  that  small  errors 
could  easily  creep  in.   This  verification  was  performed  after  each 
major  set  of  program  changes,  and  consisted  of  tests  against  pre- 
vious runs  and  comparisons  with  what  might  reasonably  be  expected 
(in  the  testing  of  new  features). 

The  first  set  of  verification  runs  was  made  after  the  modifi- 
cation for  seasonal  variation  of  materials  properties.   They  were 
made  with  a  basic  data  deck  corresponding  to  runs  made  at  ARE  in 
1975  to  test  the  expected  life/service  life  routines;  the  deck  was 
modified  as  necessary  for  the  new  features.   The  term  OLD  program 
will  refer  to  VESYS  IIM  as  received  from  FHWA  in  June  1976,  with 
the  variables  LY1  and  LY2,  used  in  the  UNLOAD  portion  of  the  static 
solutions,  returned  to  their  previous  (ARE)  values  of  6  and  8, 
respectively. 

The  sequence  of  runs  was  as  follows  (runs  are  identified  by 
the  last  two  characters  of  the  CDC  job  name): 
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RUN     DATE     PROGRAM     COMMENTS 

DM      8/30       OLD       Output  identical  to  service  life 

examples,  FH-1,  1975. 

4E      9/23       OLD       T  =  60°F  (15.6°C),  1  season. 

N3      9/30       NEW       T  =  60°F  (15.6°C),  1  season. 

Output  identical  to  4E. 

H7      10/6       NEW       T  =  60°F  (15.6°C),  4  season. 

Output  identical  to  4E  and  N3 
except  for  variance  of  cracking. 

Jobs  DM,  4E,  N3,  and  H7  demonstrate  that  for  NTEMP  =  1  the 
program  options  selected  therein  (TYPE  1  and  TYPE  3  runs)  were  not 
affected  by  the  changes,  and  further,  that,  with  the  exception  of 
the  variance  of  cracking  (a  separate  problem  in  formulation  dis- 
cussed in  Appendix  F) ,  the  results  from  the  use  of  four  identical 
"seasons"  (involving  four  static  load  solutions  and  the  correspond- 
ing complex  indexing  of  data  and  results)  are  the  same  as  those 
from  a  run  using  only  one  "season."  The  exception  mentioned  is  not 
a  result  of  the  program  changes,  as  the  particular  algorithm  involved 
was  not  modified;  in  fact,  a  later  demonstration  run  on  the  version  as 
received  from  FHWA  showed  the  same  changes  in  variance  of  cracking 
with  number  of  seasons.   The  source  of  the  changes  was  later 
demonstrated  mathematically  to  be  inherent  in  the  algorithm. 

The  second  verification  sequence  followed  the  addition  of 
multiple-radius  solutions;  here  there  were  even  more  changes  of 
indexing  and  dimensions  (commonly  lumped  under  the  term  "book- 
keeping") ;  in  addition  a  major  conceptual  change  was  made  in  the 
manner  in  which  a  distribution  of  axle  (or  wheel)  loads  was  handled. 
Again,  an  extensive  comparison  was  performed  with  the  aim  of 
exercising  all  the  new  options  and  all  of  the  new  indexing  and 
associated  loops. 

The  following  computer  runs  were  made  with  the  data  input  as 
indicated  (the  basic  data  deck  was  the  same  as  that  used  for  the 
previous  verification  runs;  comparison  is  with  job  N3  in  that  veri- 
fication sequence,  and  hence  to  the  results  from  the  original  pro- 
gram with  which  job  N3  agreed  in  every  detail): 

Name     NTEMP     NRADIUS     Comments 

B4         1  1        Same  as  N3  except  for  the  variance 

of  the  normal  deflection  and  the 
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Name     NTEMP     NRADIUS     Comments 

resulting  roughness  parameters  (very 
slight  difference) .   All  cracking 
output  and  the  rut  depth  itself  agree. 
(VCAMP  =  2644.) 

DJ         2  2        4  static  solutions,  all  for  T  =  60°F 

(15.6°C),  R  =  4.24  inches  (10.8  cm). 
Results-4  solutions  are  all  identical, 
and  agree  with  single  solution  in  B4. 
Repeated  load  gives  same  rut  depth  as 
B4.   VCAMP  -  0  so  variances  are  different 
(traffic  distribution  was  75%-25%  on  the 
two  (equal)  radii,  each  using  75  psi 
(518  kN/m^)   tire  pressure). 

UB         1  1        VCAMP  =  0.   Everything  in  repeated 

load  solution  is  identical  to  that  in 
DJ  except  the  variance  of  cracking; 
the  value  obtained  can  be  mathe- 
matically predicted  from  that  for  DJ, 
and  the  predicted  and  observed  values 
agree. 

60        1         1       VCAMP  =  2644.   Same  input  as  B4; 

run  after  locating  and  correcting  an 
error  in  subroutine  ROUGH.   Now 
gives  answers  identical  to  N3  in  all 
respects  (used  static  solutions 
generated  by  B4). 

J4        1         4       All  radii  identical-test  of  indexing 

changes.   This  job  and  the  one  fol- 
lowing (NB)  used  a  different  data 
set  from  the  others  in  this  sequence. 
Lambda  =  5000  (constant  with  time). 

NB         1  4        Same  as  J4  (used  same  static  solu- 

tions) except  LAMBDA  =  8000.   Test 
of  concept  of  saving  and  reusing 
static  solutions,  and  of  changes  in 
repeated  load  subroutines.   The  ratio 
of  "value"  (J4)  to  "value"  (NB) , 
where  "value"  is  a  predicted  behavior, 
was  in  all  cases  what  one  expected 
from  the  equations. 
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In  addition  to  this  series  of  comparisons ,  computer  runs  have 
been  made  modeling  several  AASHO  test  sections  and  are  discussed 
in  Chapter  V.   Suffice  it  here  to  say  that  using  4  seasons  and  4 
radii  for  a  total  of  16  static  solutions  in  each  case,  the  program 
generally  behaved  well  and  gave  answers  encouragingly  close  to 
those  observed. 

One  point,  however,  deserves  specific  notice  relative  to  veri- 
fication of  program  changes.   In  the  assumption  of  a  traffic  dis- 
tribution (modeling  the  graph  of  number  of  axles  vs.  axle  load  by 
a  4  bar  histogram)  the  actual  random  time  distribution  of  loads  is 
replaced  by  a  strict  sequence:   within  each  season  all  the  loads  of 
the  first  weight  are  applied,  then  all  the  loads  of  the  second 
weight,  etc.   It  was  felt  that  this  was  indeed  more  accurate  than 
representing  all  the  loads  by  a  single  radius,  a  single  amplitude, 
and  a  large  value  of  the  variance  VCAMP;  however,  if  the  order  in 
which  the  loads  were  applied  made  a  significant  change  in  the 
results,  there  would  be  a  problem  in  interpreting  the  results. 

For  fatigue  cracking,  the  cumulative  damage  (as  modeled 
in  VESYS  A)  is  a  completely  additive  process  (Miner's  law),  and 
hence  the  order  of  the  load  is  immaterial.   For  rut  depth,  the 
process  involved  is  quite  different:   at  the  beginning  of  a  new 
season/load  combination,  the  damage  (rut  depth)  to  date  is  con- 
verted in  effect  to  a  number  of  load  applications  which  under  the 
new  conditions  (load,  temperature,  etc.)  would  yield  that  rut 
depth.   This  is  added  to  the  number  of  load  applications  for  the 
new  season/load  combination,  and  a  new  rut  depth  calculated.   It 
is  not  obvious  from  the  structure  of  the  equations  involved  that 
the  order  in  which  the  loads  are  applied  will  have  little  effect 
on  the  final  rut  depth,  although  after  several  years  (many  seasons) 
the  order  within  seasons  might  be  expected  to  have  less  effect. 

In  fact,  however,  hand  calculations  using  the  same  equations 
and  static  solutions,  but  reversing  the  order  of  the  loads, 
resulted  in  very  small  changes  in  the  final  rut  depth.   Speci- 
fically, for  AASHO  section  266,  the  program  predicted  ".092009 
rut  depth  after  one  year.   A  hand  calculation  in  the  same  order 
differed  by  .03%;  a  hand  calculation  in  which  the  order  of  loads 
was  reversed  (to  heaviest  first)  differed  by  0.5%,  and  a  hand  cal- 
culation (rather  extreme)  in  which  the  order  of  seasons  was 
reversed  (to  Fall,  Summer,  Spring,  Winter)  and  the  loads  were 
applied  as  in  the  program  (lightest  to  heaviest)  differed  by  2.2% 
after  one  year.   The  first  year  is  the  one  in  which  the  rut  depth 
increases  most  rapidly.   In  addition,  hand  calculations  for  the 
second  year,  for  lightest  loads  first  and  for  heaviest  loads  first, 
yielded  rut  depths  differing  only  by  0.2%;  apparently  the  error 
becomes  less  as  time  goes  on,  as  was  expected. 
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The  last  step  in  verification  involved  the  input,  calculation 
and  printout  of  the  low  temperature  cracking;  the  values  printed 
were  compared  with  values  obtained  earlier  by  hand  and  from  a 
large  factorial  study  using  only  the  low-temperature  cracking 
routine  and  a  simple  driver  program.   The  results  were  the  same, 
indicating  again  that  the  new  input  parameters  were  being  properly 
indexed,  stored  and  used. 

Later  extensive  use  of  the  program  has  increased  our  con- 
fidence in  the  new  version;  we  have  encountered  nothing  in  six 
months  of  nearly  continuous  use  to  suggest  the  presence  of  errors 
in  the  modification. 

Summary  of  Keywords  Added,  Deleted  and  Changed 

In  the  process  of  modifying  the  VESYS  IIM  program  to  provide 
the  new  capabilities  described  above,  new  forms  of  data  were 
required,  and  the  need  for  some  old  forms  was  removed.   These 
changes  are  reflected  in  the  modifications  to  the  key  word  "dic- 
tionary," all  of  which  are  described  in  the  internal  user's  guide; 
they  will  be  listed  and  briefly  discussed  here  also. 

Table  2  lists  the  keywords  which  were  removed.   They  were 
associated  with  card  input  for  isolated  Type  3  runs,  providing  the 
data  summarizing  an  earlier,  single,  static  solution.   The  asso- 
ciated variables,  with  the  exception  of  those  associated  with 
GNURAD  and  ALFRAD  (not  used  previously)  are  still  used  in  the  same 
manner  as  before;  only  the  capability  of  card  input  for  these 
variables  has  been  removed. 


TABLE  2 

Keyword 

CREEP1 
CREEP2 
CREEP3 

GRADIAL 
GNORMAL 

ALFRAD 
ALFNORM 
GNURAD 
GNUNORM 

COEFGRAD 
COEFGNRM 


LIST  OF  KEYWORDS  REMOVED  FROM  VESYS-A 

Associated  Variable  (original  program) 

E  (1,1) ,1=1,  NDEL 
E  (1,2)  ,  I  =  1,  NDEL 
E  (1,3),  1=1,  NDEL 

GRAD(I)  ,1  =  1,  NDEL 
GNORM(I),  1=1,  NDEL 

ALF(7) 
ALF(4) 
GMU(7) 
GMU(4) 

GCVAR(7) 
GCVAR(4) 
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Table  3  lists  the  new  keywords,  the  associated  internal  vari- 
able names  and  dimensions,  the  default  values,  and  the  upper  and 
lower  limits  if  a  single  variable  and  not  an  array.   Some  of  the 
limits  have  not  been  changed  to  appropriate  values  for  keywords 
which  replaced  others  in  the  keyword  data  statement;  in  general 
these  are  cases  where  true  limits  could  be  applied  where  previously 
the  limits  were  very  nearly  ±  °°. 

It  should  be  noted  (as  it  is  in  the  internal  input  guide)  that 
the  keywords  SEASMULT,  GNU,  and  ALPHA  are  accompanied  by  an 
integral  value  expressed  as  a  decimal  digit  that  indicates  the 
layer  to  which  the  data  read  from  the  following  card  apply;  hence 
one  will  normally  have  several  occurrences  of  each  of  these  key- 
words in  a  data  deck. 

Another  point  to  note  is  that  if  NRADIUS  is  either  set  equal 
to  1  or  not  specified  (defaulting  to  1)  then  the  RADIUS  and 
AMPLITUD  keywords  are  interpreted  as  they  were  in  VESYS  IIM:   the 
value  in  the  numeric  field  immediately  following  the  keyword  con- 
tains the  single  value  of  radius  or  amplitude  required,  and  no 
additional  data-only  card  is  read.   Thus,  if  NRADIUS  =  1,  using 
the  default  value  for  TRDIST(l)  =1.,  the  program  behaves  exactly 
(with  respect  to  radius  and  tire  pressure  specification)  as  VESYS 
IIM  did. 
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CHAPTER  III 
SELECTION  OF  INTERSTATE  HIGHWAYS  FOR 
FOUR  ENVIRONMENTAL  ZONES 


The  environment  in  which  a  pavement  system  exists  has  a  signifi- 
cant effect  upon  that  pavement ?  s  response  to  the  traffic  loads  imposed 
upon  it.   In  order  to  consider  this  important  parameter  in  this  study, 
FHWA  asked  that  four  interstate  highway  sections  be  selected  for  the 
study,  one  from  a  cold  climate  and  generally  high  humidity  and  rainfall 
to  be  called  "wet-freeze",  one  from  a  cold  climate  with  little  rainfall 
and  humidity  to  be  called  "dry-freeze",  one  from  a  relatively  warm  cli- 
mate with  high  humidity  and  rainfall  to  be  called  "wet-nb  freeze",  and 
one  from  a  relatively  warm  climate  with  little  rainfall  and  humidity  to 
be  called  "dry-no  freeze". 

As  there  are  numerous  highway  sections  that  would  meet  the  environ- 
mental requirements,  the  opportunity  existed  for  making  selections  based 
on  a  number  of  practical  considerations  including:   1)  Availability  of 
data  to  describe  the  pavement  structure,  the  traffic  and  axle  load  dis- 
tribution and  the  environment  from  the  pavement  viewpoint  for  the  VESYS  A 
computer  program,  2)  The  interest  level  of  State  Department  of  Transpor- 
tation personnel  in  VESYS  and  this  project  as  this  would  impact  upon  the 
cooperation  received  in  accumulating  the  required  data,  3)  Traffic  levels 
for  commercial  vehicles  so  that  two  levels  of  traffic  may  be  established, 
one  low  and  one  high  and  4)  Pavement  structure  so  that  two  levels  could 
be  used,  one  quite  substantial  and  one  relatively  thin  section. 

The  Departments  of  Transportation  at  both  Utah  and  Florida  were  en- 
gaged in  the  implementation  of  the  VESYS  IIM  system  and  had  active  programs 
of  their  own  to  accumulate  data  and  test  materials  in  order  to  service 
the  input  requirements.   This  involvement  offered  an  opportunity  to  obtain 
rather  directly  information  significant  to  the  successful  idealization 
of  pavement  systems  from  personnel  with  a  vested  interest  in  the  VESYS 
analytical  system,  so  it  was  decided  to  utilize  a  pavement  from  Utah  as 
dry-freeze  and  one  from  Florida  as  wet-no  freeze. 

Reference  ll-*-  also  classified  pavements  studied  according  to  the 
same  environmental  designations  and  gave  valuable  information  on  a  range 
of  flexible  pavements  that  could  be  considered  for  this  study.   Review  of 
this  reference  led  to  the  selection  of  a  section  of  IH-80  from  the  reha- 
bilitated AASHO  Test  Road  as  the  wet-freeze  section.   Other  considerations 


Darter,  M.  I.,  E.  J.  Barenberg  and  J.  S.  Sarvan,  "Maintenance-Free 
Life  of  Heavily  Trafficked  Flexible  Pavements,  a  paper  presented  at  the 
55th  Annual  Meeting  of  the  Transportation  Research  Board,  January  1976. 
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leading  to  this  selection  were  that  all  aspects  of  the  AASHO  Road 
Test  were  so  thoroughly  documented  and  that  permanent  deformation 
and  creep  compliance  test  results  were  already  available  for  soils 
in  this  area  from  previous  projects.   Data  was  available  for  this 
section  of  IH-80  in  References  12^  and  13  . 

This  left  the  selection  of  an  area  to  represent  the  dry-no 
freeze  condition,  which  was  reasonably  represented  by  sections  of 
California,  Arizona,  New  Mexico  or  Texas.   On  the  basis  of  availability 
of  data,  a  section  of  IH-20  between  Midland,  Texas  and  Odessa,  Texas 
was  selected  as  it  had  been  thoroughly  studied  in  other  projects  and 
considerable  data  was  available  about  its  materials,  its  traffic  and 
its  environment. 

Selections  of  specific  sections  of  Utah  and  Florida  were  made  in 
coordination  with  Department  of  Transportation  personnel  in  those 
states.   This  led  to  the  selection  of  a  section  of  IH-80N  near  the 
border  between  Utah  and  Idaho  as  dry-freeze  and  a  section  of  IH-10 
located  approximately  5  miles  (8  km)  west  of  Madison,  Florida  to 
represent  wet-no  freeze. 

The  interstate  highway  selections  were  the  result  of  research 
study,  discussion  with  State  Department  of  Transportation  personnel, 
coordination  with  the  Federal  Highway  Administration  and  discussions 
with  other  knowledgeable  persons.   The  final  selections  were  as  follows: 

1.  Wet-Freeze  -  a  section  of  the  rehabilitated  AASHO  Test  Road 
just  west  of  Ottawa,  Illinois,  that  appears  as  "Flex  7"  in 
Reference  5  and  Section  038  in  Reference  6.   This  section  is 

an  outside  lane  12  ft.  (3.6m)  wide  with  a  12  ft.  (3.6m)  shoulder, 
The  pavement  structure  consists  of  4.5  inches  (11.4  cm)  of 
asphalt  concrete,  8.5  inches  (21.6  cm)  of  crushed  stone, 
23.0  inches  (58.4  cm)  of  gravel  and  a  silty  clay  subgrade. 

2.  Dry-Freeze  -  a  section  of  IH-80N  just  south  of  the  Utah-Idaho 
border.   This  is  a  relatively  thin  section  consisting  of 
4.25  inches  (10.8  cm)  of  asphalt  concrete,  5.0  inches  (12.7 
cm)  of  gravel  base  and  a  subgrade  of  non-plastic  sandy  silt. 

3.  Wet-No  Freeze  -  a  section  of  IH-10  located  approximately  5 
miles  (8.05  km)  west  of  Madison,  Florida.   The  pavement  struc- 
ture for  this  section  consists  of  three  inches  of  asphalt  con- 
crete, 10.5  inches  (26.7  cm)  of  lime  rock  base,  12  inches 
(30.5  cm)  of  stabilized  sand  subgrade  (natural  sand  mixed  with 
25  percent  lime  rock  base)  and  a  natural  sand  embankment. 


1  Little,  R.  J.,  L.  J.  McKenzie  and  P.G.  Dierstein,  "The  Rehabilitated 
AASHO  Test  Road,  Part  I  -  Materials  and  Construction",  Interim  Report 
IHR-28,  Illinois  Department  of  Transportation,  July  1973. 

2 "The  AASHO  Road  Test  -  Report  5  -  Pavement  Research",  Highway 
Research  Board  Special  Report  61E,  1962. 
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Dry-No  Freeze  -  a  section  of  IH-20  in  Ecotr  County  between 
Midland,  Texas  and  Odessa,  Texas.   The  pavement  structure 
for  this  section  includes  3.3  inches  (8.4  cm)  of  asphalt 
concrete,  0.6  inches  (1.5  cm)  of  double  bituminous  surface 
treatment  19.0  inches  (48.3  cm)  of  flexible  base  and  sand 
subgrade. 
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CHAPTER  IV 

DEVELOPMENT  OF  INPUT  DATA  FOR  VESYS  A 
ANALYSIS 

The  data  requirements  for  VESYS  A  include  those  input  parameters 
used  in  VESYS  IIM  and  the  new  input  parameters  required  for  the  pro- 
gram improvements  described  in  Chapter  II.   It  was  necessary  to  draw  on 
a  number  of  sources  in  order  to  arrive  at  realistic  input  values  for 
the  solutions  required  by  this  study.   The  previous  work  reported  in 
Reference  1  provided  realistic  data  for  a  number  of  the  variables  of 
general  nature,  but  the  more  specific  information  required  to  define 
the  pavement  structure,  the  traffic  loads  imposed  and  environment 
had  to  be  obtained  from  other  publications  and  from  State  Department 
of  Transportation  personnel. 

The  pavement  sections  for  which  realistic  data  were  required  in- 
cluded: 

1.  Sections  266,  272,  308  and  470  in  the  AASHO  Road  Test. 

2.  Sections  3,  4,  11  and  17  from  the  Brampton  Test  Roads. 

3.  The  four  interstate  highway  sections  described  in  the 
previous  section. 

4.  The  factorial  of  8  real  and  hypothetical  sections  used 
in  the  factorial  study  of  cost  versus  legal  axle  limits. 

The  sections  identified  in  Items  1,  2  and  3  above  were  used  for 
"calibration"  of  the  model  while  those  described  in  Section  4  were 
utilized  for  the  study. 

The  input  requirements  introduced  by  the  improvements  to  VESYS 
IIM  embodied  in  VESYS  A  necessitated  the  establishment  of  seasonal 
pavement  temperatures,  seasonal  permanent  deformation  characterizations, 
seasonal  layer  stiffness  characterizations,  low  temperature  cracking 
input  values  and  axle  load  distributions  for  the  four  legal  axle 
limits.   The  development  of  the  input  data  for  these  five  general 
parameters  are  discussed  individually  below. 

Seasonal  Pavement  Temperatures 

It  is  very  important  to  establish  the  average  pavement  temperatures 
to  be  expected  for  the  particular  environment  in  which  the  pavement 
exists,  in  terms  of  environmental  data  readily  available.  Pavement 
temperature  has  been  found  to  be  a  strong  function  of  air  temperature 
and  of  solar  radiation  as  well  as  other  variables,  but  only  the  air 
temperatures  are  generally  readily  available.  In  order  to  estimate 
pavement  temperatures  for  the  diversity  of  environmental  zones  of 
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interest,  relationships  between  air  temperature  and  pavement  tempera- 
tures available  in  the  literature  were  studied  to  develop  a  means  of 
predicting  pavement  temperature  from  air  temperature  data  for  each 
environmental  zone.   The  data  considered  included  a  test  pavement  at 
College  Park,  Maryland  (Ref .  14)  ,  a  pavement  near  Potsdam,  New  York 
(Ref.  15)2,  a  pavement  in  Michigan  (Ref. 16) 3,  a  Texas  sandstone  mix 
(Ref.  17)4,  a  pavement  in  Florida  (Ref. 18) 5,  and  the  AASHO  Road  Test 
pavements  (Ref.  13).   The  approach  taken  was  to  study  whatever  data 
was  available  to  ascertain  how  the  pavement  temperature  related  to 
air  temperatures  on  a  daily  cycle,  monthly  cycle  and  seasonal  cycle 
to  the  extent  that  the  data  was  available.   Study  results  were 
accumulated  in  the  form  of  ratios  of  pavement  temperatures  to  air 
temperatures  and,  in  the  case  of  the  AASHO  Road  Test  pavements,  the 
difference  between  pavement  temperature  and  air  temperature.   The 
results  of  this  work  were  compared  to  similar  more  general  data 
developed  by  Witczak  (Ref.  19)".   Comparisons  of  the  Witczak  plotted 
curves  and  the  data  derived  from  the  study  indicated  that  Witczak' s 
curves  approximate  the  College  Park,  Maryland  data,  predicted  2  or 
3°F  (1.1  or  1.7°C)  high  for  Texas,  predicted  4  or  5°F  (2.2  or  2.8°C) 
low  for  Michigan,  and  predicted  as  much  as  8°F  (4.4°C)  low  for  New  York. 
The  general  conclusion,  based  on  this  data  would  then  be  that  the 
Witczak  plots  underestimate  pavement  temperatures  for  the  north,  esti- 
mate very  closely  pavement  temperatures  for  College  Park,  Maryland,  and 
slightly  overestimate  pavement  temperatures  in  the  south. 

The  VESYS  analytical  system  can  consider  temperatures  for  what- 
ever increments  of  time  are  assigned.   In  previous  work,  temperatures 
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2 
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minous Pavement  Temperature  Related  to  Climate",  Highway  Research 
Record  256,  1968. 
3 
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150,  1968. 
4 
Long,  R.E.,  "Relationships  Between  Control  Tests  for  Asphalt 

Stabilized  Materials",  Texas  A  &  M  University,  PhD  Dissertation, 
May  1971. 
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termine Low  Temperature  Cracking  Potential",  Report  to  Florida 
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Witczak,  M.W.,  "Design  of  Full-Depth  Airfield  Pavement", 
Proceedings,  Third  International  Conference  on  the  Structural 
Design  of  Asphalt  Pavements,  1972. 
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were  assigned  monthly,  utilizing  the  mean  temperatures  for  that 
month.   For  VESYS  A,  the  use  of  12  "seasons"  was  not  feasible 
because  VESYS  A  generates  a  static  solution  for  each  combination 
of  season  and  radius  in  the  discretized  axle  load  distribution. 
This  means  that  a  combination  of  4  seasons  and  4  radii  in  VESYS  A 
will  require  approximately  16  times  as  much  computational  effort  as 
a  similar  problem  for  VESYS  IIM  because  the  very  large  majority  of 
the  computational  effort  is  required  for  the  static  solutions.   For 
the  reasons  previously  discussed  on  pages  6  and  7,  a  season  for  the 
VESYS  A  solutions  is  usually  taken  to  represent  a  quarter  of  a  year 
and  the  corresponding  temperature  is  arrived  at  by  averaging  the 
mean  temperatures  for  those  three  months  comprising  the  season.   It 
is  felt  that  the  use  of  these  average  air  temperatures  modified  to 
pavement  temperatures  by  appropriate  ratios  representing  the  same 
seasons  will  provide  reasonably  accurate  pavement  temperatures.   The 
ratios  of  pavement  temperature  to  air  temperature  derived  from  this 
limited  study  appear  in  Table  4  on  a  seasonal  basis.   It  was  found  that 
these  ratios  (or  multipliers)  tend  to  grow  excessively  in  magnitude  as 
the  air  temperature  approaches  zero,  so  differences  in  temperature 
were  used  for  these  lower  temperatures  to  arrive  at  the  best  estimates 
for  pavement  temperatures.   For  instance,  based  on  the  limited  data 
available,  it  was  decided  that  the  pavement  temperature  was  best  ar- 
rived at  for  the  AASHO  Road  Test  sections  by  adding  27°F  (15°C)  to  the 
air  temperature  if  the  pavement  temperature  was  less  than  40°F. 
(4.4°C).   If  the  pavement  temperature  was  more  than  40°F  (4.4°C)  only 
24°F  (13.3°C)  was  added. 

The  data  derived  from  this  study  of  pavement  versus  air 
temperatures  plus  the  curves  appearing  in  Ref  19  were  used  to  con- 
vert air  temperatures  to  pavement  temperatures  for  this  project. 
The  resulting  pavement  temperatures  appear  in  Table  5.   It  should 
be  remembered  that  these  are  mean  pavement  temperatures  and  that 
they  may  be  much  greater  around  mid-day  and  much  less  at  night. 

Seasonal  Permanent  Deformation  Characterizations 

The  ability  of  VESYS  A  to  consider  seasonal  permanent  deforma- 
tion characterizations  for  each  of  the  layers  probably  represents 
the  most  important  improvement  to  its  ability  to  simulate  a  real 
flexible  pavement  system.   It  is  necessary,  however,  in  order  to 
accrue  any  benefit  from  this  improvement,  to  define  how  the  perman- 
ent deformation  potential  of  the  various  layer  materials  vary  season- 
ally.  As  the  permanent  deformation  parameters  ALPHA  and  GNU  are 
quite  new  to  our  profession,  very  little  data  exists  as  to  either 
the  correct  values  that  should  be  used  in  general  or  how  these 
values  change  on  a  seasonal  basis.   Consequently,  it  was  necessary 
to  conduct  a  detailed  study  for  the  material  types  in  each  layer  in 
order  to  establish  realistic  input  value  for  the  parameters. 
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Table  5.   Pavement  Temperatures  by  Season 


Pavement 
Identification 

Pavement  Temperature,  °F" 

Winter 

Spring 

Summer 

Fall 

AASHO  Sections 

266,  272,  308  and  470 

and 
Illinois  IH-80 
Section  38  (Flex  7) 

51 

69 

96 

73 

Brampton  Test  Road 
Sections  1,  4,  11  and 
17 

-2 

42 

75 

51 

Utah  IH-80N 

56 

71 

95 

65 

Florida  IH-10 

and 
Wet-No  Freeze, 
Factorial  Study 

78 

95 

106 

94 

Texas  IH-20 

and 

Dry-No  Freeze, 
Factorial  Study 

54 

76 

94 

75 

Wet-Freeze, 
Factorial  Study 

35 

66 

92 

60 

Dry-Freeze, 
Factorial  Study 

38 

62 

89 

49 

On    = 


C  =  0.5554  (°F-32) 
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These  studies  and  their  results  are  described  in  detail  in 
Appendices  A,  B,  and  C  for  asphalt  concrete,  typical  base  materials 
and  typical  clay  subgrade  materials,  respectively.   As  discussed  in 
these  Appendices,  the  trends  for  seasonal  variation  in  these  input 
parameters  were  established  through  study  of  papers  available  that 
reported  the  results  of  long-term  repetitive  load  tests  on  various 
materials  in  such  format  that  ALPHA  and  GNU  could  be  determined,  from 
previous  testing  reported  in  Reference  1  and  other  testing  conducted 
for  this  project.   The  resulting  conclusions  appear  at  the  end  of 
each  of  these  Appendices  and  are  given  below: 

1.  Asphalt  concrete: 

a.  Permanent  strain  is  extremely  temperature  dependent  as 
indicated  by  project  testing  and  by  all  other  test  re- 
sults in  the  literature.   As  an  example,  the  permanent 
strain  at  100°F  (37.8°C)  was  around  30  times  that  for 
70°F  (21.1°C)  at  a  stress  of  30  psi  (207  kN/m2)  and 
100,000  load  cycles. 

b.  GNU(l)  and  ALPHA (1)  both  decreased  with  increasing 
stress.   Decrease  in  ALPHA(l)  dominates  to  cause  an 
important  increase  in  permanent  strain  with  increasing 
stress. 

c.  It  appears  that  the  mechanisms  affecting  permanent  strain 
differ  above  and  below  some  temperature  in  the  order  of 
60°F  (J5.6°C)  to  70°F  (21.1°C).   Below  15. 6°C,  GNU(l)  in- 
creases with  temperature  and  ALPHA(l)  increases  slightly 
with  temperature  for  the  higher  stress  levels,  but  de- 
creases at  the  lower  stress  levels.   Above  15.6°C,  both 
ALPHA(l)  and  GNU(l)  decrease  with  increasing  temperature. 
The  rate  of  increase  in  permanent  strain  with  increasing 
temperature  is  much  higher  after  15.6°C  than  it  is  at  the 
lower  temperatures. 

d.  The  dynamic  modulus  is  clearly  somewhat  stress  dependent. 
Above  the  15.6°C  to  21.1°C  zone,  the  dynamic  modulus  in- 
creased with  stress  level.   Below  15.6°C,  it  decreased 
markedly  with  increased  stress  level. 

e.  The  dynamic  modulus  has  a  well  defined  negative  cor- 
relation with  temperature.   It  decreases  from  over  1 
million  psi  (6,900,000  kN/m2)  at  low  temperatures  to  less 
than  100,000  psi  (690,000  kN/m2)  at  higher  temperatures. 
This  is  not  new  information,  but  has  been  noted  in  all 
papers  in  recent  years  that  report  the  results  of  re- 
petitive-load tests  on  asphalt  concrete. 

2.  Base  layer  materials: 

a.   ALPHA(2)  for  granular  base  materials  appear 

to  be  little  effected  by  varying  stress  levels. 
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b.  Crushed  limestone  will  generally  have  less  permanent 
strains  than  gravel  under  the  same  compactive  energy 
and  stress  conditions. 

c.  Permanent  strains  are  greatly  increased  by  increasing 
moisture  content.   For  example,  an  increase  from  4.7 

to  5.5  percent  moisture  content  resulted  in  an  increase 
in  permanent  strain  for  one  soil  from  0.0015  inches  per 
inch  to  0.0071  inches  per  inch. 

d.  The  resilient  modulus  for  dry  crushed  stone  may  be  ex- 
pected to  be  much  higher  than  that  for  the  same  crushed 
stone  with  capillary  saturation  or  inundation. 

e.  ALPHA(2)  appears  to  be  essentially  independent  of  both 
deviator  stress  and  lateral  stress,  but  GNU(2)  increases 
with  deviator  stress,  thus  causing  an  increase  in  per- 
manent strain.   GNU (2)  decreases  with  increasing  lateral 
pressure,  thus  causing  decreased  permanent  strain. 


3.   Silty  clay  subgrades: 

a.  GNU (3)  generally  increases  as  moisture  content  increases 
and  decreases  as  density  increases.   An  apparent  sharp 
decrease  in  GNU (3)  may  occur  if  sufficient  moisture 
exists  to  accept  pore  pressure. 

b.  ALPHA(3)  may  increase  or  decrease  somewhat  with  in- 
creasing moisture  content,  depending  on  the  nature  of 
the  soil  and  the  amount  of  moisture  present.   The  varia- 
tion due  to  moisture  appears  to  depend  on  clay  content. 
ALPHA(3)  for  the  silty  clays  seemed  to  be  almost  inde- 
pendent of  moisture  while  ALPHA(3)  for  the  "heavier" 
clay  appears  to  vary  much  more  with  moisture  content. 

c.  Both  ALPHA(3)  and  GNU (3)  increase  together  as  moisture 
increases  until  ALPHA(3)  reaches  some  level  (higher 
maxima  for  higher  density),  and  then  ALPHA(3)  continues 
to  increase  while  GNU (3)  sharply  decreases.   This  im- 
plies that  at  some  point  of  moisture  increase  the 
fraction  of  total  strain  that  is  permanent  decreases 
(fraction  of  resilient  strain  increases)  although  the 
magnitude  of  permanent  strain  is  increasing.   This 
probably  is  due  to  limitations  of  a  constant  triaxial 
pressure  to  simulate  resistance  to  large  strains 
(consistent  with  higher  moisture  content)  imposed  by 
subgrade  soils  in  the  field.   Stated  differently,  the 
passive  case  for  coefficient  of  lateral  earth  pressure 
likely  applies  after  many  wheel  loadings  in  the  field 
so  that  cell  pressure  needs  to  increase  as  the  test 
progresses  in  order  to  simulate  these  conditions  ac- 
curately. 
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d.  ALPHA(3)  appears  to  be  relatively  independent  of 
deviator  stress  if  the  clay  content  is  low,  but 
may  vary  considerably  as  clay  content  increases. 

e.  GNU (3)  appears  to  be  relatively  independent  of  de- 
viator  stress  at  lower  moisture  contents,  but  de- 
creases sharply  at  higher  moisture  contents  as  deviator 
stress  increases  (possibly  due  to  pore  water  accepting 
some  of  the  stress  increase) . 

f.  Typical  values  of  ALPHA(3)  appear  to  range  from  0.65 
to  0.95. 

g.  Typical  values  of  GNU(3)  appear  to  range  from  0.01 
to  0.40,  but  can  range  much  higher  for  unusual  cases 
or  materials. 

The  information  on  variations  of  ALPHA  and  GNU  for  the  three 
layers  were  combined  with  data  furnished  by  State  Departments  of 
Transportation  in  Utah  and  Florida,  data  available  from  Reference  1 
for  the  AASHO  Road  Test  and  typical  data  for  sandy  soils  for  the 
Texas  sections  were  used  to  arrive  at  suitable  initial  values  for 
ALPHA  and  GNU.   These  values  were  later  modified  in  some  cases  as  a 
result  of  iterative  study  of  cause  and  effect  engaged  in  during  the 
calibration  phase  of  this  project.   This  iterative  study  and  these 
modifications  will  be  discussed  in  more  detail  in  subsequent  chapters, 
The  final  values  used  for  ALPHA  and  GNU  appear  in  Table  6. 
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Seasonal  Layer  Stiffness  Characterizations 

The  use  of  seasonal  multipliers  to  modify  creep  compliance 
curves  on  a  seasonal  basis  was  discussed  briefly  in  the  introduction 
and  in  a  little  more  detail  in  Chapter  II.   The  use  of  multipliers 
to  modify  the  creep  compliance  curve  has  the  effect  of  maintaining 
the  shape  of  the  creep  compliance  curve,  but  increasing  or  decreasing 
the  magnitudes  of  the  creep  compliance  array.   While  there  is  un- 
doubtedly some  variation  in  the  shape  of  the  curve  with  seasonal 
changes  and  moisture  content  or  if  the  base  material  is  frozen,  the 
creep  compliance  for  base  and  subgrade  materials  is  not  heavily 
dependent  on  load  duration,  so  maintaining  a  constant  shape  should 
have  little  effect  on  the  values  of  creep  compliance  used. 

There  is  not  appreciable  information  in  the  literature  as  to  the 
variation  of  creep  compliance  for  subgrades  and  base  material  with 
moisture  change  or  other  variables.   However,  there  is  considera- 
ble body  of  information  available  for  the  variation  in  resilient 
modulus  with  various  changes  in  moisture  content,  density  and  freeze- 
thaw  cycles.   If  one  arrives  at  a  set  of  appropriate  multipliers  in 
terms  of  resilient  modulus,  the  inverse  of  the  multiplier  will  be 
the  appropriate  multiplier  for  the  creep  compliance  array.   Con- 
sequently, it  is  possible  to  utilize  the  data  available  for  varia- 
tion in  resilient  modulus  directly  to  arrive  at  multipliers  for 
creep  compliance  and  this  was  done  on  this  project. 

As  the  seasonal  variation  in  creep  compliance  (or  layer  stiff- 
ness) may  only  be  arrived  at  reasonably  in  terms  of  the  variation 
in  the  nature  of  the  material  due  to  some  environmental  process, 
it  is  necessary  to  define  the  variation  in  the  material  prior  to 
arriving  at  rational  values  of  the  multipliers.   The  most  important 
single  variation  in  the  material  to  be  considered  is  change  in  mois- 
ture content.   A  study  of  the  literature  was  conducted  to  arrive  at 
the  incidence  of  moisture  change  at  subgrade  materials  from  as  wide 
a  base  of  data  as  possible.   Some  of  the  references  included  in  the 
study  were  References  2Ql,  21^,  22^  and  23  .   Reference  21  was  especially 


Haliburton,  T.A. ,  "Subgrade  Moisture  Variations",  Final  Report, 
Oklahoma  Research  Program,  Project  64-01-3,  School  of  Civil  Engineer- 
ing, Oklahoma  State  University,  August  1970. 

2 
Dempsey,  B.J.,  "Climatic  Effects  on  Airport  Pavement  Systems 

State  of  the  Art",  Final  Report  FAA-RD-7 5-196,  U.S.  Army  Engineer 
Waterways  Experiment  Station,  June  1976. 
3 
Bergan,  A.T.  and  C.L.  Monismith,  "Characterization  of  Subgrade 

Soils  in  Cold  Regions  for  Pavement  Design  Purposes",  Highway  Research 
Record  431,  1973. 

Cumberledge,  G.,  Cominsky,  R.  J.  and  Bhajandas,  A.C.,  "Moisture 
Variation  Beneath  Highway  Pavements  and  the  Associated  Changes  in 
Subgrade  Support",  Commonwealth  of  Pennsylvania,  Department  of  Trans- 
portation, 1973. 
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useful  as  it  was  itself  a  detailed  summary  of  the  majority  of  the 
information  available  on  the  subject.   The  following  was  concluded 
from  this  study  and* includes  only  those  study  results  of  specific 
interest  to  this  project: 

1.  There  are  few  clear  trends  of  a  general  nature  that  may  be 
applied  to  predict  moisture  content  on  a  seasonal  basis  at  a 
high  confidence  level  for  any  pavement  subgrade  or  base  other 
than  those  upon  which  the  data  was  based.   The  moisture  content 
in  terms  of  seasonal  trends  not  only  varies  from  site  to  site, 
but  may  vary  from  one  year  to  the  next  at  the  same  site.   Con- 
sequently, a  prediction  of  moisture  content  has  to  be  based  on 
the  particular  environmental  conditions  at  the  site  of  interest, 
the  nature  of  the  subgrade  material,  the  location  of  the  water 
table  and  other  parameters. 

2.  Several  researchers  have  found  relationships  between  the 
equilibrium  moisture  content  under  the  pavement  and  the  plastic 
limit  of  the  subgrade  soil.   Conversely,  no  good  relationship 
has  been  developed  between  precipitation  and  moisture  content 
in  the  subgrade. 

3.  In-situ  moisture  content  for  clay  subgrades  will  generally 
exceed  the  optimum  moisture  content.  The  opposite  is  true  for 
sands . 

4.  Several  researchers  have  found  that  the  maximum  moisture 
contents  occur  during  the  winter  or  spring  months;  however, 
this  may  not  be  true  in  locations  such  as  Florida  where  sandy 
subgrades  prevail  and  periods  of  high  rainfall  occur  generally 
in  different  seasons  from  most  of  the  United  States. 

5.  Magnitudes  of  seasonal  moisture  changes  are  generally 

in  the  order  of  two  or  three  percent,  but  such  small  changes 

in  magnitude  of  moisture  content  can  result  in  very  considerable 

changes  in  layer  stiffness. 

6.  There  is  a  general  trend  toward  increase  in  moisture  content 
with  time. 

7.  Moisture  content  correlates  better  with  rainfall  after 

the  pavement  surface  becomes  cracked  so  that  moisture  may  enter 
the  base  and  subgrade  from  above. 

8.  It  is  important  to  know  the  depth  of  frost  penetration  as 
the  stiffness  is  greatly  increased  for  a  layer  when  it  is  frozen. 
Also,  the  moisture  content  will  generally  be  quite  high  during 
the  spring  thaw,  thus  greatly  decreasing  the  stiffness  (i.e., 
increasing  the  magnitudes  of  creep  compliance)  for  this  season. 

As  there  were  no  truly  general  trends  that  could  be  used  to  develop 
seasonal  multipliers,  it  became  necessary  for  specific  pavements  to 
draw  on  regional  data  wherever  possible.   For  example,  Ref .13  provides 
data  on  frost  penetration  depths,  seasonal  moisture  contents  of  subgrade 
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and  base  and  measured  values  of  elastic  modulus  for  several  sections 
during  the  spring  and  summer  of  1960.   There  were  also  seasonal  measure- 
ments of  elastic  modulus  using  plate  bearing  tests  on  the  base,  subbase 
and  embankment  over  a  period  of  nearly  two  years.   This  latter  gave  a 
rather  direct  variation  in  stiffness  on  a  seasonal  basis  and  was  used 
as  a  basis  for  defining  the  seasonal  multipliers.   For   the  other 
locations,  regional  information  on  changes  in  moisture  content  on  a 
seasonal  basis  were  obtained  to  the  extent  possible  and  seasonal 
multipliers  developed  on  the  basis  of  anticipated  variation  in  resil- 
ient modulus  as  a  function  of  variation  in  moisture  content. 

There  are  also  some  variations  in  dry  density  that  occur  seasonally, 
especially  where  deep  frost  penetration  and  a  serious  spring  breakup 
occurs.   While  these  variations  may  only  be  in  the  order  of  one  pcf , 
the  effect  on  the  layer  stiffness  can  be  worthy  of  consideration. 

The  values  of  seasonal  multipliers  resulting  from  this  study  are 
listed  in  Table  7.   These  seasonal  multipliers  were  applied  by  VESYS  A 
to  modify  seasonally  the  constant  arrays  of  creep  compliance  values 
appearing  in  Table  8.   (The  development  of  values  of  creep  compliance 
in  the  next  chapter) . 

Low  Temperature  Cracking 

Input  values  for  low  temperature  cracking  were  required  for  all 
pavement  sections  analyzed,  except  for  Texas  IH-20,  Elorida  IH-10,  wet- 
no  freeze  and  dry-no  freeze  cases  where  serious  frost  penetration  was 
not  expected.   As  the  age  of  pavement  and  thickness  of  the  bituminous 
layer  are  otherwise  input  to  the  computer  program,  the  only  additional 
information  required  is  a  numerical  code  representing  subgrade  soil 
type  (2  for  clay,  3  for  loam  and  5  for  sand),  the  winter  design 
temperature  in  °C_ and  the  original  asphalt  cement  stiffness. 

The  subgrade  soil  type  may  be  readily  identified  from  knowledge 
of  the  site,  but  the  other  values  are  somewhat  more  complex  and  are 
arrived  at  as  described  in  Appendix  G.   The  values  used  for  these 
studies  appear  in  Table  9. 
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Table  7.   Seasonal  Multipliers  for  Creep 

Compliance  Arrays,  Base  and  Subgrade 


Pavement 
Identification 

Layer 
Identification 

Seasonal  Multipliers 

Winter 

Spring 

•  Summer 

Fall 

AASHO  Sections 
266,  272,  308 
and  470 

Base 

0.88 

1.23 

1.00 

0.90 

Subgrade 

0.86 

1.23 

1.00 

0.88 

Brampton  Sections 
3,4,11  and  17 

Base 

0.80 

1.25 

1.00 

0.90 

Subgrade 

0.80 

1.25 

1.00 

0.90 

Illinois  IH-80, 
Section  38 
(Flex  7) 

Base 

0.88 

1.23 

1.00 

0.90 

Subgrade 

0.86 

1.23 

1.00 

0.88 

Utah  IH-80N 

Base 

0.90 

1.20 

1.00 

0.95 

Subgrade 

0.95 

1.20 

1.00 

0.95 

Florida  IH-10 

Base 

1.00 

0.77 

1.47 

1.23 

Subgrade 

1.00 

1.00 

1.00 

1.00 

Texas  IH-20 

Base 

1.10 

1.05 

1.00 

1.05 

Subgrade 

1.00 

1.00 

1.00 

1.00 

Wet-Freeze, 
Factorial  Study 

Base 

0.88 

1.23 

1.00 

0.90 

Subgrade 

0.90 

1.15 

1.00 

0.92 

Dry-Freeze, 
Factorial  Study 

Base 

0.90 

1.20 

1.00 

0.95 

Subgrade 

0.95 

1.10 

1.00 

0.95 

Wet-No  Freeze, 
Factorial  Study 

Base 

0.68 

0.52 

1.00 

0.84 

Subgrade 

1.00 

1.00 

1.00 

1.00 

Dry-No  Freeze, 
Factorial  Study 

Base 

1.10 

1.05 

1.00 

1.05 

Subgrade 

1.00 

1.00 

1.00 

1.00 
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Axle  Load  Distributions  for  the  Four 
Legal  Axle  Loads 


The  improvements  incorporated  into  VESYS  A,  the  detailed  develop- 
ment of  the  input  data  and  the  careful  construction  of  the  study 
factorial  were  all  critical  to  the  success  of  this  research  project, 
but  meaningful  results  also  depend  rather  directly  on  identification 
of  axle  load  distributions  representative  of  those  that  would  occur 
for  the  higher  legal  limits.   This  is  a  complex  question  that  is 
strongly  dependent  on  such  specifics  as  1)  allowable  dimensions  of 
equipment,  2)  economics  of  buying  new  equipment  to  take  advantage 
of  larger  weight  limits,  3)  density  of  commodities  carried,  4)  density 
mix  of  commodities  carried,  5)  applicable  bridge  formulas,  6)  user 
tax  structure  and  many  others.   Because  of  the  complexity  of  the 
question  and  its  importance,  axle  load  distributions  for  use  in  this 
study  were  developed  after  careful  consideration  of  the  literature, 
correspondence  with  agencies  representing  the  trucking  industry  and 
study  of  W-4  tables  from  a  number  of  states,  including  those  with 
relatively  higher  legal  axle  limits. 

Fortunately,  this  question  had  been  considered  in  detail  pre- 
viously by  Winfrey  and  others  (References  21-*-  and  22^)  and  methods 
developed  for  predicting  these  axle  loads.   Typical  axle  load  dis- 
tributions grouped  by  single  and  tandem  axles  (See  Figure  1)  appear 
in  Reference  24  and  include  data  from  states  having  18,  20  and  22-kip 
(8100,  9000,  and  9900  kg)  single  axle  limits  and  related  tandem  axle 
limits  of  32,  35  and  38-kips  (14440,  15750  and  17100  kg),  respectively. 
Predicted  axle  load  distributions  for  24-kip  (10,900  kg)  single  and 
41-kip  (18,600  kg)  tandem  axles  are  also  included.   Further  study 
appeared  to  reinforce  the  validity  of  these  axle  load  distributions, 
so  the  decision  was  made  to  combine  these  single  and  tandem  axle  load 
distributions  on  a  rational  basis  and  to  use  them  for  this  study. 

In  order  to  combine  the  single  and  tandem  axle  load  distributions, 
it  was  necessary  to  arrive  at  the  ratio  of  each  to  the  total  number  of 
axle  loads  in  a  typical  truck  traffic  stream.   Study  of  W-4  tables 
from  California,  New  York,  Maryland,  Pennsylvania,  New  Jersey,  South 


Winfrey,  R.,  and  others,  "Economics  of  the  Maximum  Limits  of 

Motor  Vehicle  Dimensions  and  Weights",  Volume  1,  Report  No.  FHWA-RD- 

73-69,  September  1968. 

? 
Whiteside,  R.  E.,  T.Y.  Chu,  J.  C.  Cosby,  R.  L.  Whitaker,  and 

R.  Winfrey,  "Changes  in  Legal  Vehicle  Weights  and  Dimensions,  Some 

Economic  Effects  on  Highways",  NCHRP  Report  141,  1973. 
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Dakota,  Kansas  and  Texas  led  to  a  typical  axle  mix  of  52  percent 
tandem  axles  to  48  percent  single  axles.   By  applying  this  ratio 
and  multiplying  the  number  of  tandem  axles  by  2.25  (in  lieu  of  2.0 
to  consider  approximately  the  superposed  effects  from  the  adjacent 
axle)  to  convert  them  to  approximate  equivalent  single  axles,  the 
axle  load  distributions  appearing  in  Table  11  were  derived.   These 
values  were  used  for  the  study. 

The  axle  load  distributions  shown  in  Table  11  were  each  divided 
into  four  discrete  segments  for  separate  consideration  in  VESYS  A 
as  shown  in  Table  10. 


TABLE  10-   Discretized  Representations  of  Legal  Axl< 

a 

Load  Distributions  by  Percent  of  Total 

Legal  Axle 

1  kg  =  0.0022  kips 

Load  Limit 

(Kips) 

2 

to  8  kips 

8  to  12  kips   12  to  18  kips 

18  to  28  kips 

18 

34.1 

41.9          23.5 

0.5 

20 

30.4 

41.6          26.0 

2.0 

22 

27.2 

40.7          29.1 

3.0 

24 

23.7 

41.0          31.5 

3.8 

It  was  assumed  that  tire  pressures  for  the  axle  load  increments 
would  be  75  psi  (5.2  kg/cm2)  for  the  2  to  8-kip  increment,  80  psi 
(5.5  kg/cm2)  for  the  8  to  12-kip  increment,  85  psi  (5.9  kg/cm2)  for 
the  12  to  18-kip  increment  and  90  psi  (6.3  kg/cm2)  for  the  18  to 
28-kip  increment.   These  estimated  values  were  based  on  general  data 
obtained  from  queries  to  agencies  representing  the  trucking  industry. 
Calculations  were  made  for  each  discrete  increment  in  Table  10, 
using  all  included  increments  in  Table  11  to  arrive  at  a  weighted 
mean  value  of  radius  of  area  loaded  by  a  single  wheel  load  (dual 
wheels  considered  as  a  single  wheel)  for  each  legal  axle  load  and 
each  discrete  increment  of  the  axle  load  distribution.  A  mean  value 
of  the  four  mean  radii  calculated  for  each  axle  load  increment  in 
Table  11  was  then  obtained  as  a  common  radius  for  all  legal  axle 
load  limits.   This  was  done  to  restrict  the  number  of  internal 
static  solutions  required  without  loss  of  accuracy.   The  values  of 
AMPLITUD  (approximate  equivalent  of  tire  pressure)  were  the  modified 
slightly  to  insure  that  the  correct  weighted  wheel  load  resulted  for 
each  legal  axle  load.   The  values  of  radii,  AMPLITUD,  and  resulting 
values  of  wheel  load  used  appear  in  Table  12. 

Fatigue  Potential  Characterizations 

Realistic  prediction  of  fatigue  cracking  damage  is  necessary  as 
the  fatigue  cracking  mode  of  distress  is  the  one  that  most  frequently 
creates  the  need  for  overlaying  pavements.   This  means  that  the  materials 
characterizations  used  in  the  classical  fatigue  equation  must  be  realistic 
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It  has  been  quite  common  practice  to  base  fatigue  cracking 
predictions  on  materials  characterizations  derived  from  laboratory 
fatigue  testing.   The  inadequacies  of  standard  laboratory  fatigue 
tests  for  simulating  the  stress  state  in  and  the  fatigue  response 
of  a  material  in  a  real  supported  surface  layer  have  been  discussed  in 
detail  (Refs.  26^  and  272).   In  general,  laboratory  test  such  as  stand- 
dard  beam  tests  relate  to  cracking  at  the  bottom  of  a  pavement  whereas 
the  cracking  distress  of  interest  is  the  appearance  .of  cracks 
the  surface.  Wheel  Tracking  Tests  reported  by  Van  Dijk  (Ref.  27), 
indicate  that  it  takes  around  8  to  10  times  as  many  wheel  loads 
to  propagate  a  crack  to  the  surface  as  it  takes  to  initiate  a 
crack  at  the  bottom  of  the  surface  layer.   Consequently,  fatigue 
characterizations  generally  must  be  related  to  studies  of  real 
pavements  or  careful  simulations  in  the  laboratory  to  be  meaning- 
ful for  predictive  use. 

Kingham  (Ref  28)  and  Austin  Research  Engineers  (Ref  29)  have 
independently  conducted  multiple  regressions  on  AASHO  Road  Test  data 
for  cracking  appearing  at  the  surface.   The  values  for  the  exponent 
K2  were  essentially  identical,  but  Kingham fs  coefficient  K^,  was 
somewhat  larger.  Values  of  K^  =  2.6  x  10~13  and  K2  -  5.0  were  used 
to  approximate  the  fatigue  characterizations  from  the  AASHO  Road  Test 
at  a  temperature  of  70°F  (21.1°C). 

The  seasonal  values  of  K,  (T)  and  K~(T)  were  arrived  at  by 
using  the  procedures  described  in  References  1  and  21  for  convert- 
ing K- (70°F)  and  K  (70  F)  to  other  temperatures,  using  the  pave- 
ment temperatures  shown  in  Table  5.   The  values  of  K.(T)  and  K.(T) 
used  appear  in  Table  13. 


Pavement  Structure,  Truck  Traffic  and  Other  Input  Values 

It  was  necessary  to  select  typical  low  and  high  levels  for  both 
pavement  structure  and  truck  traffic  as  these  were  factorial  variables, 
It  was  decided  to  use  the  Texas  IH-20  pavement  structure  and  truck 
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TABLE  13.  Fatigue  life  potential  characterizations 
for  the  four  environmental  zones  and  the 
four  seasons 


Pavement 

Environmental 

Temperature 

Zone 

Seasons 

(°F) 

Kl« 

:t> 

K2(T) 

Winter 

35 

2.6 

X 

io-" 

5.18 

Spring 

66 

1.7 

X 

io-" 

5.02 

Wet-Freeze 

Summer 
Fall 

92 
60 

5.2 
4.4 

X 
X 

10'11 

-14 

10  q 

-15 

4.89 
5.05 

Winter 

38 

3.1 

X 

10  r. 

5.16 

Spring 

62 

7.3 

X 

10 

-14 
10 

5.04 

Dry-Freeze 

Summer 
Fall 

89 
49 

2.6 
1.0 

X 
X 

4.91 
5.11 

Winter 
Spring 

78 
95 

1.0 
1.6 

X 
X 

10  a 

io:!o 

4.96 
4.88 

Wet-No  Freeze 

Summer 

106 

2.6 

X 

4.82 

Fall 

94 

1.2 

X 

10 
-14 

4.88 

Winter 
Spring 

54 

76 

1.8 
5.7 

X 
X 

10-13 

10-10 

10-13 
10 

5.08 
4.97 

Dry-No  Freeze 

Summer 
Fall 

94 
75   ' 

1.2 

4.7 

X 

X 

4.88 
4.98 

°C  =  0.5554  (°F  -32) 
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traffic  as  low  levels  and  those  for  Illinois  IH-80  as  high  values. 
This  resulted  in  a  "thin"  pavement  section  having  3.9  inches  (9.91 
cm)  of  surface  and  19.0  inches  (48.3  cm)  of  base,  and  a  "thick"1  pave- 
ment section  having  4.5  inches  (11.4  cm)  of  surface  and  31.5  inches 
(80.0  cm)  of  base.   The  truck  traffic  values  used  appear  in  Table  14. 
These  values  are  based  on  real  data  extrapolated  rationally  to  20 
years.   Where  overlays  were  added,  the  extrapolation  was  continued 
past  20  years  to  allow  prediction  of  expected  life. 

The  same  distributions  of  truck  traffic  in  time  were  used  for 
all  legal  axle  loads,  which  implies  an  increase  in  tonnage  with  in- 
creasing legal  axle  loads.   An  alternative  approach  could  have  been 
decreases  in  traffic  with  increases  in  axle  load  distributions  to 
maintain  constant  tonnage,  but  this  was  the  less  severe  condition 
and  the  relative  realism  of  either  extreme  is  moot. 

Creep  compliance  arrays  for  the  asphalt  concrete  were  developed 
for  each  of  the  pavement  sections  used  in  the  comparative  studies 
and  for  the  factorial  study.   The  values  in  these  arrays  were  based 
on  laboratory  creep  testing  in  some  cases,  resilient  modulus  testing 
results  plus  laboratory  creep  testing  or  judgment  for  a  few  cases 
based  on  comparisons  of  mix  characteristics  to  others  for  which 
test  data  was  available.   The  creep  compliance  arrays  used  appear 
in  Table  15. 

Other  fixed  values  used  as  problem  input  and  not  discussed 
separately  above  appear  in  Table  16.   These  values  were  generally 
obtained  from  engineering  experience  and  previous  studies  aimed 
at  developing  realistic  input  values  for  VESYS  IIM. 
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TABLE  14.   Truck  traffic  in  Axles  per  Day 


Time  in  Years 
from  Construction 
(TRANDOM) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Axles  per  Day 

(LAMBDA) 

Low  Traffic 

High  Traffic 

1770 

2340 

1810 

4000 

1860 

5220 

1910 

6470 

1970 

7660 

2030 

8550 

2100 

8750 

2170 

9260 

2250 

9450 

2330 

9700 

2410 

9800 

2480 

10000 

2590 

10150 

2690 

10300 

2790 

10450 

2890 

10600 

2995 

10750 

3100 

10900 

3210 

11050 

3340 

11200 
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Table  16.   Fixed  Input  Values 


Input  Name 


Description 


Input  Value 


BETA 
C0EFK1 

C0EFK2 
CORLCOEF 

CORLEXP 

DURATION 

K1K2C0RL 

NTEMPS 

NTRADOM 
NTSTATIC 

PSIFAIL 

QUALITYO 
REFTEMP 


Time-temperature 
shift  function 

Coefficient  of  variation 
for  fatigue  coefficient 

Kl 

Coefficient  of  variation 
for  fatigue  exponent  K2 

Materials  coefficient  B 
for  the  spacial  auto- 
correlation function 

Materials  exponent  C  for 
the  spatial  auto-correla- 
tion function 


0.054 


0.30 


0.04 


1.00 


0.06 


Wheel  load  duration  at  a 

point,  seconds  0.015 

Correlation  coefficient       -0.867 
for  K^  and  K2 

Number  of  points  in 

temperature  array  4 

Number  of  years  considered     20 

Number  of  points  in  creep 
compliance  arrays  11 

Level  of  serviceability 
index  considered  unaccep- 
table 2.5 

Initial  serviceability  index   4.2 

Reference  temperature  for 

master  creep  compliance 

curves,  F.  70 
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TABLE  16. 


Fixed  Input  Values  (continued) . 


Input  Name 


Description 


Input  Value 


STDEVO 


Standard  deviation  of 
initial  serviceability 
index 


0.25 


TOLERNCE 


VARC0EF1 


VARC0EF2 


VARC0EF3 


VCAMP 


VCDUR 


Minimum  acceptable 
reliability,  percent 

Coefficient  of  variation 
for  surface  layer  creep 
compliance  array 

Coefficient  of  variation 
for  base  course  creep 
compliance  array 

Coefficient  of  variation 
for  subgrade  creep  com- 
pliance array 

Variance  of  AMPLITUD  or 
tire  pressure 

Variance  of  DURATION  or 
wheel  load  duration 


50 


0.25 


0.30 


0.30 


5.4  x  10 


-6 


58 


CHAPTER  V 
CALIBRATION  OF  MODEL  FOR  REALISTIC  PREDICTIONS 

Several  possibilities  exist  for  calibration  of  a  mathematical 
model  such  as  VESYS  A  so  that  it  produces  calculated  responses  consis- 
tent with  the  real  responses  of  the  engineering  system  that  it  models. 
Three  of  these  considered  were:   1)  the  introduction  of  multipliers 
that  will  cause  calculated  responses  to  equal  measured  responses  of 
the  real  system,  2)  the  introduction  of  functions  that  will  transform 
the  calculated  responses  to  those  measued,  and  3)  development  of  input 
data  sufficiently  consistent  with  the  real  conditions  modeled  and  with 
the  nature  of  the  model  itself  such  that  the  calculated  responses  are 
realistic.   The  introduction  of  either  multipliers  or  functions  is  the 
most  direct  approach,  but  the  multipliers  are  too  simplistic  and  can- 
not be  expected  to  be  consistent  for  a  range  of  conditions  and  the 
development  of  appropriate  functions  would  represent  a  very  consider- 
able research  project  of  its  own.   Consequently,  the  procedure  employed 
for  this  project  was  an  iterative  procedure  involving  exercising  the 
model  to  arrive  at  predictions,  comparing  the  predictions  to  measured 
performance,  analyzing  the  differences  to  assess  their  cause  and 
rationally  modifying  the  problem  input  where  indicated  to  improve  pre- 
dictions.  Only  revisions  to  input  values  were  made  that  could  be 
justified  rationally  when  short-comings  of  previous  values  were  identi- 
fied through  analysis. 

The  initial  step  in  this  process  was  to  develop  input  data  for  a 
number  of  cases  for  which  measured  responses  were  available.   This 
included  four  sections  from  the  AASHO  Road  Test,  four  sections  from 
the  Brampton  Test  Road,  one  section  of  IH-80N  in  Utah,  and  one  section 
of  IH-10  in  Florida.   The  sections  of  IH-20  in  Texas  and  of  IH-80  in 
Illinois  were  also  included  for  study  although  measured  responses  were 
not  available.   Solutions  were  obtained  for  these  sections  and  compared 
to  measured  data.   Studies  were  then  initiated  to  determine  what  input 
variables  could  have  caused  differences  between  calculated  and  measured 
responses,  whether  the  input  values  represented  as  accurate  a  character- 
ization as  was  possible  and  whether  review  of  the  basis  for  arriving  at 
the  input  values  indicated  logical  modifications  that  might  be  expected 
to  give  more  accurate  calculated  responses.   Comparisons  were  also  made 
between  the  calculated  responses  and  the  input  data  for  the  entire 
group  of  calculations  in  order  to  determine  if  the  trends  and  variations 
between  calculated  and  measured  responses  were  consistent,  if  the  rela- 
tive values  of  inputs  were  consistent  for  the  various  cases  considered 
and  if  these  comparisons  identified  any  obvious  anomolies  that  should 
be  corrected. 

As  the  capability  for  considering  seasonal  characterizations  of 
materials  was  new,  it  also  became  necessary,  in  order  to  understand 
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differences  in  predicted  and  measured  rutting,  to  run  a  factorial  of 
separate  solutions  using  an  elastic  layered  program  called  ELSYM5. 
Deflections  were  calculated  for  the  various  seasonal  conditions  and 
general  equations  from  VESYS  A  applied  to  estimate  separately  the 
rutting  that  occurred  in  each  season  and  in  each  layer.   This  gave 
valuable  insight  that  was  used  to  improve  the  values  of  ALPHA  and  GNU 
for  the  four  seasons. 

The  effort  conducted  to  arrive  at  a  "calibrated"  test  factorial 
is  described  in  subsequent  sections  in  more  detail. 

Comparisons  to  AASHO  Road  Test  Measurements 

The  four  test  sections  from  the  AASHO  Road  Test  used  for  this 
study  were  those  previously  used  for  similar  comparisons  in  Reference 
1.   These  were  Sections  266,  272,  308  and  470.  Actual  materials  from 
the  AASHO  Road  Test  site  were  tested  during  the  research  culminating 
in  Reference  1  to  arrive  at  both  the  creep  compliance  arrays  in  Tables 
9  and  15  and  the  permanent  deformation  parameters.   Therefore,  a  very 
considerable  data  base  existed  already  for  these  four  sections,  but 
this  data  was  in  terms  of  the  requirements  for  VESYS  IIM  and  did  not 
include  the  sophistication  of  characterizations  in  terms  of  seasonal 
variations,  except  for  the  creep  compliance  arrays  for  the  asphalt 
concrete  surface  materials.  Also,  the  temperatures  used  for  the 
months  of  the  year  in  the  analyses  of  Reference  1  were  air  temperatures 
in  lieu  of  pavement  temperatures. 

The  pavement  temperatures  were  established  as  previously  described 
and  these  temperatures  then  became  the  basis  for  establishing  seasonal 
values  of  ALPHA(l)  and  GNU(l)  and  for  the  fatigue  parameters  K^  and 
K2.   These  values  plus  others  already  available  provided  the  materials 
characterizations  for  the  surface  layer. 

Fortunately,  Reference  13  includes  very  detailed  information  on 
all  pavement  layers  over  a  considerable  period  of  time.   Figure  97 
in  Reference  13  displays  the  results  in  terms  of  elastic  moduli  of 
plate  bearing  tests  on  the  base  course,  subbase  and  embankment  or 
subgrade  layers.   This  provides  very  directly  the  modifications  that 
should  be  made  seasonally  to  the  creep  compliance  arrays.  By  assuming 
the  elastic  moduli  during  summer  as  the  base  values  (i.e.,  having 
seasonal  multipliers  of  1.00),  mean  values  of  elastic  moduli  for 
other  seasons  may  be  divided  by  the  summer  modulus  to  arrive  at 
multipliers  to  convert  to  the  appropriate  stiffnesses  during  other 
seasons.   This  value  can  then  be  inverted  to  arrive  at  the  appropriate 
seasonal  multipliers  for  the  creep  compliance  arrays.   This  was  the 
approach  taken  to  arrive  at  the  seasonal  multipliers  for  the  AASHO 
Road  Test  sections  appearing  in  Table  7.  As  VESYS  A  is  a  three  layer 
model,  it  was  necessary  to  combine  on  a  weighted  basis  the  multipliers 
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for  the  base  and  subbase  to  obtain  the  seasonal  multipliers  for  the 
base  as  shown. 

Table  34  of  Reference  13  provided  variations  in  moisture  and 
density  between  Spring  and  Summer  seasons  for  the  base,  subbase  and 
subgrade.   Depths  of  frost  penetration  also  appeared  on  Page  11  as  a 
function  of  time.   Review  of  this  data  indicated  that  the  change  in 
moisture  content  for  the  subgrade  was  negligible  and  that  the  density 
variation  was  only  about  0.7  pcf  (11.2  kg/m^) .   A  change  in  moisture 
content  for  the  subbase  from  5.5  percent  to  4.8  percent  is  quite  signi- 
ficant, however,  as  is  a  change  from  4.3  percent  to  3.6  percent  for  the 
base.   These  variations  in  moisture  content  seasonally  and  the  presence 
of  frost  in  the  winter  were  considered  in  arriving  at  the  permanent  de- 
formation parameters  ALPHA  and  GNU  for  the  base  and  subgrade. 

The  most  difficult  task,  of  course,  was  the  establishment  of 
the  permanent  deformation  parameters  for  the  three  layers  as  functions 
of  seasons.   Because  no  precedents  existed  for  doing  this,  the  results 
of  the  studies  described  in  detail  in  Appendices  A,  B,  and  C  for  the 
surface  layer,  base  and  subgrade,  respectively,  were  used  for  estab- 
lishing these  values.   The  values  of  ALPHA  and  GNU  appearing  in 
Reference  1  were  used  as  a  starting  point  and  all  the  information 
available  was  applied  to  arrive  at  the  values  of  ALPHA  and  GNU  appearing 
in  Table  6. 

Figure  2  provides  measured  rut  depths  and  those  predicted  by 
VESYS  IIM  during  the  research  leading  to  Reference  1  for  AASHO  Sections 
266  and  308.   As  can  readily  be  seen,  the  previous  rut  depth  predictions 
were  very  much  lower  than  those  measured;  the  new  predictions  using 
VESYS  A  are  also  generally  lower,  but  are  much  closer  to  those  predicted, 
As  air  temperatures  were  used  in  lieu  of  pavement  temperatures  for  the 
AASHO  sections  previously,  this  does  not  represent  an  entirely  fair 
comparison.   Therefore,  a  new  run  was  made  for  Section  470  using  pave- 
ment temperatures  for  VESYS  IIM  and  these  appear  for  rut  depth,  slope 
variance  and  present  serviceability  index  predictions  in  Figure  3. 
Note  that  the  predictions  using  pavement  temperatures  in  VESYS  IIM  for 
rut  depth  are  considerably  better  than  those  for  air  temperatures.   A 
very  significant  additional  improvement  was  made  by  including  the 
seasonal  materials  characterizations  in  VESYS  A,  and  the  addition  of 
the  discretized  representation  of  the  axle  load  distribution  to  finish 
out  the  new  VESYS  A  capabilities  resulted  in  almost  perfect  predictions 
for  the  measured  rut  depth  for  this  section.   Looking  now  at  slope 
variance,  the  addition  of  pavement  temperature  for  VESYS  IIM  made  very 
little  difference,  but  the  improvements  in  VESYS  A  gave  a  considerably 
higher  prediction  of  slope  variance.  As  might  be  expected,  the  addi- 
tional rutting  and  slope  variance  predicted  by  VESYS  A  resulted  in  a 
considerable  reduction  in  serviceability. 
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62 


.6  i- 


a. 
O 


3 


Me 


Predicted   YESYS    A 


ed  by  Seasonal 

als   Characterizations 


cted   VESYS  HM 

Pavement   Temp. 

ted    VESYS  HM 
Air  Temp. 


1  cm  =   0.394   inch 


a. 


I  I  I  I 


1963     65       67       69       71        73 
Year,  19_ 
a.  Rut  depth  predictions   versus  measured  values 


3 

2 

fll 

<£> 

O 

O 

c 

— 

o 

X 

v. 

I 

o 

V) 

i . 

> 

c 
o 

a> 

T) 

a. 

o 

o 

or 

CO 

— 

& 


& 


ement    Temperature    VESYS   Ufa 


Temperature    VESYS  1M 


1963     65       67       69      71         73 

Year    I9_ 
b  Slope  variance  predictions  (no  measured  values  available) 


Predicted   VESYS  II    M 
Pavement  Temperature 


Improved  by  Seasonal 
Materials  Characteristics 


Predicted  VESYS    A 


Figure    3 


Year,  I9_ 
c.  Serviceability  predictions  (no  measured  values  available) 

Comparisons  of   measured    responses   and    those 

predicted   by   VESYS    IIM  and    VESYS   A,    AASHO 

Section   470.  r ., 

63 


The  comparisons  are  not  plotted  for  AASHO  Section  272  because 
the  predictions  using  air  temperature  in  VESYS  IIM  reported  in 
Reference  1  exceeded  the  measured  rut  depths  and  those  for  VESYS  A 
far  exceeded  the  rut  depths.   This  was  attributed  previously  to  a 
very  high  test  value  of  GNU(l)  and  this  is  still  believed  to  be  the 
cause  of  the  overpr edict ion. 

These  studies  for  the  AASHO  Road  Test  Section  provided  a  base 
for  continued  study  and  development,  indicated  clearly  that  the 
VESYS  A  model  gave  improved  predictions  as  expected,  showed  that  the 
VESYS  model  had  no  inherent  limitations  to  predicting  reasonable 
values  if  appropriate  materials  and  traffic  characterizations  were 
provided  and  provided  initial  indications  as  to  what  changes  in  input 
values  might  be  appropriate  subject  to  further  study. 

Comparisons  to  Brampton  Test  Road  Measurements 

Sections  3,  4,  11  and  17  of  the  Brampton  Test  Road  were  also 
included  in  the  previous  studies  and  reported  in  Reference  1.   The 
input  used  for  the  current  study  was  a  combination  of  that  from  the 
previous  study  and  carefully  derived  data  from  References  301,  31  , 
323,  334,  345,  356,  and  367. 
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The  construction  of  the  sections  including  depths  and  density 
of  the  layers  and  their  moisture  content  were  discussed  in  Reference 
29.   Section  3  was  basically  a  "full-depth"  asphalt  concrete  pavement 
with  11.5  inches  (29.2  cm)  of  asphalt  concrete  and  no  base  material. 
Section  4  could  be  called  "deep-strength"  with  9.5  inches  (24.1  cm) 
of  asphalt  concrete  and  6  inches  (15.2  cm)  of  select  granular  base. 
Sections  11  and  17  were  conventional  pavements  with  3.5  inches  (8.9 
cm)  of  asphalt  concrete.   Section  11  had  four  inches  (10.2  cm)  of 
crushed  rock  base  and  18  inches  (45.7  cm)  of  select  granular  subbase. 
Section  17  had  2  inches  (5.1  cm)  of  crushed  gravel  base  and  18  inches 
45.7  cm)  of  select  granular  base.   Traffic  information  and  performance 
data  for  the  first  four  years  of  operation  of  the  Brampton  Test  Road 
were  provided  by  Reference  31.   Further  performance  information  was 
provided  by  Reference  32.   Reference  33  provided  additional  performance 
data,  pavement  temperature  distribution  data,  resilient  moduli  of  the 
granular  base,  resilient  moduli  of  the  sand  subbase,  resilient  moduli 
of  the  subgrade  material,  stiffness  modulus  of  the  asphalt  concrete 
material  and  of  the  asphalt  cements,  stress  distributions  with  depth 
in  asphalt/concrete  layer  predicted  by  elastic  layer  theory  and  other 
very  useful  information  for  characterization  of  materials  for  VESYS  A. 
Reference  34  provided  information  in  detail  on  resilient  modulus 
testing  for  the  asphalt  concrete  surface  material,  asphalt  concrete 
base  materials,  crushed  gravel  base  course,  crushed  stone  base  course 
and  sand  cushion  subbase  course.   Reference  35  provided  resilient 
modulus  data  in  detail  on  the  subgrade  materials.   Reference  36  pro- 
vided information  on  a  substantial  test  effort  to  arrive  at  permanent 
deformation  characteristics  for  the  asphalt  base  materials  used  at 
the  Brampton  Test  Road.   The  wealth  of  information  available  in  the 
references  for  the  Brampton  Test  Road  was  considered  in  detail  in 
developing  the  seasonal  material  characteristics. 

The  first  solution  obtained  was  for  Section  3  and  the  results  for 
rut  depths  are  plotted  in  Figure  4  versus  measured  depth  and  the 
serviceability  in  Figure  8  versus  measured  serviceability.   The 
values  of  ALPHA (1)  and  GNU(l)  for  the  four  seasons  are  also  shown  on 
Figure  4  and  this  solution  is  identified  as  Run  1.   As  can  be  seen, 
rut  depth  is  very  considerably  overestimated  and  serviceability  under- 
estimated because  of  the  high  predictions  of  rut  depths.   There  ap- 
peared to  be  reason  to  reduce  GNU(l)  considerably  for  the  summer 
and  fall  seasons  and  to  increase  ALPHA (1)  for  all  seasons.   The 
results  of  this  appear  as  Run  2  in  Figures  4  and  8  and  amounted  to 
gross  underpredictions  of  rutting  and  overpredictions  of  service- 
ability.  For  Run  3,  ALPHA(l)  was  decreased  and  GNU(l)  was  slightly 
decreased  and  the  results  were  very  little  different  from  Run  2. 
Consequently,  it  appeared  that  GMJ(l)  was  the  variable  to  be  revised 
and  higher  values  were  used  for  the  summer  and  fall  seasons  as  indi- 
cated in  Figure  4  for  Run  4.   These  increases  in  GNU(l)  had  the  effect 
of  producing  much  improved  estimates  of  rut  depth  and  considerably 
improved  estimates  of  serviceability. 
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With  the  information  from  Section  3,  only  two  runs  each  were 
required  for  Sections  4,  11  and  17  and  plots  for  these  solutions  appear 
in  Figures  5,  6,  7,  9,  10,  and  11.   Higher  values  of  ALPHA(l)  wore 
generally  used  for  the  full-depth  and  deep-strength  sections  than 
for  the  conventional  sections  11  and  17  because  of  the  differences 
in  stress  levels  in  the  surface  layer  and  their  effects  on  permanent 
deformations.   For  the  same  reason,  generally  lower  values  of  GNIT(1) 
were  used  for  the  conventional  sections.   In  general,  rutting  was 
overpredicted  for  the  first  six  years  of  study  for  the  full-depth  and 
deep-strength  pavements  and  underpredicted  for  the  conventional  pave- 
ments. 

The  plots  of  Figure  12  were  made  after  the  four  runs  on  Section  3 
to  gain  as  much  insight  as  possible  as  to  the  effects  of  varying 
ALPHA(l)  and  GNU(l).   Figure  12  indicates  that  relatively  small  de- 
creases in  ALPHA(l)  did  not  greatly  increase  rut  depth,  except  when 

accompanied  by  substantial  increases  in  GNU(l).   As  the  greatest  por- 
tion of  the  rutting  occurs  in  the  summertime,  the  values  of  ALPHA(l) 

and  GNU(l)  for  the  summer  season  and  the  four  solutions  for  Section 
3  are  plotted  in  Figure  13.   Comparison  of  Figures  12  and  13  indicates 
that  there  was  very  little  difference  in  rutting  predictions  for  Runs 
2  and  3,  although  a  substantial  reduction  was  made  in  ALPHA(l)  for 
Run  3.   For  essentially  the  same  value  of  ALPHA(l)  but  grossly  dif- 
ferent values  of  GNU(l),  Run  1  gave  many  times  as  much  rut  depth 
predicted  as  Run  3.   Run  4  at  an  intermediate  value  of  GNU(l)  produced 
an  intermediate  value  of  rut  depth.   It  was  concluded  that  for  this 
particular  material,  variation  in  GNU(l)  was  the  most  feasible  way  of 
improving  the  permanent  deformation  characterization. 

It  should  be  borne  in  mind  that  all  of  the  entironmental  and 
mechanistic  parameters  that  control  actual  rutting  in  the  field  are 
not  modeled  in  VESYS.   It  should  also  be  noted  that  the  increase  in 
rut  depth  measured  on  the  Brampton  Sections  after  a  certain  time  is 
not  generally  representative  of  the  form  of  measured  rut  depth  curves. 
In  this  case,  rutting  increased  beginning  1969  at  a  rate  higher  than 
it  had  previously  and  this  is  very  probably  due  to  some  mechanism 
unmodeled  in  VESYS  A  or  other  mechanistic  models. 

While  the  results  gained  from  the  study  of  the  Brampton  Test 
Sections  were  useful,  no  really  strong  conclusions  could  be  reached 
on  the  basis  of  this  study  for  direct  application  to  the  four  environ- 
mental zones  being  considered. 

Analysis  of  Solutions  of  Four  Interstate  Highway  Pavement  Sections 

Information  was  available  for  development  of  input  data  for  these 
four  interstate  sections  from  a  variety  of  sources  that  included 
direct  contact  with  representatives  of  the  State  Departments  of 
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Figure  12.  Rut  depth  predicted  after  seven  years  versus  the  surface 
layer  permanent  deformation  parameters  for  Section  3, 
Brampton  Test  Road. 
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Transportation,  References  12  and  13  for  the  section  of  IH-20  in 
Illinois,  Reference  37 1  for  IH-80N  in  Utah,  Reference  382  for  IH-10 
in  Florida  and  Reference  393  for  IH-20  in  Texas.  Measured  rut  depths 
were  available  for  the  Utah  and  Florida  sections  and  were  used  on  a 
comparative  basis  as  for  previously  discussed  sections.  Although 
measured  responses  were  not  available  for  the  Texas  and  Illinois 
sections,  solutions  were  obtained  for  those  sections  and  they  were 
included  in  the  iterative  comparisons  and  input  improvements. 

Figures  14,  15,  and  16  show  measured  responses  and  those  pre- 
dicted with  the  first  solution  for  a  section  of  IH-80N  in  Utah 
between  Snowville  and  the  Idaho  border.   The  agreement  between  pre- 
dicted and  measured  responses  is  considered  excellent  for  rut  depth 
and  present  serviceability  index  and  is  quite  good  for  fatigue  crack- 
ing.  However,  there  was  apparently  no  low  temperature  cracking  noted 
whereas  the  VESYS  model  predicted  a  moderate  amount.   The  quality  of 
these  predictions  indicated  that  the  materials  characterizations  and 
other  input  were  probably  fairly  accurate,  so  their  values  were  given 
strong  consideration  in  subsequent  studies. 

Simultaneous  iterative  studies  were  conducted  for  the  pavement 
sections  in  Illinois,  Florida  and  Texas.  Some  of  these  studies  in- 
cluded : 

1.  Multiple  solutions  with  variations  in  layer  stiffnesses  and 
permanent  deformation  parameters  as  indicated  by  magnitudes 
of  predictions. 

2.  Separate  solutions  using  program  ELSYM5  in  order  to  learn 
what  the  strains  were  in  the  three  layers. 

3.  Solutions  using  ELSYM5  and  the  season  values  of  stiffness 
for  the  three  layers  in  order  to  arrive  at  deflections  at 
the  surface  for  each  of  the  different  seasons.  These  deflec- 
tions could  then  be  used  with  equations  from  VESYS  A  to  learn 
what  the  contributions  were  to  the  rutting  for  each  season. 


Anderson,  D.  I.,  J.  C.  McBride  and  D.  E.  Peterson,  "Field 
Verification  of  the  VESYS  IIM  Structural  Subsystem  in  Utah,"  Pro- 
ceedings, Fourth  International  Conference,  Structural  Design  of 
Asphalt  Pavements,  Volume  1,  August  1977. 

2 
Sharma,  J.,  L.  L.  Smith  and  B.  E.  Ruth,  "Implementation  and 

Verification  of  Flexible  Pavement  Design  Methodology,"  Proceedings, 

Fourth  International  Conference,  Structural  Design  of  Asphalt  Pavements, 

Volume  1,  August  1977. 

3 

Written  data  collected  previously  by  Texas  Transportation 

Institute  and  furnished  ARE  for  this  project. 

77 


S3 

o 

00 

1 

1 

33 

M 

U-l 

o 

C      • 

O    H 

•H    CD 

■U    T3 

O    H 

CD    O 

co  pq 

cd    o 

Xi 

u   ctj 

O  T3 

M-i    H 

CO     0) 

Xi  Xi 

4J    4-1 

a 

<u  -r) 

T)-    C 

CO 

a^ 

•u 

!H 

3    CD 

>-l    rH 

n 

■H 

2j 

T3    -H 
CD    > 

w 

U     3 

>4 

3    O 

CO    C 

cO  c/) 

CD 

B    G 

CD 

T3    0) 

C    S 

CO     4J 

cu 

"O    Xi 

at 

■U    .G 

O    cO 

•H    -U 

T3   tD 

CU 

u  a 

CU    -H 

cu 

60 
•H 


Cni)  Hidaa  imi 

78 


A 


55 

o 

00 

I 

SB 


m 

C 

r^ 

o 

•H 

■u 

a     • 

cu   »-i 

CO    CU 

T3 

cd    M 

O 

u  pq 
o 

M-l    o 

x: 

>,    CtJ 

+J    T3 

•H   M 

H 

•H     CD 

co 

^3  ,£ 

t^ 

cC    -P 

O   T3 

•H    C 

>   rt 

H 

cu   a) 

CO    .-1 

rH 

CM 

T3    -H 

r^ 

CD     > 

M    3 

3    O 

co    C 

cd  c/3 

cu 

e  C 

CU 

rH 

X)    01 

r^~ 

C    S 

cO    -UJ 

CU 

T3    43 

CU 

■u  s: 

o   cd 

•H    -u 

-a  £> 

o 

cu 

r~» 

^  a 

P«    -H 

cu 

oO 
•H 


CO 


CM 


O 


ISd 


79 


<x> 


o 


(^ 

vD 


M 

d 

0) 

o 

T) 

•H 

J-J 

4J 

0 

a 

M 

CD 

CO 

o 

-G 

cO 

CO 

td 

U 

H 

o 

M-l 

di 

>C 

0) 

4-1 

M 

cd 

X) 

E 

d 

cti 

cO 

X) 

a) 

60  t- 

c 

H 

•H 

•rl 

^ 

fe 

CJ 

3 

CO 

o 

u 

e 

o 

w 

"O 

c 

CD 

CD 

U 

01 

3 

S 

CO 

4-1 

CO 

QJ 

a) 

rO 

e 

(**| 

T3 

CO 

(3 

4-1 

rt 

1= 

TD 

c 

a) 

•rl 

—i 

CJ 

53 

•H 

O 

T3 

00 

0) 

1 

rl 

ffi 

!X 

M 

vO 


a) 

rl 

•H 

Pn 


O 

-3- 


O 
CO 


o 

CN 


O 


('3.1    "bS   0001/ '^.i    *bs)    Sut>iobj3 


80 


The  initial  solutions  or  predictions  for  these  three  sections  all 
yielded  predicted  rutting  that  was  too  high.   Analysis  indicated  that 
the  values  for  ALPHA(l)  were  in  some  cases  too  low,  that  ALPHA(3)  and 
GNU (3)  were  too  high  for  Illinois  and  that  the  creep  compliance  for 
both  the  base  and  subgrade  layers  were  too  high  for  Illinois.   It  was 
also  determined  that  too  much  effect  was  being  assigned  to  spring  break- 
up in  the  two  northern  environments  because  the  spring  season  was  three 
months  in  duration  and  the  period  for  breakup  was  much  shorter  than 
that. 

It  can  be  seen  from  Table  6  that  the  values  of  ALPHA (1)  for 
Illinois  IH-80  are  quite  a  lot  higher  than  for  the  other  pavements. 
It  is  interesting  to  consider  a  plot  of  the  various  values  of  ALPHA(l) 
versus  the  pavement  temperatures  for  each  of  the  interstate  sections 
considered.   Figure  17  shows  this  data  plotted  and  indicates  that 
most  of  the  values  fell  on  a  fairly  smooth  line,  with  the  exception 
of  all  Illinois  values  and  values  for  the  Winter  season.   The  Texas 
and  Illinois  Winter  values  appear  to  be  low  so  the  curve  for  Illinois 
was  drawn  as  indicated  without  consideration  of  the  Winter  value 
(the  values  for  ALPHA  for  Illinois  were  extrapolated  from  testing  of 
AASHO  pavement  cores  that  were  reported  in  Reference  1) .   The  higher 
values  of  ALPHA(l)  for  Illinois  may  be  due  to  the  fact  that  these 
pavements  had  undergone  years  of  traffic  prior  to  their  being  cored 
and  tested. 

Figure  18  is  indicative  of  the  type  of  information  obtained  from 
separate  layer  solutions  to  obtain  strain  levels  in  the  separate 
layers  and  subsequent  hand  calculations  to  convert  these  elastic 
strains  to  permanent  deformations  within  each  layer.   The  distribution 
of  permanent  deformation  between  the  layers  varies  considerably  be- 
tween pavement  sections  considered.   For  Texas  IH-20,  the  permanent 
deformation  was  similar  in  magnitude  for  each  layer.   For  Utah,  the 
base  material  contributed  almost  nothing  to  the  rutting.   For  Illinois, 
only  the  base  material  was  contributing  to  the  rutting.   For  Florida 
also,  the  base  material  was  contributing  most  of  the  rutting  with  the 
subgrade  and  surface  layers  making  some  contribution.   From  this  kind 
of  analysis  it  appear  obvious  that  the  input  data  for  certain  layers 
is  in  error,  but  it  is  not  at  all  obvious  from  comparisons  of  the  in- 
put data  alone.   As  a  typical  result  of  this  analysis,  it  was  obvious 
immediately  that  the  materials  characterizations  for  the  base  layer 
for  Illinois  IH-80  were  in  error  and  should  be  carefully  scrutinized. 

Figure  19  provides  data  in  terms  of  seasonal  contributions  to 
rutting  much  like  that  in  terms  of  layer  contributions  in  Figure  18. 
It  can  be  seen  that  the  rutting  in  Texas  was  predicted  primarily  to 
be  in  the  summer  and  that  the  same  was  true  for  Florida  with  more 
occurring  in  the  spring  and  fall  than  for  Texas,  consistent  with 
Florida's  higher  pavement  temperatures  in  those  seasons.   While  the 
results  for  Texas  and  Florida  appear  rational,  the  occurrence  of 
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virtually  all  of  the  rutting  in  the  springtime  for  Illinois  appears 
to  indicate  that  much  too  much  emphasis  was  being  given  to  the  spring 
breakup  in  establishing  materials  characterizations  for  the  base 
course  and  subgrade. 

As  the  study  described  with  typical  cases  above  gave  increasingly 
better  results  with  the  iterative  trials  and  comparisons,  it  became 
more  and  more  clear  what  input  values  would  be  rational  and  effective 
for  the  factorial  study. 

Final  Input  Values  for  Solution  Factorial 

The  iterative  scheme  of  predictions  (comparisons  to  measured  or 
logical  values,  analysis  of  the  comparisons  to  isolate  possible  causes 
of  differences  and  analyses  of  the  probable  input  deficiencies  to  ar- 
rive at  rational  corrections)  applied  to  input  values  described  above 
for  the  AASHO  Road  Test,  Brampton  Test  Road  and  interstate  highway 
sections  led  to  the  input  values  for  the  factorial  study  defined  in  the 
previous  chapter.   These  input  values  represented  the  best  that  could 
be  developed  by  the  authors  for  producing  rational  and  reasonable  pre- 
dictions of  pavement  distress  and  performance.   The  results  of  using 
these  input  values  are  discussed  in  subsequent  chapters. 
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CHAPTER  VI 
FACTORIAL  SOLUTIONS  FOR  INITIAL  CONSTRUCTION  AND  OVERLAYS 

Results  From  Initial  Construction  Solutions 

A  factorial  of  64  solutions,  using  input  data  developed  as 
described  in  Chapters  III,  IV,  and  V,  was  run  to  determine  the  simu- 
lated behavior  of  initial  construction.   The  solutions  were  obtained 
in  two  parts.   First   the  static  solutions  required  for  a  given  en- 
vironment and  pavement  structure  were  computed  and  saved  and  then  the 
eight  repeated  load  solutions  for  that  environment  and  structure  were 
run.   This  procedure  was  made  possible  by  the  addition  into  VESYS  A 
of  the  capability  to  write  out  the  results  of  the  static  solutions 
for  future  use;  it  proved  extremely  helpful  and  economical  in  both 
computer  and  analyst  time. 

The  results  of  these  factorial  solutions  are  presented  in  Tables 
17,  18,  19,  and  20  for  rut  depth,  PSI,  Fatigue  Damage  Index,  and  low 
temperature  cracking  index  (where  applicable),  respectively.   Values 
are  included  for  the  ends  of  years  1,  5,  10,  15  and  20.   Note  that 
this  related  only  to  the  distress  and  performance  of  the  initial 
construction;  no  overlays  are  included  in  this  set  of  tables.   Note 
also  that  the  predicted  low  temperature  cracking  index  values  are 
functions  only  of  the  pavement  structure,  age,  and  environment;  there 
is  no  dependence  on  either  the  quantity  or  load  distribution  of  traffic. 

Overlays  to  Restore  Pavements  and  Their  Simulations 

in  VESYS  A 

The  following  overlay  strategy,  based  on  common  practice,  was 
chosen:   For  pavements  which  failed  in  cracking  (prediction  of  total 
square  yards  cracked/1000  sq. yards  >  999.0),  a  three-inch  (7.6  cm) 
overlay  was  added  for  the  low  traffic  cases  and  a  four-inch  (10.2  cm) 
overlay  was  added  for  the  high  traffic  cases.   Where  a  pavement  failed 
in  rutting  (PSI  <  2.5,  found  only  in  the  wet-no  freeze  climate),  a 
two-inch  (5.1  cm)  overlay  was  added  for  both  high  and  low  traffic. 

As  we  now  have  a  four-layer  system  and  VESYS  is  only  capable  of 
analyzing  a  three-layer  system,  it  was  necessary  to  combine  the  two 
bituminous  layers  into  one  layer  of  some  equivalent  or  "effective" 
thickness.   For  the  rutting  failure  cases,  the  basic  structure  was 
felt  to  be  sound  (low  accumulated  fatigue)  so  the  effective  thickness 
adopted  was  simply  the  sum  of  the  overlay  and  original  pavement  thick- 
nesses.  The  creep  compliance  curve  used  for  the  new  combined  layer 
was  the  same  as  for  the  original  asphalt. 
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For  those  cases  where  the  original  pavement  failed  in  cracking, 
a  separate  layered-theory  program  (LAYER)  capable  of  handling  up  to 
five  layers  was  used  to  calculate  strains  under  the  top  layer  in  two 
different  layered  systems.   The  first  system  consisted  of  base  and 
subgrade  as  used  in  the  VESYS  solutions  (moduli  obtained  by  inverting 
the  constant  values  of  creep  compliance  used  for  these  layers),  the 
original  pavement  with  a  reduced  "cracked"  modulus  of  70,000  psi 
(483,000  kN/m2) (adopted  for  Class  2   cracked  asphalt  in  Reference  4Q1), 
and  an  overlay  of  either  3  or  4  inches  (7.6  to  10.2  cm)  of  asphalt  with 
a  modulus  of  752,000  psi  (52,061  kg/cm2)  derived  from  the  creep  compli- 
ance curve  at  a  load  duration  of  .03  seconds.   The  second  system  had 
the  same  base  and  subgrade,  but  had  a  single  bituminous  layer  with  the 
same  modulus  as  that  adopted  for  the  overlay  and  a  thickness  varying 
from  2.5  to  6  inches  (6.3  to  15.2  cm).   For  all  of  these  solutions, 
a  Poisson's  ratio  of  0.5  was  used  for  all  layers,  as  is  used  in  VESYS. 

Solutions  of  both  kinds  were  made  for  both  the  thin  and  the  thick 
original  pavement  structures  and  for  both  3-inch  (7.6  cm)  and  4-inch 
(10.2  cm)  overlays.   The  radial  strains  calculated  for  the  three-layer 
system  are  plotted  in  Figure  20  versus  the  four  levels  of  surface  layer 
thickness  considered  (See  solid  lines).   The  radial  strains  calculated 
for  the  four-layer  solutions  were  plotted  as  horizontal  dashed  lines 
and  their  intersections  with  the  solid  lines  represent  the  equivalent 
single-layer  thickness  for  use  in  three-layer  solutions  to  produce 
radial  strain  levels  expected  in  the  four-layer  system.   The  equivalent 
thicknesses  are  given  in  Table  21. 

Table  21 

Equivalent  Thicknesses  for  Single  Layer 
Approximation  to  Overlay  +  Original  Pavement, 
Inches  (1  cm  =  0.394  inches) 

Overlay                Original  Pavement  Structure 
Thickness  Thin Thick 

3  3.9  4.1 

4  4.8  4.9 

In  the  factorial  runs  on  the  original  pavements,  there  were  no 
cases  where  an  overlay  became  necessary  on  a  thick  pavement  under  low 


Austin  Research  Engineers  Inc,  "Asphalt  Concrete  Overlays  of 
Flexible  Pavements,  Volume  1  -  Development  of  New  Design  Criteria, 
FHWA  Report  No.  FHWA-RD-75-75,  August  1975. 
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traffic,  corresponding  to  a  3-inch  (7.6  cm)  overlay  requirement,  so 
solutions  with  an  equivalent  thickness  of  4.1  inches  (10.4  cm)  were 
not  required. 

New  static  solutions  were  performed  for  the  climate/pavement 
system  combinations  and  top  layer  thicknesses  listed  in  Table  22. 
As  the  equivalent  thickness  of  3.9  inches  (9.9  cm)  is  identical  to 
one  thickness  of  original  pavement,  existing  static  solutions  could 
be  used  for  the  thin  pavement  systems  for  wet-freeze,  dry-freeze  and 
dry-no  freeze  environments  where  the  low  traffic  levels  were  to  be 
simulated. 


Table  22.   New  Static 

Solutions  for 

Overlay 

s  (1  ( 

zm   =  0.3 

Climate 

Pavement  System 

AC  Thickness 

Wet-Freeze 

Thin 

4.8 

inches 

Wet-Freeze 

Thick 

4.9 

inches 

Dry-Freeze 

Thin 

4.8 

inches 

Dry-Freeze 

Thick 

4.9 

inches 

Wet-No  Freeze 

Thin 

5.9 

inches 

Wet-No  Freeze 

Thick 

6.5 

inches 

Dry-No  Freeze 

Thin 

4.8 

inches 

Dry-No  Freeze 

Thick 

4.9 

inches 

With  this  array  of  static  solutions  now  available,  simulations 
were  run  using  values  of  average  daily  traffic  (LAMBDA)  picking  up 
at  the  year  following  the  year  in  which  the  original  pavement  failed. 
All  other  data  remained  the  same.   Where  it  was  necessary  to  extend 
the  simulations  beyond  the  point  (20  years)  to  which  traffic  data  had 
been  given,  a  linear  extrapolation  was  used,  adding  100  to  the  value 
of  LAMBDA  for  each  year  for  the  low  traffic,  and  adding  150  each  year 
for  high  traffic.   These  extensions  were  necessary  to  obtain  the  es- 
timated lifetime  of  the  original  pavement  with  any  overlays  when  it 
exceeded  20  years. 

In  the  cases  where  more  than  one  overlay  was  required,  the  cost 
of  the  second  overlay  was  simply  added  at  the  appropriate  year  and  the 
lifetime  estimated.   The  equivalent  thickness  approach  was  not  applied 
for  a  second  overlay  because  the  extrapolation  of  this  concept  to  a 
second  overlay  did  not  appear  justified,  and  because  the  estimated  life 
based  on  established  trend  data  for  the  initial  pavement  and  one  overlay 
offered  a  basis  for  a  reasonably  accurate  estimate  of  pavement  life  for 
the  second  overlay. 
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Results  From  Overlay  Construction  Solutions 

Repeated  load  solutions  were  obtained  for  the  47  individual 
sets  of  conditions  requiring  at  least  one  overlay.   Since  these 
solutions  start  at  many  different  years  with  different  traffic  values 
over  their  lifetime,  tabulation  of  rut  depth,  PSI,  and  other  per- 
formance indicators  for  side-by-side  comparison  is  not  appropriate. 
The  performance  of  the  individual  overlays  is  indicated  on  the  graphi- 
cal displays  for  each  case,  Figures  21  through  36. 

Table  23  gives  for  all  64  individual  sets  of  conditions  the  age 
at  first  overlay  (if  any),  the  age  at  second  overlay  (if  any),  and  the 
total  life  of  the  pavement  system  as  overlaid  to  achieve  at  least  a 
20-year  lifetime.   All  of  the  total  lifetimes  in  cases  where  a  second 
overlay  was  required  are  marked  as  estimates  since,  as  mentioned  in 
the  foregoing  section,  no  attempt  was  made  to  simulate  a  second  overlay 
with  the  VESYS  program. 

It  should  be  noted  that,  as  for  the  original  construction,  all  the 
overlay  failures  were  in  cracking  except  for  the  wet-no  freeze  environ- 
ment, where  PSI  <   2.5  occurred  usually  long  before  cracking  became 
significant. 

Overall  Performance  Results 

The  predicted  performance  for  all  16  cases  including  four  levels 
of  environmental  zone,  two  levels  of  truck  traffic  volume  and  two 
levels  of  pavement  section  thickness  are  plotted  separately  in 
Figures  21  through  36.   Each  figure  includes  a  plot  of  cracking  versus 
time,  a  plot  of  rut  depth  versus  time  and  a  plot  of  serviceability 
index  versus  time.   Each  plot  includes  separate  graphs  for  each  of  the 
four  axle  load  limits.   While  there  is  value  in  separate  study  of 
these  plots,  the  major  results  and  trends  noted  from  such  a  study 
follow: 

1.  The  expected  increases  in  cracking  and  rutting  with  in- 
creased traffic  and  increased  axle  loads  were  predicted. 

2.  The  expected  important  decreases  in  cracking  and  rutting 
with  a  thicker  pavement  section  also  were  predicted. 
These  can  be  readily  seen  by  comparing  the  performance 
plots  for  the  same  environment  and  traffic  level,  but 
with  differing  pavement  section  thicknesses.   For  ex- 
ample Figure  21  indicates  a  need  for  two  overlays  to 
attain  20  years  service  for  low  traffic  with  a  legal 
axle  limit  of  24-kips  (10,900  kg)  and  a  thin  pavement, 
but  Figure  25  indicates  only  one  for  the  thick  pavement. 
Figure  31  shows  a  need  for  two  overlays  to  restore  ride 
quality  for  high  traffic  with  legal  limits  of  22  and  24 
kips  (9,980  and  10,900  kg)  and  a  thin  pavement  section 
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in  the  wet-no  freeze  zone,  while  only  one  overlay  is  re- 
required  for  the  thicker  pavement  as  shown  in  Figure  35. 

3.  Predicted  failures  were  due  to  cracking  in  all  the  environ- 
mental zones  except  the  hotter  wet-no  freeze  environment 
represented  by  Florida.   Overlays  for  the  wet-no  freeze 
environment  were  necessitated  by  excessive  rutting  and 
consequent  loss  of  serviceability. 

4.  Predicted  cracking  failure  occurred  much  more  rapidly  in  the 
colder  wet-freeze  and  dry-freeze  climates  than  in  the  dry- 
no  freeze  climate  represented  by  Texas.   Low  temperature 
cracking  also  added  to  the  higher  levels  of  the  fatigue 
cracking  in  the  colder  climates.   The  higher  temperatures  in 
the  wet-no  freeze  zones  reduced  the  incidence  of  cracking 
such  that  vertually  no  cracking  was  predicted,  except  for  a 
minor  amount  for  the  22  and  24-kip  (9,980  and  10,900  kg) 
axle  loads  for  high  traffic  on  thin  pavements. 

5.  Although  the  expected  increase  in  incidence  of  higher  axle 
loads  for  higher  legal  axle  limits  is  numerically  minor,  it 
can  be  seen  that  failures  are  predicted  considerably  earlier 
for  the  higher  legal  axle  limits  than  for  the  more  usual  18- 
kip  (8,160  kg)  limit.   This  is,  of  course,  due  to  the  in- 
creased rate  of  damage  caused  by  the  higher  axle  loads. 

6.  PSI  remained  higher  where  cracking  was  the  cause  of  overlay 
and  overlays  were  spaced  close  together  in  time.   More 
time  between  the  overlays  allowed  more  rut  development  and 
loss  of  serviceability.   (Changes  in  PSI  due  to  cracking  as 
predicted  by  VESYS  A  or  any  other  mechanistic  model  in  use 
today  are  very  small.   However,  these  changes  are  generally 
underestimated  because  moisture  entering  through  the  cracks 
will  generally  cause  loss  of  stiffness  in  underlying  layers 
that  will  result  in  an  increasing  rate  of  rutting.) 

7.  Twice  as  many  cases  (8  cases  to  4)  with  legal  limits  of  22 
and  24  kips  (9,980  and  10,900  kg)  required  two  overlays  than 
their  corresponding  18-kip  (8,160  kg)  cases.   Also,  five 
cases  of  18-kip  (8,160  kg)  legal  limits  required  no  overlays 
whereas  only  three  at  24  kips  (10,900  kg)  survived  without 
an  overlay. 
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Figure  21.   Predicted  performance,  wet-freeze  environment,  low 
traffic  and  thin  pavement  sections. 
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Figure  22-   Predicted  performance,  dry-freeze  environment, 
low  traffic  and  thin  pavement  section. 
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Figure  23.   Predicted  performance,  wet-no  freeze  environment,  low 
traffic  and  thin  pavement  section. 
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Figure  24.   Predicted  performance,  dry-no  freeze  environment, 
low  traffic,  thin  pavement  section. 
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Figure  25.  Predicted  performance,  wet-freeze  environment, 
low  traffic  and  thick  pavement  section. 
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Figure  26.   Predicted  performance, dry-freeze  environment,  low 
traffic  and  thick  pavement  section. 
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Figure  27.   Predicted  performance,  wet-no  freeze  environment, 
low  traffic,  thick  pavement  section. 
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Figure  28.  Predicted  performance,  dry-no  freeze  environment,  low 
traffic  and  thick  pavement  section. 

104 


4      6      8      10     12      14     16     18 

Time  in  Years  Since  Initial  Construction 


20 


06f- 


a.      Cracking  vs.    time 


1   cm  =   0.394   inches 
1   kg  =   0.0022   kip; 

1   sq.    yd.    =   0.8  in' 


4  6  8  10  12  14  16  18 

Time  in  Years  Since  Initial  Construction 


20 


b.   Rut  depth  vs.  time 


ai 

C 


•H 

0) 
CJ 
•H 
> 
U 

as 


6  8  10  12  14  16 

Time  in  Years  Since  Initial  Construction 


Serviceability  Index  vs.  time 


Figure  29-   Predicted  performance,  wet-freeze  environment,  high 
traffic  and  thin  pavement  section. 
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Figure  30.   Predicted  performance,  dry-freeze  environment,  high 
traffic  and  thin  pavement  section. 
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Figure  31.  Predicted  performance,  wet-no  freeze  environment, 
high  traffic  and  thin  pavement  section. 
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Figure  32.   Predicted  performance,  dry-no  freeze  environment,  high 
traffic  and  thin  pavement  section. 
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Figure  33.  Predicted  performance,  wet-freeze  environment,  high 
traffic  and  thick  pavement  section. 
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Figure  34.   Predicted  performance,  dry-freeze  environment,  high 
traffic  and  thick  pavement  section. 

110 


JOOOr 

cr 
to 

o  800 

o 
o 

r-l 

"T600 

(0 

T5 

.  400 

cr 
to 


co  200 

C 

•H 
CJ 

J-i 

o 


No  Cracking  Predicted 


1  cm  =  0.394  inches 
1  kg  =  0.0022  kips 
1  sq.  yd.  =  0.8  m2 


4      6      8      10     12      14     16 
Time  in  Years  Since  Initial  Construction 


18 


20 


a.      Cracking  vs.    time 


06f- 


4  6  8  10  12  14  16  18 

Time  in  Years  Since  Initial  Construction 


b.   Rut  depth  vs.  time 


5r 


c 


•H 

« 

O 

•H 
> 

0) 
CO 


S 


8 


10 


12 


14 


16 


18 


20 


Time  in  Years  Since  Initial  Construction 

c.   Serviceability  Index  vs.  time 

Figure  35.  Predicted  performance,  wet-no  freeze  environment, 
high  traffic  and  thick  pavement  section. 
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CHAPTER  VII 

RELATIONSHIPS  BETWEEN  LEGAL  AXLE  LOADS  AND  LIFE 
CYCLE  COSTS  FOR  FLEXIBLE  PAVEMENTS 

The  work  described  in  previous  chapters  provides  a  realistic 
data  base  in  terms  of  pavement  sections  initially  constructed  and  of 
major  maintenance  requirements  (considered  to  be  overlays  only)  to  at- 
tain the  desired  20-year  service  life.   This  chapter  deals  with  the 
final  product  of  this  research  effort,  the  relationships  between 
costs  to  attain  a  20-year  service  life  and  four  levels  of  legal  axle 
loads. 

The  data  base  on  which  the  cost  analysis  depends  consists  of 
a  factorial  of  solutions  for  original  pavement  and  overlay  lives  for 
two  pavement  thicknesses,  two  levels  of  truck  traffic,  and  four  load 
distributions  corresponding  to  four  values  of  legal  axle  load  limit. 
This  data  may  now  be  used  to  develop  cost  estimates  and  to  relate 
these  estimates  in  several  ways  to  the  legal  axle  load  limits.   It  is 
hoped  that  these  cost  comparisons  will  be  useful  in  understanding  the 
study  results  and  that  at  least  one  of  the  four  types  of  costs  consider- 
ed will  fit  the  particular  interests  or  needs  of  individual  pavement 
engineers. 

The  unit  costs  for  use  in  estimating  costs  of  initial  construction 
and  of  overlays  were  obtained  through  studies  of  Reference  41-*-  and 
correlation  with  Reference  42^.   Unit  costs  per  square  yard  of  pavement 
developed  for  this  study  appear  in  Table  24  for  initial  construction. 
The  costs  for  overlays  were  estimated  to  be  $2.39  for  2-inch  (5.1  cm) 
overlays,  $3.55  for  3-inch  (7.6  cm)  overlays  and  $4.71  for  4-inch  (10.2  cm) 
overlays.   In  order  to  limit  the  number  of  variables  in  the  study  factorial 
as  much  as  possible  without  loss  of  important  detail,  the  varying  features 
of  embankment,  shoulders,  ditches,  other  drainage  facilities,  and  other 
features  that  might  differ  between  various  environmental  zones  were  ig- 
nored.  Regional  cost  factors  were  also  ignored.   These  variables  were 
not  included  in  the  study  factorial  because  they  would  tend  to  obscure 
the  meanings  of  the  comparative  results  more  than  they  would  improve 
them  through  realism. 

It  was  originally  intended  to  carry  out  the  cost  estimates  in 
terms  of  present  worth,  but  the  unpredictable  variability  in  time  of 
interest  rates  and  of  inflation  plus  the  fact  that  interest  rates  and 


A  statistical  summary  of  current  bid  data  for  all  highway 
projects  in  the  State  of  Texas. 

2 
Means  Cost  Data,  1977. 
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inflation  rates  have  been  similar  and  off-setting  in  recent  years 
resulted  in  the  decision  to  base  the  cost  estimates  on  1977  costs 
for  both  initial  construction  and  subsequent  overlays.   This  repre- 
sents in  effect  the  present  worth  based  on  the  assumption  of  equal 
interest  and  inflation  rates.   Stated  differently,  the  probability 
of  inaccuracy  is  no  higher  through  assuming  equal  interest  and  infla- 
tion rates  than  through  attempting  to  predict  their  values  in  time. 

Cost  Evaluations  and  Relating  Costs  to  Legal  Axle  Loads 

One  of  the  very  difficult  decisions  made  during  the  research 
effort  was  how  to  meaningfully  express  the  relationships  between 
costs  and  legal  axle  loads.   One  type  of  cost  that  must  be  considered 
is  the  total  cost  to  obtain  twenty  years  of  acceptable  pavement  service, 
While  meaningful,  this  does  not  discriminate  with  much  sensitivity 
between  the  real  costs  of  operating  real  pavements  that  are  certainly 
not  to  be  abandoned  after  twenty  years.   As  can  be  seen  later,  the 
same  number  of  overlays  may  suffice  for  several  axle  load  distributions 
although  the  relative  distresses  and  performances  may  differ.   Conse- 
quently, it  was  decided  that  "value  of  remaining  pavement  life"  or 
"salvage  value"  past  twenty  years  should  also  be  included  as  a  part  of 
the  analysis  and  presentation.   It  also  was  apparent  that  costs  to  re- 
tain an  existing  pavement  in  service  under  increased  axle  load  distri- 
butions were  at  least  as  important  as  the  costs  to  build  and  maintain 
new  pavements  under  various  axle  load  distributions. 

From  the  numerous  types  of  costs  that  might  have  been  used  for 
this  comparative  study,  four  were  selected  and  are  identified  and 
defined  below: 

1.  "Total  Cost"  consisting  of  the  initial  construction 
costs  and  costs  of  necessary  overlays  to  maintain 
satisfactory  service  for  20  years. 

2.  "Average  Annual  Cost"  calculated  by  dividing  the 
"Total  Cost"  by  the  total  life  of  the  pavement. 

3.  "Total  Cost  of  Overlays"  or  the  cost  for  the  over- 
lays necessary  to  attain  20  years  of  satisfactory 
service  (same  as  "Total  Cost"  less  initial  construc- 
tion cost) . 

4.  "Average  Annual  Cost  for  Overlays"  or  the  "Total  Cost 
of  Overlays"  divided  by  the  life  of  the  pavement  after 
the  first  overlay. 

"Total  Cost"  and  "Average  Annual  Cost"  include  initial  construc- 
tion costs  and  therefore  relate  to  relative  costs  for  new  pavements 
subjected  to  truck  traffic  representative  of  different  legal  axle 
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loads,  while  the  other  two  cost  types  relate  in  an  approximate 
fashion  to  major  maintenance  costs  for  existing  pavements  subjected 
to  different  legal  axle  loads.   Total  Cost,  Average  Annual  Cost  and 
Total  Cost  of  Overlays  appear  to  be  rather  straightforward  and  their 
relation  to  real  situations  direct.   The  cost  type  labelled  "Average 
Annual  Cost  of  Overlays"  is  not  as  direct,  however,  but  does  give  a 
measure  of  sorts  for  relating  axle  loading  to  cost  of  maintaining 
pavement  serviceability.   There  are  other  cost  terms  that  can  be  de- 
fined and  considered  for  existing  pavements,  but  none  have  any  more 
general  applicability.   Both  the  "Average  Annual  Cost"  and  "Average 
Annual  Cost  of  Overlays"  bring  the  value  of  the  pavement's  remaining 
life  after  twenty  years  of  service  into  consideration. 

The  "Total  Costs"  for  the  sixty-four  combinations  of  conditions 
studied  may  be  developed  by  taking  the  appropriate  initial  cost  ($8.05 
for  a  thin  or  $10.37  for  a  thick  pavement  section)  and  adding  to  that 
the  cost  of  the  overlays.   Table  23  indicates  the  number  of  overlays. 
Their  thicknesses  and  costs  as  previously  described  were  as  follows: 

1.  Pavements  for  all  environmental  zones  except  wet-no  freeze 
received  three-inch  (7.6  cm)  overlays  at  $3.55  per  sq.  yd. 
for  low  traffic  and  four-inch  (10.2  cm)  overlays  at  $4.71 
per  sq.  yd.  for  high  traffic.   (1  sq.  yd.  =  0.8  m^) 

2.  Two-inch  (5.1  cm)  overlays  at  $2.39  per  sq.  yd.  were  placed 
on  pavements  in  the  wet-no  freeze  environmental  zone. 

The  estimated  costs  developed  in  this  manner  are  tabulated  as  "Total 
Cost"  in  Table  25  for  low  truck  traffic  and  Table  26  for  high  truck 
traffic.   The  other  three  cost  types  are  also  tabulated  in  the  same 
tables. 

With  all  the  basic  data  now  tabulated,  the  costs  relative  to 
legal  axle  loads  will  be  considered  hereafter  in  terms  of  new  pave- 
ments (initial  construction  costs  considered)  and  existing  pavements 
(only  overlay  costs  considered) . 

Cost  Relationships  for  New  Pavements 

It  can  be  readily  seen  from  review  of  Tables  25  and  26  that  the 
total  cost  was  not  dependent  on  a  legal  axle  load  for  nine  of  the 
sixteen  cases  that  were  considered  for  each  legal  axle  load.   This 
was  because  the  pavement  either  did  not  require  an  overlay  to  reach 
a  20  year  design  life  or  an  equal  number  of  overlays  were  required 
for  each  of  the  legal  axle  loads  considered.   This  does  not  imply, 
however,  that  damage  to  the  pavements  and  salvage  value  after  20  years 
were  not  affected  by  the  variations  in  axle  loads.   Average  annual  costs 
always  reflect  variations  in  cost  due  to  variations  in  total  service  life, 

Because  of  the  higher  frequency  of  cracking  in  colder  climates 
and  the  lower  frequency  of  overlays  required  where  the  failure  mode 
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is  rutting  rather  than  cracking,  the  costs  for  the  northern  environ- 
mental zones  were  higher  than  for  the  southern  environmental  zones 
for  thin  pavements.   For  the  thick  pavements,  the  differences  in 
costs  among  environmental  zones  was  minor. 

Another  useful  way  to  consider  the  data  developed  is  by  normaliz- 
ing all  results  presented  in  Tables  25  and  26  by  dividing  each  value 
by  its  corresponding  value  for  the  18-kip  (8,180  kg)  legal  axle  load. 
These  normalized  values  appear  in  Tables  27  and  28  and  indicate  rather 
directly  what  increases  in  costs  result  from  increases  in  legal  axle 
loads.   This  information  has  been  summarized  statistically  in  Table 
29  for  new  pavements  by  obtaining  mean  values  for  the  various  group- 
ings of  cases  and  environmental  zones  as  well  as  for  all  cases  combined. 

For  ease  of  observing  trends  and  identifying  some  unexpected  re- 
sults, the  average  annual  cost  (probably  more  meaningful  than  total 
cost  because  it  considers  pavement  life  after  twenty  years)  from 
Tables  25  and  26  have  been  plotted  in  Figures  37  through  40.   These 
figures  indicate  that  the  average  annual  costs  generally  increase  as 
expected  with  increasing  legal  axle  loads.   However,  an  unexpected 
result  occurs  in  some  cases  (  e. g. ,  dry-freeze  plot  in  Figure  38)  where 
cost  is  shown  decreasing  with  increasing  legal  axle  load.   This  is  not 
generally  consistent  with  expectations  from  field  experience  and  is 
partially  due  to  the  fact  that  accelerating  deterioration  as  moisture 
percolates  through  surface  cracks  into  underlying  layers  and  reflection 
cracking  are  not  modeled  by  mechanistic  models  such  as  VESYS  A. 

The  cause  of  the  decreasing  cost  (in  terms  of  average  annual  costs 
over  the  total  predicted  service  life  of  the  pavement)  with  increasing 
legal  axle  load  is  that  increasing  damage  predictions  necessitate  an 
overlay,  which  in  turn  increases  the  predicted  pavement  life  sufficiently 
that  its  effect  in  decreasing  average  annual  costs  is  greater  than  the 
increase  caused  by  the  added  overlay  cost. 

It  should  be  remembered,  however,  that  the  modeling  limitations 
mentioned  above  apply  not  only  to  overlays,  but  also  to  initial  construc- 
tion predictions  as  well.   That  is  to  say  that  pavement  life  should  al- 
ways be  overpredicted  if  the  modeling  limitation  is  serious.   As  the 
predictions  do  not  appear  unrealistic  for  either  original  pavements  or 
overlays,  it  can  be  concluded  that  the  model  limitations  may  not  be  the 
primary  cause  of  this  unexpected  result.   It  may  simply  mean  that  an 
overlay  can  not  only  meet  an  immediate  need  to  maintain  service,  but 
may  be  cost  effective.   The  Florida  Department  of  Transportation  has  in 
fact  found  it  cost  effective  to  extend  the  life  of  a  pavement  through  a 
relatively  thin  overlay  when  cracking  begins  at  the  bottom  of  the  surface 
layer  rather  than  to  apply  a  thicker  overlay  after  the  cracks  have  propa- 
gated to  the  surface. 


119 


CU 

XI 

• 

u 

CO 

CU 

>> 

CO 

.Q 

CO 

o 

CO 

CU 

o 

3 

•H 

r— i 

i+-i 

ta 

14-1 

> 

CO 

;-< 

rH 

4-1 

rH 

cO 

£ 

o 

60  r- 

C 

•H 

^ 

t: 

CO 

•h 

4-) 

> 

CO 

•H 

O 

X) 

CJ 

Po 

4-! 

,Q 

•H 

q 

T3 

•H 

0) 

H 

N 

•H 

QJ 

iH 

rH 

co 

X 

6 

cO 

£ 

O 

rH 

c 

CO 

toO 

m 

OJ 

CN 

rH 

QJ 

a 

rH 

•H 

^> 

^! 

co 

I 

H 

CO 

rH 

e 

o 

toO 

n 

C 

»4-t 

•H 

-o 

CO 

c 

4-J 

o 

(fl 

tu 

T3 

CO 

cu 

4-1 

1-4 

CO 

5-1 

O 

o 

o   o 


CN 


JO 

CO 

H 


<j- 

in 

. 

•<r 

O 

CN 

O 

O 

o 

m       o 

C 

CO 

o 

co 

rH 

o 

rH 

o 

o 

o 

co 

o 

o 

<3- 

■X 

■K 

• 

* 

* 

£3 

CO 

<M 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

T-i 

rH 

O 

H 

H 

Pu 
•H 

^ 

' 

O 

W 

o 

CO 

O 

O 

sj 

CO 

O 

<Ti 

C 

o 

o 

Vt 

o 

o 

en 

T3 

CN 

CN 

o 

CN 

o 

O 

CO 

O 

* 

# 

O 

O 

o 

o 

o 

CN 

o 

o 

H 
53 

cO 
O 

rH 

rH 

rH 

iH 

r-\ 

rH 

rH 

r-\ 

rH 

rH 

r-H 

rH 

rH 

rH 

H 

h-1 
CU 

| 

> 

t-i 

o 

CO 

o 

o 

O 

10 

o 

c 

o 

^D 

o 

o 

o 

m 

o 

o 

< 

P4 

< 

O 

o 

rH 

o 

o 

o 

rH 

o 

o 

o 

c 

o 

o 

o 

rH 

o 

o 

H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

r-\ 

rH 

rH 

rH 

rH 

iH 

CO 
toO 

cu 

H 

(-4 

o 

o 

O 

o 

o 

O 

O 

O 

o 

O 

c 

o 

c 

o 

o 

O 

00 

rH 

o 

o 

O 

O 

O 

O 

O 

O 

o 

o 

o 

O 

o 

O 

o 

O 

rH 

H 

rH 

r—i 

rH 

" 

T-i 

rH 

rH 

rH 

r-i 

rH 

rH 

rH 

rH 

rH 

T-t 

iH 

r*- 

O 

00 

rH 

m 

O 

rH 

o 

in 

O 

o 

o 

CO 

O 

CO 

<3 
CNJ 

en 

rH             O 

O 

CO 

CM 

O 

CO 

o 

rH 

c 

o 

o 

CO 

O 

CO 

* 

CO 

ex 

rH 

rH"             CN" 

rH 

rH 

rH 

CN 

rH 

rH 

rH 

rH 

r-t 

rH 

rH 

rH 

rH 

S3 

O 

r 

. 

M 

•H 

H 

W 

O 

O 

m 

O 

CO 

rH 

-vT 

o 

rH 

o 

O 

o 

O 

o 

00 

O 

CO 

w 

•> 

(M 

o 

OJ 

O 

o 

CO 

rH 

o 

CN 

o 

rH 

o 

o 

o 

CM 

O 

CO 

en 

T3 

CN 

H 
S3 
W 

CO 
O 

rH 

rH 

rH 

rH 

rH 

rH 

CN 

-1 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

T-i 



*si 

cu 

0 

rH 

> 

3 

o 

•<t 

O 

00 

T~i 

oo 

O 

o 

O 

CM 

o 

o 

O 

Sf 

O 

<y\ 

< 

O 
CN 

o 

rH 

o 

o 

CO 

o 

O 

o 

O 

O 

o 

o 

o 

rH 

O 

rH 

53 

M 

H 

cO 
toO 
cu 

1-H 

T-i 

<-i 

rH 

rH 

o 

CM 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

W 

___, 

H 

t-3 

1 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

O 

o 

O 

O 

o 

o 

I 

CO 

o 

O 

o 

o 

o 

o 

o 

o 

o 

c 

O 

o 

O 

O 

O 

o 

rH 

H 

rH 

rH 

rH 

rH 

r-\ 

rH 

r-\ 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

J 

^-^ 

CO 
S  r-H    CC 

!  <-^ 

CO 
> 

It-H  cO 

rH 

CO 

> 

rH  CO 

rH 

CO 

r*1 
rH     CO 

i 

cO 

<4-l 

CO  i— 

\  ^ 

14H 

cOrH 

CO 

CH 

COrH 

CO 

4H 

CO  rH 

H  W 

4-1 

3 

c 

O 
4-J 

3     r4 

pj  cu 

4J 

O 
4-1 

3    V4 

a  cu 

4-1 

3 

a 

o 

4J 

3   U 

a  cu 
s  > 

4-1 

3 

O 

•u 

3     r-l 

C   cu 
c  > 

<J  O 

en  P-i 

CO 

CO 

<i  o 

CO 

\< 

CO 

<0O 

CO 

< 

CO 

<iO 

CO 

< 

co 

O  >-> 

O 

o 

O  cc 

O    > 

1 

4    Cl)  U-l 

0) 

O    CO 

QJ  <J-4 

O 
CJ> 

0) 

O  CO 

CU<-M 

O 
C_> 

CU 

O  CO 

CU  MH 

toO 

Ct 

1  00  o 

toO 

CO 

toOO 

toO                CO 

toOO 

00 

cO 

toO  O 

rH 

cO 

rH  t— 

:;   cri 

rH 

CO 

rH  r-\ 

cO 

rH 

cO 

HH 

cO 

rH 

CO 

rHrH 

cO 

CO 

M4J 

CU  CO 

CO    J- 

S    AJ     Cl 

1    U  4-1 

5  cu  co 

1      4-1 

J-l  4J 
CU   CO 

CO     r4 
JJ     CU 

M  4-1 

0)   CO 

ccj 

4J 

U   4-1 

9  cu  co 

CO  U 
4-1  CU 

QJ  CO 

CO 
4-1 

U  4-J 

CU   CO 

CO  r4 
4-1  CU 

!-4    4-1 

QJ    CO 

5 

1   >  OS   o   > 

1  <:o|  H  c 

>  o 

H 

>  o 

o  > 
H  O 

5   ° 

O 
H 

^  ° 

O   > 

HO 

>  c 

O 

H 

>  o 

O  > 
HO 

>  o 
<!  u 

h-1 

<! 

w 

p- 

J 

w 

H 

CN 

^ 

3 

W 

55 

W 
W 
PS 
ft, 
1 
H 

P 

P 
p. 

3 

-< 

CN 

W 

1     W 

H  pe5 

CN 

W 

!     W 

>-<  Pi 

o  o 

> 

H 

W  Pn 

Pi  fe 
°o 

25 

M 

g 

P 

4 
3 

O 

523 

s 

P^OSiiT-.H 

, 

CO 

a. 

•H 

CN 
CN 

o 

o 


toO 


O 
U 

cu 

N 
!>> 

tO 

QJ 

•H 
> 
•H 
n3 


O 
C! 
C 
CO 

o 

CO 

CO 

QJ 
j-J 
cO 
C 
•iH 

e 

rJ 
i> 
4-1 

QJ 
T3 
C 


120 


' 

"™" '! 

<t 

rH 
CO 

CM 

o 

O 

o 

ON 

CO 

o 

o 

CO 

o 
o 

o 

<3- 

o 
o 

rH 

O 
O 

1-1  I 

o 
o 

CO 

o 
o 

CM 

in 

CM 

rH 

• 

rH 

CM 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

1-1  1 

rH 

rH 

r-i 

r-i 

z 

CO 

o 

M 

•H 

H 

^ 

o 
o 

U 

W 
en 

CM 

O 
O 

00 

O 

o 

rH 
CO 

o 
o 

CM 

o 
o 

CM 
CM 

O 
O 

CM 

rH 

rH 
rH 

O 
O 

CO 

o 
o 

ON 
CO 

H 

"3 

(T3 

CM 

rH 

o 

CM 

• 

iH 

rH 

r-i 

rH 

rH 

rH 

r-{ 

rH 

rH 

rH 

rH 

rH 

J-\ 

53 
W 

O 

0) 

rH 

a 

o 

O 

-* 

O 
O 

rH 

o 

O 

o 

o 
o 

O 

o 
o 

00 

o 

o 
o 

rH 
rH 

o 
o 

in 

rH 

o 
o 

rH 

PU 

CM 

rH 

iH 

rH 

rH 

<-i 

rH 

rH 

rH 

rH 

rH 

^ 

rH 

rH 

rH 

w 

rH 

60 

w 

(U 

H 

l-J 

00 

O 

o 

O 

o 

o       o 
o       o 

O 

o 

O 
O 

o 
o 

o 
o 

O 

o 

O 

o 
o 

O 
° 

o 
o 

O 
O 

O 

o 

O 

o 

i-H 

rH 

H 

H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

r-i 

rH 

r-\ 

rH 

<t 

O 
O 

m 

O 
O 

CO 

o 
o 

m 

CM 

o 
o 

00 
rH 

CO 
CM 

00 

o 

O 

o 

CM 

co 

CM 

rH 

o 
o 

CM 

CO 

CO 

ex 

CM 

.H 

H 

H 

rH 

rH 

rH 

rH 

l-H 

rH 

rH 

CM 

rH 

rH 

rH 

CM 

rH 

53 
o 

M 

•H 

H 

CJ 

w 

Ni 

CM 

O 
O 

rH 
CM 

o 
o 

rH 

o 
o 

ON 

rH 

o 
o 

00 
rH 

CO 
CM 

CO 
O 

o 
o 

CO 
CO 

CO 

m 
o 

O 
O 

<-*    I 
CM    1 

*    I 

C/3 

cO 

CM 

H 

H 

iH 

i-H 

^A 

rH 

r-i 

rH 

rH 

rH 

CM 

rH 

rH 

rH 

CM 

rH    § 

H 

O 

i 

53 

W 

►J 

_J 

1            1 

s 

cu 

J 

<1 

O 

O 
O 

o 

O 
O 

m 

O 

o 
o 

ON 

o 

O 

o 

00 
O 

O 

o 

O 

o 
o 

CM 

rH 

O 
O 

m 

rH 

o 

o 

m 

rH 

fa 

i-H 

CM 

i-H 

rH 

rH 

rH 

1-i 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

f~\ 

rH 

rH 

rH 

53 

tfl 

M 
PC 

00 
CI) 

H 

.J 

00 

O 
O 

O 

o 

O 
O 

O 

o 

O 

O 

o 
o 

O 

o 

o 
o 

o 
o 

o 
o 

o 
o 

O 

O 

o 
o 

O 
O 

O 
O 

O 
O 

l-\ 

H 

H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

r-\ 

rH 

rH 

rH 

rH 

rH 

CO 

CO 

CO 

CO 

r-\ 

>■ 

rH     CO 

rH 

> 

rH  CO 

i-^ 

rH    Ct 

rH 

r— 1    CO 

CO 

4-1 

CtfrH 

cO 

14-4 

CO-H 

CO 

U-l 

c0» — 

cO 

U-l 

COrH 

3 

O 

3    J-l 

3 

O 

3  u 

3 

o 

3  J- 

3 

O 

3  U 

c 

3    CU 

c 

a  cu 

C 

c  a 

C 

3  &> 

4-i 

e 

4-1 

3     > 

4-J 

3 

4-1 

a  > 

4-J 

C 

4-1 

c  > 

4-) 

3 

4-1 

3  > 

E- 
0 
C 

H     PJ 

CO 

< 

CO 

<:o 

CO 

< 

CO 

<co 

CO 

< 

CO 

<;c 

CO 

< 

CO 

n  fa 

o 
u 

0) 

O    CO 

o  > 

,    CD4H 

o 
o 

0) 

O    CO 

0)4-4 

o 
u 

o 

O  CO 

CUU- 

o 

CU 

O  CO 

u  > 

,    0)4-1 

0C 

CO 

00  o 

00 

CO 

OOO 

00 

CO 

00  c 

00 

cO 

ooo 

r- \ 

cO 

i-H  rH 

CO 

rH 

cO 

rH  rH 

cO 

rH 

cO 

i— li— 1 

cO 

rH 

cO 

1 — IrH       CD 

CO 

Jj   4-J 

CO  u 

U    4-1 

cO 

rl  4J 

CO     rJ 

J-l  4-1 

CO 

rJ     4-1 

cO  U 

r-4    4- 

CO 

J-I4J 

CO  J-l  |  U  4-J 

j_> 

Q)    CO 

4J   CU 

0J    CO 

4-1 

CU  CO 

+J     CU 

CU  CO 

4-) 

CU    CO 

4J   CU 

cu  cr 

4-1 

CU  CO 

4-J  CU  I  CU  CO 

o 

>  c 

O    > 

>  o 

o 

>  o 

O     > 

>  o 

o 

>  o 

o  > 

>  c 

O 

>  o 

o  >    >  o 

H 

!<:  cj 

HO 

<so 

H 

1  <o 

H  O 

<:o 

H 

<  CJ 

HO 

<c_ 

H 

<u 

HO!<fiO 

► 

J 

d 

3 

p 

J 

3 

[s 

3 

t- 

0 

j 

■H 

pt 

3 

[3 

J 

w 

pa 

3 

pi 

3 

p 

J 

N 

NJ 

i 

t| 

p 

/ 

c 

2 

1 

W 

1 

w 

§   S3 

2  ° 

3  N 

A 

p 

=4 

H 

w 

w 
PS 

>* 

P5 

m 
pi 

i 

E- 

H 

> 

-I 

3= 

fa 

p 

fa 

i 

i 

4 

p 

c 

4 

o 

o 

z 

53 

i 

a 

1 

CO 

a. 

•H 

CM 

CM 
O 
O 

d 
ii 
oo 

rM 

rH. 


121 


Table  29.   Mean  Values  of  Costs  Versus  Legal  Axle  Load  Normalized 
by  18-Klp  Legal  Axle  Limit  Costs. 


Data  Sample  from 
which  mean  values 


Legal  Axle  Loads,  Kips* 


are  derived 

18 

20 

22 

24 

Total  Costs: 

All  thin  pavements 

1.00 

1.04 

1.12 

1.16 

All  thick  pavements 

1.00 

1.00 

1.05 

1.12 

All  low  traffic  volumes 

1.00 

1.02 

1.09 

1.14 

All  high  traffic  volumes 

1.00 

1.04 

1.09 

1.16 

Wet-Freeze  Environment 

1.00 

1.00 

1.08 

1.12 

Dry-Freeze  Environment 

1.00 

1.00 

1.00 

1.23 

Wet-No  Freeze  Environment 

1.00 

1.08 

1.16 

1.16 

Dry-No  Freeze  Environment 

1.00 

1.00 

1.06 

1.06 

All  pavements 

1.00 

1.00 

1.09 

1.09 

Average  Annual  Costs 


All  thin  pavements 
All  thick  pavements 
All  low  traffic  volumes 
All  high  traffic  volumes 
Wet-Freeze  Environment 
Dry-Freeze  Environment 
Wet-No  Freeze  Environment 
Dry-No  Freeze  Environment 
All  pavements 


1.00 

1.08 

1.16 

1.22 

1.00 

1.12 

1.15 

1.20 

1.00 

1.10 

1.17 

1.20 

1.00 

1.10 

1.14 

1.23 

1.00 

1.12 

1.14 

1.18 

1.00 

1.07 

1.16 

1.23 

1.00 

1.01 

1.03 

1.04 

1.00 

1.05 

1.09 

1.13 

1.00 

1.10 

1.15 

1.21 

*  1  kg  =  0.0022  kips 
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Figure  37.   Average  annual  cost  by  environmental  zone  for 
low  traffic  and  thin  pavement  section. 
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Figure  38.   Average  annual  cost  by  environmental  zone  for 
low  traffic  and  thick  pavement  section. 
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Figure  39, 


Average  annual  cost  by  environmental  zone  for 
high  traffic  and  thin  pavement. 
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Figure  A0.  Average  annual  cost  by  environmental  zone  for 
high  traffic  and  thick  pavement. 
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It  should  be  noted  that  the  increase  in  costs  predicted  for  the 
increased  axle  load  distributions  (See  Table  10)  that  result  from  a 
change  in  legal  axle  limit  are  not  very  severe.   There  is  no  real 
similarity  at  all,  for  instance,  to  a  road  test  such  as  the  AASHO  Road 
Test  where  an  increase  in  test  axle  load  means  that  essentially  all 
subsequent  traffic  has  that  axle  load. 

While  it  may  not  be  appropriate  to  claim  quantitative  accuracy 
for  the  predicted  costs,  it  does  appear  reasonable  to  expect  that  all 
the  results  will  be  biased  in  the  same  way;  the  trends  would  thus  be 
reliable  and  the  costs  erring  (if  at  all)  toward  the  low  side. 

Continuing  with  analysis  of  the  cost  predictions  developed, 
the  following  general  conclusions  may  be  drawn  from  the  summary 
values  of  Table  29: 

1.  The  total  cost  for  a  20  year  service  life  may  be  expected 
to: 

a.  Increase  nominally  for  an  increase  from  18  to  20  kips 
(8,180  to  9,000  kgs)  in  legal  single  axle  load  limits. 

b.  Increase  up  to  16  percent  for  an  increase  in  legal  axle 
limit  from  18  to  22  kips  (8,180  to  10,000  kgs)  (de- 
pendent primarily  on  pavement  thickness  and  environment), 
with  a  general  expected  cost  increase  less  than  10 
percent. 

c.  Increase  from  6  to  23  percent  for  an  increase  in  legal 
axle  limit  from  18  to  24  kips  (8,180  to  10,900  kgs), 
with  an  average  increase  of  9  percent  for  all  the 
pavements. 

2.  The  average  annual  cost  (considering  the  value  of  the  pave- 
ments remaining  life  after  20  years)  may  be  expected  to: 

a.  Increase  from  one  to  12  percent  for  an  increase  in 
legal  axle  limit  from  18  to  20  kips  (8,180  to  9,000 
kgs),  with  an  average  increase  of  10  percent. 

b.  Increase  from  three  to  seventeen  percent  for  an  increase 
in  legal  axle  limit  from  18  to  22  kips  (8,180  to  10,000 
kgs),  with  an  average  increase  to  15  percent. 

c.  Increase  from  4  to  23  percent  for  an  increase  in  legal 
axle  limit  from  18  to  24  kips  (8,180  to  10,900  kgs), 
with  an  average  increase  of  21  percent. 

3.  The  overall  effect  of  increased  axle  loads  is  more  severe 
for  thin  than  for  thick  pavements  and  for  colder  than  for 
warmer  climates. 

4.  While  the  effect  on  damage  due  to  different  levels  of  truck 
traffic  volume  on  a  pavement  section  is  obvious,  the  rates 

of  increase  in  costs  for  accommodating  the  traffic  at  different 
legal  axle  limits  are  not  very  sensitive  to  levels  of  truck 
traffic. 
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Cost  Relationships  for  Existing  Pavements 

The  cost  types  previously  described  as  "Total  Costs  of  Overlays" 
and  "Average  Annual  Costs  for  Overlays"  relate  to  an  analysis  of 
effects  of  increased  legal  axle  limits  on  existing  highways  because 
the  initial  construction  costs  are  not  included  in  the  analysis. 
Table  30  provides  mean  costs  for  overlays  arrived  at  as  previously 
described  for  new  highways  except  that  initial  construction  costs 
were  omitted.   These  values  were  normalized  by  dividing  through  by 
the  corresponding  18-kip  (8,180  kg)  axle  load  costs  and  the  normalized 
values  appear  in  Table  31. 

As  would  be  expected,  the  percentage  increases  in  overlay  costs 
for  increases  in  legal  axle  limits  are  much  greater  when  a  relatively 
large  fixed  initial  construction  cost  is  not  included.   Analysis  of 
the  results  in  Tables  30  and  31  indicates  that: 

1.  The  predicted  overlay  costs  for  the  higher  legal  axle 
limits  in  the  wet-no  freeze  environmental  zone  (or  warm 
temperature)  represented  by  Florida  IH-10  are  nominal. 

2.  As  would  be  expected,  the  total  cost  for  maintaining  the 
highway  is  much  less  for  low  than  for  high  traffic 
volumes,  but  the  percentage  increase  in  total  overlay 
costs  with  increase  in  axle  loads  is  much  higher  for  the 
low  traffic  case.   For  instance,  a  mean  predicted  total 
overlay  cost  for  a  24-kip  (10,900  kg)  legal  axle  load  on 
all  low  traffic  volumes  is  2.34  times  as  much  as  for  an 
18-kip  (8,180  kg)  legal  axle  load  compared  to  only  20 
percent  more  for  all  high  traffic  volumes. 

3.  The  predicted  overlay  costs  themselves  are  much  higher  for 
northern  than  for  southern  environments,  but  the  percentage 

increase  in  overlay  cost  with  increasing  legal  axle  load 

are  similar  to  those  for  the  southern  highways. 

4.  The  cost  of  overlaying  pavements  was  logically  predicted  to 
be  greater  for  thin  than  thick  pavements,  but  the  percentage 
increase  in  overlay  cost  with  legal  increasing  axle  load  was 
greater  for  the  thick  pavements. 

5.  In  general,  total  overlay  costs  may  be  expected  to  increase 
about  6  percent  for  an  increase  in  legal  axle  limit  from 

18  to  20  kips  (8,180  to  9,000  kg),  33  percent  for  an  in- 
crease from  18  to  22  kips  (8,180  to  10,000  kg)  and  45  per- 
cent for  an  increase  from  18  to  24  kips  (8,180  to  10,900  kg). 

6.  An  increase  in  legal  axle  limit  from  18  to  20  kips  (8,180  to 
9,000  kg)  increases  the  average  annual  cost  of  overlays  about 
10  percent,  an  increase  from  18  to  22  kips  (8,180  to  10,000 
kg)  increases  it  about  30  percent  and  an  increase  from  18  to 
24  kips  (8,180  to  10,900  kg)  increases  it  about  46  percent. 
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Table  30.  Mean  Values  of  Costs  (Dollars) 
for  Overlays  versus  Legal  Axle 
Load. 


Data  Sample  from 
which  mean  values 
are  derived 


18 


Legal  Axle  Loads,  Kips* 

20        21        22 


Total  Costs  of  Overlays: 


All  thin  pavements 
All  thick  pavements 
All  low  traffic  volumes 
All  high  traffic  volumes 
Wet-Freeze  Environment 
Dry-Freeze  Environment 
Wet -No  Freeze  Environment 
Dry-No  Freeze  Environment 
All  pavements 


4.57 

5.02 

5.90 

6.35 

2.65 

2.65 

3.69 

4.13 

1.33 

1.78 

2.22 

3.11 

5.90 

5.90 

7.07 

7.07 

4.42 

4.42 

5.60 

7.37 

5.60 

6.49 

7.37 

7.37 

1.20 

1.20 

1.79 

1.79 

3.24 

3.24 

4.42 

4.42 

3.61 

3.84 

4.80 

5.24 

Average  Annual  Costs  of  Overlays: 


All  thin  pavements 
All  thick  pavements 
All  low  traffic  volumes 
All  high  traffic  volumes 
Wet-Freeze  Environment 
Dry-Freeze  Environment 
Wet-No  Freeze  Environment 
Dry-No  Freeze  Environment 
All  pavements 


0.316 

0.344 

0.380 

0.404 

0.163 

0.183 

0.240 

0.296 

0.127 

0.138 

0.186 

0.229 

0.351 

0.389 

0.434 

0.471 

0.288 

0.318 

0.340 

0.433 

0.335 

0.350 

0.470 

0.513 

0.105 

0.118 

0.130 

0.135 

0.230 

0.268 

0.300 

0.320 

0.239 

0.263 

0.310 

0.350 

*  1  kg  =  0.0022  kips 
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Table  31.  Mean  Values  of  Costs  for  Overlays  Normalized  by  18-Kip 
Axle  Load  Costs  versus  Legal  Axle  Load. 


Data  Sample  from 
which  mean  values 
are  derived 


18 


Legal  Axle  Loads,  Kips  * 

20        22       24 


Total  Costs  of  Overlays: 


All  thin  pavements  1.00 

All  thick  pavements  1.00 

All  low  traffic  volumes  1.00 

All  high  traffic  volumes  1.00 

Wet-Freeze  Environment  1.00 

Dry-Freeze  Environment  1.00 

Wet-No  Freeze  Environment  1.00 

Dry-No  Freeze  Environment  1.00 

All  pavements  1.00 


1.10 

1.29 

1.00 

1.39 

1.34 

1.67 

1.00 

1.20 

1.00 

1.27 

1.16 

1.32 

1.00 

1.49 

1.00 

1.36 

1.06 

1.33 

1.40 
1.56 
2.34 
1.20 
1.67 
1.32 
1.49 
1.36 
1.45 


Average  Annual  Costs  of  Overlays: 


All  thin  pavements 
All  thick  pavements 
All  low  traffic  volumes 
All  high  traffic  volumes 
Wet-Freeze  Environment 
Dry-Freeze  Environment 
Wet-No  Freeze  Environment 
Dry-No  Freeze  Environment 
All  pavements 


1.00 

1.09 

1.20 

1.28 

1.00 

1.12 

1.47 

1.82 

1.00 

1.09 

1.46 

1.80 

1.00 

1.11 

1.24 

1.34 

1.00 

1.10 

1.18 

1.50 

1.00 

1.05 

1.40 

1.53 

1.00 

1.12 

1.24 

1.29 

1.00 

1.17 

1.30 

1.39 

1.00 

1.10 

1.30 

1.46 

*  1  kg  =  0.0022  kips 
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It  is  hoped  that  the  distress  and  performance  predictions,  the 
predicted  costs  related  to  legal  axle  load  and  the  analysis  and  in- 
terpretation of  the  data  will  provide  valuable  insight  into  the 
effects  of  raising  legal  axle  limits  on  our  highways.   The  data  has 
been  presented  in  detail  so  that  other  more  specific  analyses  may  be 
made  to  serve  special  interests  as  required. 
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CHAPTER  VIII 
SUMMARY 

The  important  results  of  the  reasearch  effort  reported  herein  are: 

1)  Relationships  between  costs  for  twenty  years  of  acceptable 
pavement  service  and  established  legal  single  axle  limits. 
These  relationships  are  based  on  distress  and  performance 
predicted  by  Computer  Program  VESYS  A,  using  the  best  traffic, 
environment  and  materials  characterizations  possible  in  a 
factorial  study  of  sixteen  cases  with  legal  single  axle  loads 
of  18,  20,  22  and  24-kips  (8,180,  9,000,  10,000  and  10,900  kg) 
for  each. 

2)  Improvements  to  VESYS  IIM  resulting  in  an  improved  computer 
code  called  VESYS  A.  VESYS  A  includes  all  the  capabilities 
of  VESYS  IIM  plus  capabilities  for: 

a)  discretized  representation  of  axle  load  distribution, 

b)  seasonal  characterizations  of  stiffness  and  permanent 
deformation  properties  of  materials,  and 

c)  low-temperature  cracking  predictions. 

These  improvements  allow  more  accurate  modeling  of  the  flexi- 
ble pavement  structure  and  more  accurate  predictions  of  dis- 
tress and  performance. 

3)  Increased  understanding  of  the  input  parameters  for  the 
VESYS  system,  especially  those  permanent  deformation 
parameters  that  are  relatively  new  to  the  engineering 
profession.   Primary  emphasis  in  this  effort  was  on 
definition  of  trends  in  variation  of  the  permanent 
deformation  parameters  ALPHA  and  GNU  for  asphalt  concrete 
as  functions  of  stress  and  temperature  and  of  base  and 
subgrade  materials  as  functions  of  moisture  content  and 
density. 

Relationships  between  costs  for  twenty  years  of  acceptable  pavement 
service  and  established  legal  axle  loads  are  developed  and  enumerated  in 
detail  in  the  previous  chapter  for  both  new  and  existing  pavements. 
These  results  are  summarized  below: 

1.   New  Pavements: 

a.   The  total  cost  for  a  20  year  service  life  may  be  expected 
to: 

1)   Increase  nominally  for  an  increase  from  18  to  20  kips 
(8,180  to  9,000  kgs)  in  legal  single  axle  load  limit. 
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2)  Increase  up  to  16  percent  for  an  increase  in  legal  axle 
limit  from  18  to  22  kips  (8,180  to  10,000  kgs)  (dependent 
primarily  on  pavement  thickness  and  environment),  with 
an  average  expected  cost  increase  of  9  percent  for  all 
pavements. 

3)  Increase  from  6  to  23  percent  for  an  increase  in  legal 
axle  limit  from  18  to  24  kips  (8,180  to  10,900  kgs), 
with  an  average  increase  of  9  percent  for  all  the  pave- 
ments. 

b.  The  average  annual  cost  (considering  the  value  of  the  pave- 
ments remaining  life  after  20  years)  may  be  expected  to: 

1)  Increase  from  1  to  12  percent  for  an  increase  in  legal 
axle  limit  from  18  to  20  kips  (8,180  to  9,000  kgs), 
with  an  average  increase  of  10  percent. 

2)  Increase  from  3  to  17  percent  for  an  increase  in  legal 
axle  limit  from  18  to  22  kips  (8,180  to  10,000  kgs), 
with  an  average  increase  to  15  percent. 

3)  Increase  from  4  to  23  percent  for  an  increase  in  legal 
axle  limit  from  18  to  24  kips  (8,180  to  10,900  kgs), 
with  an  average  increase  of  21  percent. 

c.  The  overall  effect  of  increased  axle  loads  is  more  severe 
for  thin  than  for  thick  pavements  and  for  colder  than  for 
warmer  climates. 

2.   Existing  Pavements: 

a.  The  predicted  overlay  costs  for  the  higher  legal  axle 
limits  in  the  wet-no  freeze  environmental  zone  (or  warm 
temperature)  represented  by  Florida  IH-10  are  nominal. 

b.  As  would  be  expected,  the  total  cost  for  maintaining  the 
highway  is  much  less  for  low  than  for  high  traffic  volume 
cases,  but  the  percentage  increase  in  total  overlay  in- 
crease in  axle  loads  is  much  higher  for  the  low  traffic 
case.   For  instance,  a  mean  predicted  total  overlay  cost 
for  a  24-kip  (10,900  kg)  legal  axle  load  on  all  low  traffic 
volumes  is  2.34  times  as  much  as  for  an  18-kip  (8,180  kg) 
legal  axle  load  compared  to  only  20  percent  more  for  all 
high  traffic  volumes. 

c.  The  predicted  overlay  costs  themselves  are  much  higher  for 
northern  than  for  southern  environments,  but  the  percentage 
increases  in  overlay  costs  with  increasing  legal  axle  load 
are  similar  to  those  for  the  southern  highways. 

d.  The  cost  of  overlaying  pavements  was  logically  predicted 
to  be  greater  for  thin  than  thick  pavements,  but  the  per- 
centage increase  in  overlay  cost  with  legal  increasing  axle 
load  was  greater  for  the  thick  pavements. 
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e.  In  general,  total  overlay  costs  may  be  expected  to  in- 
crease about  6  percent  for  an  increase  in  legal  axle 
limit  from  18  to  20  kips  (8,180  to  9,000  kgs),  33  percent 
for  an  increase  from  18  to  22  kips  (8,180  to  10,000  kgs) 
and  45  percent  for  an  increase' from  18  to  24  kips  (8,180 
to  10,900  kgs). 

f.  An  increase  in  legal  axle  limit  from  18  to  20  kips  (8,180 
to  9,000  kgs)  increases  the  average  annual  cost  of  over- 
lays about  10  percent,  an  increase  from  18  to  22  kips 
(8,180  to  10,000  kgs)  increases  it  about  30  percent  and 
an  increase  from  18  to  24  kips  (8,180  to  10,900  kgs)  in- 
creases it  about  46  percent. 

It  is  clear  that  the  accuracy  of  these  cost  relationships  is  de- 
pendent on  the  ability  of  the  flexible  pavement  structural  model  VESYS 
A  to  simulate  real  pavements  and  on  the  accuracy  of  the  characterization 
of  traffic,  axle  load  distribution,  environment  and  material  input  to  the 
model.   VESYS  A  is  certainly  one  of  the  most  complete  mechanistic  models 
of  a  flexible  pavement,  but  it  shares  the  limitation  of  most,  if  not  all, 
other  such  models  in  that  it  does  not  model  reflection  cracking  or  the 
accelerating  deterioration  of  its  structural  components  as  surface  crack- 
ing progresses  and  moisture  enters  the  base  and  subgrade  layers,  nor  does 
it  consider  loss  of  stiffness  in  the  asphalt  concrete  surface  due  to 
cracking  nor  changes  in  stiffness  and  visco-elastic  properties  as  the 
asphalt  hardens  in  time.   Also,  the  distress  and  performance  predictions 
are  very  sensitive  to  materials  characteristics  that  are  difficult  to 
accurately  quantify.   Consequently,  exceptional  quantitative  accuracy 
may  not  be  claimed,  but  all  solutions  are  subject  to  the  same  biases  so 
that  the  trends  established  may  be  expected  to  be  reliable  and  the  quanti- 
tative relationships  sufficiently  accurate  to  offer  valuable  insight  until 
the  results  of  other  much  more  comprehensive  studies  are  available. 
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APPENDIX  A 

PERMANENT  DEFORMATION  PARAMETERS  ALPHA(l)  AND  GNU(l) 
FOR  ASPHALT  CONCRETE  AS  FUNCTIONS  OF  TEMPERATURE 
AND  DEVIATOR  STRESS 

The  new  capabilities  developed  in  this  project  for  VESYS  A 
include  seasonal  variations  in  the  permanent  deformation  para- 
meters ALPHA(l)  and  GNU(l)  in  order  to  more  rationally  model 
the  contributions  of  the  surface  layer  to  rutting  by  season.   A 
limited  amount  of  information  existed  as  to  realistic  values  of 
ALPHA(l)  and  GNU(l),  most  of  which  was  collected  and  reported  in 
Reference  1.   It  is  commonly  known  that  rutting  or  permanent  de- 
formation in  asphalt  concrete  is  heavily  dependent  on  temperature 
and  stress  level,  but  insufficient  information  existed  to  define 
sufficiently  the  seasonal  variation  of  these  input  quantities  if 
values  were  known  for  one  set  of  conditions.   The  purpose  of  this 
study  was  to  learn  how  ALPHA (1)  and  GNU(l)  vary  with  pavement 
temperature  and  deviator  stress.   The  resulting  information  was 
used  as  a  partial  basis  for  selected  input  values  appearing  in 
Table  6  and  applied  to  the  study  factorial  of  VESYS  A  solutions. 

Two  major  sources  of  information  were  utilized  in  the  study. 
The  first  was  a  small  factorial  of  eighteen  tests  on  a  typical 
asphalt  concrete  mix  conducted  by  Austin  Testing  Engineers  Inc. 
The  second  was  analysis  of  data  of  data  appearing  in  the  literature 
(References  1,  431,  442,  and  453) . 


McLean,  D.  B.,  and  C.  L.  Monismith,  "Estimation  of  Permanent 
Deformation  in  Asphalt  Concrete  Layers  Due  to  Repeated  Traffic 
Loading,"  Transportation  Research  Record  510,  1974. 

2 
Morris,  Jack,   The  Prediction  of  Permanent  Deformation  in 

Asphalt  Concrete  Pavements,"  The  Transport  Group,  Department  of 

Civil  Engineering,  University  of  Waterloo,  Waterloo,  Ontario, 

Canada,  September,  1973. 

3 
Kenis,  W.  J.,  and  M.  G.  Sharma,  "Prediction  of  Rutting  and 

Development  of  Test  Procedures  for  Permanent  Deformation  in  Asphalt 
Pavements,"  paper  prepared  for  presentation  in  the  symposium  en- 
titled "Deformation  Characteristics  of  Pavement  Sections"  held 
during  the  55th  Annual  Meeting  of  the  TRB,  January  1976. 
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Permanent  Deformation  Testing  for  a  Typical 
Asphalt  Concrete  Mix  by  Austin  Testing  Engineers 

Asphalt  Concrete  Mix  Description  and  Sample  Preparation 

An  AC-10  asphalt  cement  with  an  original  penetration  of  92 
and  a  penetration  at  the  time  of  mixing  of  83  was  used  with  crushed 
limestone  aggregate  graded  as  shown  in  Figure  41.   The  asphalt  ce- 
ment content  was  6  percent  and  the  resulting  specimens  had  a  rather 
high  mean  density  of  151.0  pcf  (2416  kg/m3)  and  varied  from  148.7  to 
153.9  pcf  (2379  to  2468  kg/m3). 

Aging  was  simulated  to  some  extent  by  mixing  longer  than  usual 
during  preparation  of  the  samples.   The  mix  was  heated  to  250°F  (121°C) 
and  then  compacted  by  a  motorized  gyratory-shear  molding  press  in  molds 
6  in.  (15.24  cm)  in  diameter  and  8  in.  (20.32  cm)  in  height.   The 
compacted  specimens  were  then  cured  at  room  temperature. 

As  the  diameter  of  specimens  desired  was  4  inches  (10.16  cm), 
the  specimens  were  cooled  to  40°F  (4.44°C)  to  stiffen  them  and  then 
cored  out  to  produce  specimens  4  inches  (10.16  cm)  in  diameter  and  8 
inches  (20.3  cm)  in  height.   A  special  split  tube  (firmly  secured 
around  the  specimen)  and  base  plate  were  used  to  strengthen  and  pro- 
tect the  specimens  during  coring. 

Permanent  Deformation  Testing 

The  test  factorial  included  duplicate  tests  of  specimens  at 
three  levels  of  temperature,  40,  70  and  10Q°F  (4.4,  21.1  and  37.8°C), 
and  three  levels  of  deviator  stress,  15,  30  and  45  psi  (104,  207 
and  310  kN/m^).    Specimens  were  both  brought  to  temperature  and 
tested  in  a  special  controlled-temperature  chamber. 

The  equipment  used  and  the  test  setup  was  essentially  that 
now  considered  common  for  dynamic  modulus  testing  of  asphalt  con- 
cretes, except  that  loads  were  applied  as  "square  wave"  pulses. 
Each  load  was  applied  for  0.1  second  duration  and  with  a  0.9 
second  rest  period.   Vertical  strain  was  measured  periodically 
with  two  LVDT's  mounted  in  clamps  near  the  top  and  bottom  of  the 
specimen  and  recorded  by  an  oscillographic  recorder  such  that  both 
resilient  and  permanent  strain  could  be  read. 

All  tests  were  conducted  to  100,000  cycles  of  loading.   Loading 
was  in  unconfined  compression  such  that  the  deviator  stress  was 
that  applied  vertically. 

Permanent  strains  were  plotted  versus  number  of  load  cycles 
on  log-log  paper  in  order  to  determine  ALPHA(l)  and  GNU(l)  as 
described  in  Reference  1. 
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Test  Results 

The  results  of  this  test  program  appear  in  Table  32.  The 
sample  numbers  are  coded  such  that  the  first  digit  is  the  test 
number  for  a  specific  temperature  and  the  digits  after  the  dash 
mark  give  the  test  temperature  in  °F.   These  results  have  also 
been  plotted  in  various  ways  in  Figures  42  through  47  for  ease  of 
studying  the  data  and  trends  in  variability. 

The  general  trend  information  gained  from  this  test  program 
was: 

1.  Permanent  deformation  is  extremely  temperature  dependent 
as  indicated  by  all  other  testing  in  the  literature.   For 
example,  the  permanent  strain  at  100°F  (37.8°C)  and  100,000 
load  cycles  was  around  30  times  that  for  70°F  (21.1°C)  at 

a  stress  of  30  psi  (207  kN/m2) . 

2.  GNU(l)  and  ALPHA(l)  both  decrease  with  increasing  stress. 
The  decrease  in  ALPHA (1)  dominates  to  cause  an  important 
increase  in  permanent  deformation  with  increasing  stress. 

3.  It  appears  that  the  mechanisms  affecting  permanent  defor- 
mation differ  above  and  below  some  temperature  in  the 
order  of  60  to  70°F(15.6  to  21.1°C): 

a.  Below  60°F  (15.6°C): 

1)  GNU(l)  increases  with  temperature. 

2)  ALPHA(l)  increases  slightly  with  temperature 
for  the  higher  stress  levels,  but  decreased 
at  the  lower  stress  level. 

b.  Above  60  to  70°F  (15.6  to  21.1°C): 

1)  GNU(l)  decreases  with  increasing  temperature. 

2)  ALPHA(l)  decreases  with  increasing  temperature. 

The  rate  of  increase  in  permanent  deformation  is  much 
higher  after  60  to  70°F  (15.6  to  21.1°C)  than  it  is  at 
the  lower  temperature. 

4.  The  dynamic  modulus  is  clearly  somewhat  stress  dependent. 
Above  the  60  to  70°F  (15.6  to  21.1°C)  zone,  ER  increased 
with  stress  level,  while  below  it  decreased  markedly  with 
increasing  stress  level. 

5.  The  dynamic  modulus  has  a  well-defined  negative  correla- 
tion with  temperature.   It  decreases  from  greater  than 
1,000,000  psi  (6,900,000  kN/m2)  at  low  temperatures  to 
less  than  100,000  psi  (690,000  kn/m2)  at  higher  tempera- 
tures. 
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concrete  versus  temperature. 
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Figure  45 
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Results  of  Reviewing  Other  Permanent 
Deformation  Data  in  the  Literature 

Review  of  the  permanent  deformation  test  data  in  References 
39,  40  and  41  served  to  corroborate  the  trends  noted  from  the 
ATE  test  factorial  and  previously  described.   No  additional  or 
conflicting  trends  were  noted. 

Combining  the  New  Test  Data 
With  Data  Reported  in  Reference  1 

It  is  useful  to  add  the  test  data  produced  by  the  ATE  test 
factorial  to  Figures  52,  61  and  62  of  Reference  1.   These  augmented 
figures  appear  in  this  appendix  as  Figures  48,  49  and  50,  respectively. 
As  can  be  seen,  the  additional  data  at  lower  and  higher  temperatures 
add  considerable  insight  to  the  nature  of  these  parameters  ALPHA(l) 
and  GNU(l). 
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Figure  50.   Plots  of  ALPHA(l)  and  GNU(l)  versus  temperature. 
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APPENDIX  B 

PERMANENT  DEFORMATION  PARAMETERS  ALPHA (2)  AND  GNU (2) 

FOR  GRANULAR  BASE  MATERIALS  AS  FUNCTIONS  OF 

DENSITY  AND  MOISTURE  CONTENT 


The  new  capabilities  developed  in  this  project  for  VESYS  A 
include  seasonal  variations  in  the  permanent  deformation  parameters 
ALPHA(2)  and  GNU (2)  in  order  to  take  into  account  seasonal  changes 
in  moisture  content  and  possible  changes  in  density  (such  as  might 
occur  during  spring  breakup) .  It  was  decided  that  the  limited 
testing  budget  available  could  best  be  committed  to  asphalt  concrete 
and  subgrade  testing,  so  the  results  of  literature  research  and  past 
experience  reported  in  Reference  1  were  the  basis  for  this  study. 

Sources  other  than  Reference  1  used  for  this  study  were  References 
3,  4,  461  and  472. 

Permanent  Deformation  Testing  Reported 
by  Kalcheff  and  Hicks 

This  test  regime  included  gravel  with  soil  and  stone  dust  fines, 
samples  with  varying  amounts  of  fines,  crushed  limestone  and  a  ma- 
terial called  "VA21A  Diabase"  at  95  percent  and  100  percent  of  T180 
compaction.  The  results  of  this  study  have  been  expanded  to  calculate 
values  of  ALPHA(2)  and  GNU(2)  in  some  cases  and  general  conclusions 
of  interest  to  this  research  project  drawn  and  discussed  below: 

1.  There  is  a  considerable  difference  in  permanent  strain  for 
different  materials.  Crushed  limestone  will  generally 
undergo  less  permanent  strain  than  gravel  when  the  same 
compactive  energy  is  applied  to  prepare  specimens  and  the 
same  stress  conditions  applied. 

2.  The  amount  of  permanent  strain  is  greater  when  the  fines 
are  soil  as  compared  to  fines  of  stone  dust.  It  also  in- 
creases with  percent  of  fines  passing  the  No.  200  mesh 
sieve. 


Chisolm,  E.  E.  and  F.  C.  Townsend,  "Behavioral  Characteristics 
of  Gravelly  Sand  and  Crushed  Limestone  For  Pavement  Design,"  Report 
No.  FAA-RD-75-177,  U.  S.  Army  Waterways  Experiment  Station,  September 
1976. 

2 
Raymond,  G.  P.,  "Research  oh  Rail  Support  Improvement,"  Paper 

presented  at  Session  73,  Transportation  Research  Board,  January  1977. 
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3.  The  material  called  VA21A  Diabase  underwent  much  more 
permanent  strain  than  gravel  when  compacted  to  a  similar 
density. 

4.  The  permanent  strain  for  the  VA21A  diabase  compacted  to 
146  pcf  was  approximately  five  times  that  compacted  to 
154  pcf. 

5.  Table  6  of  Reference  4  shows  nearly  5  times  as  much  per- 
manent strain  for  a  base  at  5.5  percent  moisture  content 

as  for  the  same  base  at  4.7  percent  moisture  content.   This 
is  demonstrated  by  Figure  51. 

6.  The  resilient  modulus  for  dry  crushed  limestone  may  be  ex- 
pected to  be  much  higher  than  the  same  material  after  ca- 
pillary saturation  or  inundation. 

Permanent  Deformation  Testing  Reported 
by  Chisolm  and  Townsend 

Tests  were  conducted  to  100,000  cycles  of  loading  on  a  crushed 
limestone  and  a  gravelly  sand.  General  results  of  expanding  the 
test  data  to  obtain  ALPHA(2)  and  GNU(2)  and  further  study  are: 

1.  The  resilient  modulus  for  the  crushed  limestone  was 
generally  higher  than  that  for  the  gravelly  sand. 

2.  Gravelly  sand  results: 

a.  ALPHA(2)  is  essentially  independent  of  deviator 
stress,  so  variations  in  permanent  strain  with 
variations  in  deviator  stress  relate  to  varia- 
tions in  GNU(2).   (See  Figure  52). 

b.  Permanent  deformation  increased  almost  linearly 
for  the  test  results  appearing  in  Figure  52a,  but 
are  probably  not  typical. 

3.  Crushed  limestone  results: 

a.  Figure  53b  shows  that  ALPHA(2)  was  essentially 
independent  of  lateral  stress  when  deviator 
stresses  were  held  constant. 

b.  Slight  variations  in  GNU (2)  accompanied  great 
differences  in  permanent  strain. 

c.  Increasing  lateral  stress  (or  confinement) 
greatly  decreased  permanent  strain. 

d.  The  samples  tested  at  the  very  high  deviator 
stress  of  116  psi  (800  kN/m2)  showed  decreases 
in  both  ALPHA(2)  and  GNU(2).   This  is  not  con- 
sidered typical. 

149 


I 


o 

i-H 

c 


w 

0) 
H 

a 
>^ 
u 

cfl 
o 

►J 

o 
o 
o 


+J 
CO 

c 

•H 

cd 
j-i 

■p 

C/3 


C 

aj 
C 

CD 

6 
u 

0) 


%  Moisture 


Figure  51.   Permanent  strain  versus  percent  moisture 
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4.  ALPHA(2)  appears  in  general  to  be  fairly  independent 
of  stress,  but  GNU(2)  varies  to  account  for  variations 
in  permanent  strain. 

5.  Permanent  strain  increases  with  increasing  deviator 
stress  and  decreases  with  increasing  lateral  stress. 

Permanent  Deformation  Testing 
by  G.  P.  Raymond 

Rail  ballast  is  by  virtue  of  its  usual  gradation  and  degree  of 
compaction  a  different  material  from  a  highway  base,  but  it  was 
felt  that  some  information  of  value  might  result  from  studying  the 
results  of  Raymond's  tests.   The  primary  information  gained  was 
that  ALPHA(2)  was  also  for  this  material  essentially  independent  of 
deviator  stress  and  that  GNU(2)  values  were  high  for  relatively  loose, 
open  graded  ballast  of  low  density.  Also,  permanent  strain  increased 
with  deviator  stress.  All  of  these  conclusions  are  consistent  with 
the  results  of  testing  of  highway  bases  reported  above. 

General  Trends  for  Variations  in 
ALPHA(2)  and  GNU (2) 

1.  Crushed  limestone  will  generally  have  less  permanent 
strains  than  gravel  under  the  same  compactive  energy 
and  stress  conditions. 

2.  Permanent  strains  are  greatly  increased  by  increasing 
moisture  content.   For  example,  an  increase  from  4.7 

to  5.5  percent  moisture  content  resulted  in  an  increase 
in  permanent  strain  for  one  soil  from  0.0015  inches  per 
inch  to  0.0071  inches  per  inch. 

3.  The  resilient  modulus  for  dry  crushed  stone  may  be  ex- 
pected to  be  much  higher  than  that  for  the  same  crushed 
stone  with  capillary  saturation  or  inundation. 

4.  ALPHA (2)  appears  to  be  essentially  independent  of  both 
deviator  stress  and  lateral  stress,  but  GNU(2)  increases 
with  deviator  stress,  thus  causing  an  increase  in  per- 
manent strain.   GNU (2)  decreases  with  increasing  lateral 
pressure,  thus  causing  decreased  permanent  strain. 
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APPENDIX  C 

PERMANENT  DEFORMATION  PARAMETERS  ALPHA (3)  AND  GNU (3)  FOR 
SUBGRADE  AS  FUNCTIONS  Of  MOISTURE  CONTENT,  DENSITY 
AND  DEVIATOR  STRESS 

The  new  capabilities  developed  in  this  project  for  VESYS  IIM  in- 
clude seasonal  variations  in  materials  characterizations.   While  much 
information  exists  in  the  literature  as  to  variations  in  stiffness  of 
subgrade-type  soils  with  moisture  change  and  density,  little  exists 
for  permanent  deformation  parameters.   The  purpose  of  this  study  was  to 
learn  how  the  permanent  deformation  parameters  ALPHA(3)  AND  GNU(3)  vary 
with  moisture  content  and  density.   The  resulting  information  was  used 
as  the  basis  for  seasonal  variation  of  input  values  in  predictions  of 
pavement  responses  of  axle  loads,  traffic,  pavement  designs  and  environ- 
mental conditions. 

Two  major  sources  of  information  were  utilized  in  the  study.   The 
first  was  a  small  factorial  of  eighteen  tests  on  a  typical  silty  clay 
subgrade  conducted  by  Austin  Testing  Engineers,  Inc.   The  second  was 
analysis  of  data  appearing  in  the  literature  during  the  course  of  the 
project  References  481,  492,  503  and  51^. 


J-Monismith,  C.  L.,  N.  Ogawa  and  C.  R.  Freeme,  "Permanent  Deforma- 
tion Characteristics  of  Subgrade  Soils  in  Repeated  Loading,"  draft  of 
a  prepared  paper. 

2 
Townsend,  F.  C,  and  E.  E.  Chisolm,  "Plastic  and  Resilient 

Properties  of  Heavy  Clay  Under  Repetitive  Loadings,"  Report  No. 

FAA-RD-76-107 ,  U.  S.  Army  Engineer  Waterways  Experiment  Station, 

Soils  and  Pavements  Laboratory,  November  1976. 

3 

Sharma,  J.,  L.  L.  Smith  and  B.  E.  Ruth,  "Implementation  and 

Verification  of  Flexible  Pavement  Design  Methodology,"  paper  prepared 

for  presentation  at  the  Fourth  International  Conference  on  Structural 

Design  of  Asphalt  Pavements,  October  1976. 

M.  G.  Sharma,  "Prediction  of  Rutting  and  the  Development  of 
Test  Procedures  for  Permanent  Deformation  in  Asphalt  Pavements,"  paper 
prepared  for  presentation  at  the  55th  Annual  Transportation  Research 
Board,  January  1976. 
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Permanent  Deformation  Testing  for  a  Typical  Silty 
Clay  Subgrade  Soil  by  Austin  Testing  Engineers 

Soil  Description 

The  subgrade  soil  selected  for  testing  is  a  brown  silty  clay  having 
a  liquid  limit  of  37%  and  a  Plasticity  Index  of  19%.   Its  gradation  is 
as  follows: 


Sieve  No. 

Percent  Passing 

4 

99 

10 

97 

40 

90 

200 

76 

The  classifications  for  the  soil  selected  are  CL  for  the  Unified 
Soil  Classification  System  and  A-6  (11)  for  the  AASHO  Classification 
System. 

Approach  to  Development  of  a  Test  Factorial  and  Specimen  Preparation 

It  was  desired  to  test  over  the  practical  range  of  density  and 
moisture  conditions.   The  approach  used  to  gather  data  for  consideration 
was  to  conduct  standard  moisture-density  compaction  testing  at  three 
levels  of  compactive  energy.   The  compactive  energies  selected  approxi- 
mated a  range  of  standard  to  modified  Proctor  and  included  7.16,  13.26 
and  32.5  ft. -lbs/in3  (0.38,  0.71  and  1.74  m-N/m3) .   Using  a  standard 
4-inch  mold  with  a  volume  of  1/30  cu.  ft. (.00094  m3)  and  the  10-lb.(4.54 
kg)  hammer  dropping  18  inches,  the  number  of  blows  for  each  of  the  five 
layers  was  5.5,  10.2  and  25,  respectively.   Fractions  of  blows  were 
obtained  by  one  more  blow  on  such  layers  as  necessary  to  obtain  the 
average.   (1  inch  =2.54  cm) 

The  data  from  the  compaction  tests  were  then  plotted  in  the  usual 
moisture-density  curves  with  one  curve  for  each  compactive  energy,  and 
subsequently  with  dry  density  as  a  function  of  compactive  energy  in 
blows  per  layer  with  a  curve  for  each  moisture  content.   Study  of  the 
first  set  of  curves  allowed  selection  of  feasible  ranges  of  moisture 
content  and  dry  density  to  consider  and  the  second  how  many  blows  to 
obtain  the  dry  densities  desired.   The  number  of  blows  were  then  con- 
verted to  compactive  energy  and  this  applied  to  the  split  mold  having 
2.8-inch  diameter  and  5.8-inch  height  (7.11  cm  diameter  and  14.73  cm 
height)  used  to  prepare  the  triaxial  specimens.   The  relationship 
developed  was  that  the  number  of  blows  per  layer  in  the  2.8-inch 
(7.11  cm)  specimen  mold  was  that  for  the  4-inch  (10.16  cm)  compaction 
mold  multiplied  by  0.642.   Both  were  compacted  in  five  layers.   As 
might  be  expected,  some  adjustment  in  compactive  energy  was  required 
to  obtain  the  desired  density  in  the  triaxial  specimens. 
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The  plots  for  dry  density  versus  moisture  content  appear  in  Figure 
54  and  those  for  dry  density  versus  number  of  blows  in  Figure  55.  As  can 
be  seen,  the  optimum  moisture  content  varies  from  about  14%  to  20%, 
depending  on  compactive  energy.   Similarly,  the  maximum  dry  density 
varied  from  about  100  to  112  pcf  (1600  to  1792  kg/m3) . 

The  objective  was  to  sumilate  a  typical  subgrade  for  an  interstate 
highway.   This  is  believed  to  imply  in  most  cases  the  higher  range  of 
comparative  efforts,  so  emphasis  was  placed  on  the  13.26  to  32.5 
(modified  AASHO  density)  in-lb/in3  (1.92  to  4.71  m-N/m3)  plots.   The 
region  of  densities  and  moisture  contents  considered  appears  in  Figure  54 
as  a  hatched  reactangle.   The  eight  specific  combinations  to  be  consider- 
ed in  the  test  series  are  indicated  by  the  squared  points  on  the  periph- 
ery and  the  center  of  the  rectangular  region.  The  resulting  factorial 
experiment  is  shown  in  Figure  56.   The  sample  at  22%  moisture  content 
and  a  density  of  110  pcf  (1760  kg/m3)  could  not  be  obtained  as  it  fell 
over  the  zero  air  voids  curve. 

In  order  to  obtain  the  maximum  simulation  of  field  conditions  and 
to  minimize  variability  between  specimens,  an  attempt  was  made  to  com- 
pact all  specimens  at  12%  moisture  and  to  obtain  other  moisture  contents 
through  capillarity.   It  was  not  possible,  however,  to  obtain  constant 
moisture  contents  throughout  the  sample  and  this  was  abandoned.   Instead, 
the  samples  were  dried  and  moisture  added  to  approximate  as  closely  as 
possible  the  moisture  content  desired.   The  number  of  blows  indicated  by 
Figure  55  were  applied  as  a  starting  point  and  revised  as  necessary  to 
arrive  at  specimens  having  the  desired  density  and  moisture  content. 
Duplicate  specimens  were  prepared  for  each  point  in  the  factorial. 

Test  Procedures 

The  equipment  used  and  test  setup  was  essentially  that  now  considered 
common  for  resilient  modulus  testing  for  soils.   In  order  to  obtain  maxi- 
mum uniformity  except  in  the  control  variables  moisture  content  and  dry 
density,  all  specimens  were  tested  in  a  triaxial  cell  to  100,000  cycles 
of  loading  with  a  deviator  stress  of  2.0  psi   (14  kN/rn^)  and  a  lateral 
pressure  of  3.5  psi   (24  kN/m^).   Each  load  was  applied  as  a  "square 
wave"  of  0.1  second  duration  and  with  a  0.9  second  rest  period.   Vertical 
strain  was  measured  periodically  with  two  LVDT's  mounted  in  clamps  in 
the  usual  manner  and  recorded  by  an  oscillographic  recorder  such  that 
both  resilient  and  permanent  strain  could  be  read. 

Permanent  strains  were  plotted  versus  number  of  load  cycles  on  log- 
log  paper  as  previously  described  for  the  asphalt  concrete  tests  and 
ALPHA(3)  and  GNU (3)  determined  in  the  same  manner. 

Test  Results 

The  results  of  this  test  program  appear  in  Table  33.  The  sample 


156 


112  - 


108  - 


S  104 

PL. 


•H 
CO 

c 
o 


1-1 

Q 


100  - 


96- 


92 


—  ■      1                i                i                i            — r 

Compactive                              No.    of   Blows 
Legend                        Energy,    In-Lb/ln3                       Per  Layer 

O                                          32.5                                          25.0 

X                                          13.26                                       10.2 

•                                            7.16                                          5.5 

□    Combinations  of  density  and  moisture 
content   for   factorial 

- 

1   kg/m3  =  16  pcf 

1  m-kN/m3  =   0.145   in-lb/in3 

/    X 

/    \ 

y////x 

- 

-     / 

I 

i  ///^//X/0 

i 

A 

^            - 

/ 

- 

i                         i                         I                         i                         I 

10 


14  18         22 

Moisture  Content,  % 


26 


Figure  54.  Plots  of  dry  density  versus  moisture  content 
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numbers  are  coded  such  that  the  first  two  digits  represent  the  desired 
moisture  content,  the  next  two  the  desired  dry  density  in  PCF  and 
the  letters  A  and  B  discriminate  between  the  duplicate  samples. 

As  can  be  seen,  there  is  some  unavoidable  variability  in  the  actual 
moisture  contents  and  densities  obtained  in  the  prepared  specimens,  but 
this  was  taken  into  account  by  plotting  actual  values  during  analysis. 
The  results  appear  in  the  form  of  plots  in  Figures  57,  58  and  59.  As 
can  be  seen,  there  was  considerable  variability  in  results  and  only  very 
general  trend  information  was  obtained. 

It  is  believed  that  the  variability  was  due  primarily  to: 

1.  Changes  in  length  of  the  measuring  rods  between  clamps  for 
the  LVDT's  over  the  approximately  23  hours  of  test  due  to  inadequate 
temperature  control.   The  lesson  learned  here  was  that  invar  steel 
rods  should  be  used  due  to  their  insensitivity  to  temperature  change 
and  that  close  control  of  temperature  is  likely  required  for  long-term 
repetitive-load  soil  tests  as  well  as  for  asphalt  concrete. 

2.  It  is  believed  that  pore  water  pressures  came  into  play  at 
some  point  during  the  tests  at  the  higher  moisture  content-density 
combinations  and  affected  results. 

3.  Some  natural  variability  between  samples  (despite  care  in 
preparation)  that  cause  important  variations  in  permanent  deformation 
responses  to  long-term  loading. 

The  general  trend  information  gained  was  that  for  a  silty  clay  such 
as  that  tested: 

1.  ALPHA(3)  is  fairly  constant  with  moisture  content  and  increases 
with  increasing  density.   There  is  also  some  indication  that  it  increases 
near  optimum  moisture  content. 

2.  GNU(3)  increases  with  increasing  moisture  content  and  decreases 
with  increasing  density. 

As  a  general  note,  GNU  (3)  is  dependent  on  resilient  strain,  so 
the  stress  level  used  during  testing  needs  to  be  simulative  of  that  in 
the  field. 

Permanent  Deformation  Testing  for  a  Silty  Clay 
Reported  by  Monismith,  Ogawa  and  Freeme 

Monismith  et  al  (Reference  48)  reported  test  results  for  repetitive 


161 


3    0. 


i     k     r 


▲  • 


LEGEND 


Dry  Density  of   98  PCF 
Dry  Density  of  104  PCF 
«  Dry  Density  of  110  PCF 


L 


L 


A 


10 


12 


14       16        18       20 
Moisture  Content,  % 
a.   ALPHA(3)  as  a  function  of  moisture  content. 


22 


24 


m 


■J 

< 


r 

i          i 

A  1      i      i  A  #  i 

•  - "■"       '  "~*^rr  mi.  - 

— 

■" 

- 

• 

• 

- 

- 

• 

"■ 

LEGEND 
%          %       127o   Moisture 

■  -A — —  A   17%  Moisture 

- 

- 

Hft— *—  22%  Moisture 

- 

' 

i            i 

1       1       1       1 

98 


100 


b. 


102     104     106 
Dry  Density,  PCF 
ALPHA(3)  as  a  function  of  dry  density. 


108 


110 


1  kg/m-*  =  16  pcf 


Figure  57.   Variations  in  ALPHA(3)  with  moisture  content  and  dry  density. 
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loads  on  a  silty  clay  subgrade  from  Contra  Costa  County,  California 
having  a  liquid  limit  of  35  and  plasticity  index  of  15.   Plots  of  these 
results  in  the  referenced  paper  were  used  in  this  project  to  obtain 
values  of  ALPHA(3)  and  GNU (3)  for  three  cases  and  these  are  replotted 
in  Figures  60  and  61  at  the  same  scales  as  previous  Figures  57  and  58. 
Study  of  these  test  results  led  to  the  following  conclusions: 

1.  ALPHA(3)  at  higher  stresses  did  not  appear  to  be  greatly 
affected  by  moisture  content  or  density  for  this  limited 
sampling. 

2.  While  this  limited  sampling  does  not  provide  clear  trends, 
it  may  be  roughly  concluded  that  GNU (3)  increases  with  in- 
creasing moisture  content  and  decreases  with  increasing 
density. 

3.  Lower  stresses  generally  gave  lower  permanent  strains  and 
lower  values  of  ALPHA(3)  and  GNU(3).   (Lower  values  of 
ALPHA(3)  imply  more  of  the  total  strain  is  permanent,  al- 
though total  strains  are  smaller  for  lower  stresses) . 

Permanent  Deformation  Testing  for  a  Heavy 
Clay  Reported  by  Townsend  and  Chisolm 

Townsend  and  Chisolm  (Reference  49)  reported  test  results  for 
repetitive  loads  on  a  heavy  clay  having  a  liquid  limit  of  73  and  a 
plasticity  index  of  48.   Tabulations  of  these  results  allowed  values 
of  ALPHA(3)  and  GNU(3)  to  be  obtained  for  three  batches  of  four 
specimens  each.   Moisture  contents  and  densities  varied  between  batches 
and  vertical  stress  was  varied  within  a  batch.  All  tests  were  run  at 
a  lateral  pressure  of  2  psi  (13.8  kN/m  )  within  the  triaxial  cell  and 
were  carried  to  a  maximum  of  50,000  load  cycles. 

The  results  of  these  tests  in  terms  of  ALPHA(3),  GNU(3)  and  per- 
manent strain  after  50,000  cycles  appear  in  Figures  62,  63,  and  64. 
As  the  range  of  deviator  stresses  are  much  higher  than  that  usually 
experienced  in  a  highway  subgrade,  this  data  was  extrapolated  to  a 
deviator  stress  of  3  psi  (20.7  kN/m^)  and  plots  of  this  extrapolated 
data  appear  in  Figure  62. 

The  following  was  concluded  for  the  heavy  clay  tested: 


1.   (ALPHA(3)  varies  considerably  with  moisture  content  (0.7  to 
0.9)  between  24%  and  30%  of  moisture.   The  general  trend  was 
increasing  ALPHA (3)  as  moisture  content  increased  up  to 
some  "optimum"  after  which  it  decreased  with  additional  mois- 
ture . 
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Figure  60.  Variations  in  ALPHA(3)  with  moisture  content  and  dry  density 
from  data  in  Reference  48. 


166 


CO 


O 


2.0 


1.5 


1.0 


0.5 


10 


Yd 

= 

20  psi 

= 

107.2  PCF 

w 

= 

19.3%   . 

°d 
Yd 

= 

5  psl 

= 

112.1  PCF 

w 

= 

16.7% 

ad 

YH 

= 

20  psi 

= 

112.5  PCF 

w 

= 

16.4% 

12 


14 


16 


18 


20 


22 


Moisture  Content,  % 
a.   GNU (3)  as  a  function  of  moisture  content. 


24 


2.0 


1.5 


CO 


w  1.0 


0.5 


96 


ad 

Yd 

= 

20  psi 

= 

107.2  PCF 

w 

= 

19.3%  i 

ad 
Yd 

= 

5  psi 

= 

112.1  PCF 

w 

= 

16.7%   . 

°d 
Yd 

= 

20  psi 

= 

112.5  psi 

w 

= 

16.4% 

2  _ 


JL 


J= 


si 


100 


104         108 
Dry  Density,  PCF 
b.   GNU (3)  as  a  function  of  dry  density. 


112 


Figure  61.  Variations  in  GNU (3)  with  moisture  content  and 
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Figure  62.  Variations  in  ALPHA(3),  GNU(3)  and  permanent  strain  at 
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and  dry  density  from  data  in  Reference  49  . 
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Figure  63.  Permanent  strain  at  10,000  cycles  of  load  versus  deviator  stress. 
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2.  GNU(3)  increased  with  increasing  moisture  and  decreased 
with  increasing  density. 

3.  Permanent  deformations  increased  with  moisture  content 
increases  and  decreased  as  density  increased. 

4.  ALPHA(3)  was  erratic  with  a  deviator  stress  increase  at 
medium  or  higher  optimum  moisture  content,  but  generally 
increased  with  increasing  deviator  stress. 

5.  GNU (3)  was  rather  independent  of  deviator  stress  at  lower 
moisture  content,  but  decreased  sharply  with  increasing 
deviator  stress  at  the  higher  moisture  content. 

6.  Resilient  strain  increased  with  deviator  stress  and  at 
an  increasing  rate  as  deviator  stress  increased. 

Permanent  Deformation  Testing  for  Penn  State  Clay 
Subgrade  Reported  by  Sharma,  Kenis,  Larson  and  Grambling 

Sharma  et  al  (Reference  50)  reported  test  results  for  repetitive 
loads  on  a  clay  subgrade  soil  from  the  Pennsylvania  State  University 
Test  Track.   The  data  available  was  for  two  moisture  contents.   For  a 
moisture  content  of  20.5%,  the  resilient  modulus  was  approximately 
83,000  psi  (572,410  kN/m2),  ALPHA(3)  was  0.71,  GNU(3)  was  0.0199  and 
the  permanent  strain  after  100,000  cycles  of  load  was  1.44  x  10"^  in. /in. 
For  a  moisture  content  of  23.5%,  the  resilient  modulus  was  approximately 
15,300  psi  (105,520  kN/m2),  ALPHA(3)  was  0.64,  GNU(3)  was  0.036  and  per- 
manent strain  after  100,000  cycles  of  load  was  8.40  x  10~3  in. /in. 

Based  on  this  limited  data,  the  following  apparent  trends  were 
noted: 

1.  ALPHA (3)  decreased  with  increasing  moisture  content,  but 
only  based  on  two  pieces  of  data. 

2.  GNU (3)  increased  with  increasing  moisture  content. 

3.  Permanent  deformation  at  100,000  cycles  increased 
58-fold  for  an  increase  in  moisture  from  20.5%  to  23.0%. 

Permanent  Deformation  Testing  for  Sand  Embankment 
Reported  by  Sharma,  Smith  and  Ruth 

Sharma  et  al  (Reference  51)  reported  test  results  for  100,000 
cycles  of  repetitive  loads  on  a  sand  embankment  material  from  the 
Chief land  Test  Road  Experiment  in  Florida.   The  data  available  was  for 
samples  prepared  at  four  mqisture  contents  and  three  levels  of  dry  den- 
sity. One  set  of  values  for  ALPHA(3)  and  GNU (3)  appeared  in  the  paper 
and  others  were  calculated  from  plots  of  test  results.  The  resulting 
data  are  tabulated  top  of  Page  172. 
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Permanent  Deformation  Test  Results  on 
Florida  Sand  Embankment  Material 


Moisture 

Content 

(Percent) 

Dry 

Density 

(PCF) 

109.1 

ALPHA (3) 
0.917 

GNU(3) 
0.492 

Permanent 
Strain  @  10  Cycles 
(In. /In.) 

3.4 

4.02  x  10"3 

3.4 

111.5 

0.920 

0.592 

3.84  x  10"3 

5.9 

107.8 

0.861 

0.718 

5.22  x  10"3 

12.1 

107.5 

0.913 

1.317 

7.70  x  10"3 

1  kg/m3  = 

16 

pcf 

Comparisons  may  be  made  between  the  two  specimens  at  equal  moisture 
content  of  3.4%  and  differing  dry  densities  and  between  the  two  specimens 
having  essentially  the  same  density  and  differing  moisture  contents. 
The  results  of  these  limited  comparisons  are: 

1.  ALPHA(3),  GNU(3)  and  permanent  strain  all  three 
increased  with  moisture  content. 

2.  ALPHA(3)  was  almost  independent  of  dry  density,  but 
increased  slightly  with  the  moisture  content  increase. 

3.  GNU (3)  increased  with  increasing  dry  density. 

4.  Permanent  strain  after  100,000  cycles  decreased  with 
increasing  dry  density. 

Possible  Effects  of  Pore  Water  Pressure  on  Permanent 
Deformation  Testing  of  Clay  Subgrade  Material 

Repeated  load  testing  was  conducted  by  Brown,  Lashine  and  Hyde 
(Reference  52^-)  on  a  silty  clay  (liquid  limit  of  32%  and  plasticity  index 
of  14%)  reduced  to  slurry  and  reconsolidated  at  various  overconsolida- 
tion  ratios  to  produce  different  saturated  moisture  contents  and  densities, 
The  results  of  this  testing  indicated  decreasing  permanent  strains  with 
increasing  density  as  expected  and  noted  previously. 


1  Brown,  S.  F.,  A.K.F.  Lashine  and  A.F.L.  Hyde,  "Repeated  Load 
Triaxial  Testing  of  a  Silty  Clay,"  Geotechnique  25,  No.  1,  95-114,  1975. 
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Pore  water  pressures  during  the  tests  were  measured  and  offer  some 
valuable  insight  as  to  the  effective  stresses  applicable  during  the 
test.   Some  of  these  results  follow: 

1.  Pore  water  pressure  increases  as  loading  continues  until 
some  level  is  reached  at  which  it  stabilizes. 

2.  The  pore  water  pressures  stabilized  at  around  30%  to  35% 
of  the  deviator  stresses,  so  effective  stresses  were  about 
two-thirds  of  the  deviator  stress. 

3.  Pore  pressures  stabilized  much  earlier  in  the  test  for 
denser  samples  with  lower  moisture  contents  than  for 
those  of  less  density  and  higher  moisture  content. 

These  tests  were  carried  out  to  several  million  cycles  if  they  did 
not  fail.   It  is  important  to  note  that  ALPHA(3)  and  GNU (3)  were  still 
increasing  after  a  million  cycles  of  load  if  no  failure  was  imminent. 
For  instance,  ALPHA(3)  for  one  sample  at  a  moisture  content  of  19.7%  in- 
creased from  0.83  at  100,000  cycles  to  0.867  at  1,000,000  cycles,  while 
GNU (3)  increased  from  0.24  to  0.305.   The  indication  is  that  a  set  of 
ALPHA(3)  and  GNU (3)  taken  at  100,000  cycles  as  in  current  practice  may 
be  expected  to  overpredict  long-term  permanent  deformations  (say  100,000 
or  more  cycles) . 

Brown,  et  al  also  concluded  that  resilient  strains  were  a  function 
of  applied  stress,  while  permanent  strain  was  heavily  dependent  on  stress 
history  as  well  as  applied  stress. 

While  most  subgrade  soils  are  not  fully  saturated,  it  should  be  ex- 
pected that  pore  pressures  may  develop  during  repetitive  load  tests  of 
long  duration  due  to  consolidation  of  the  specimen  with  repeated  loads. 
The  effect  may  be  to  produce  erroneous  values  of  ALPHA(3)  and  GNU(3). 

General  Conclusions  on  ALPHA (3)  and  GNU(3)  as 
Functions  of  Moisture  Content,  Density  and  Deviator  Stress 

After  studying  the  available  test  results  individually  and  collectively, 
the  following  general  conclusions  were  reached: 

1.  GNU (3)  generally  increases  as  moisture  content  increases  and 
decreases  as  density  increases.   An  apparent  sharp  decrease  in 
GNU (3)  may  occur  if  sufficient  moisture  exists  to  accept  pore 
pressure. 

2.  ALPHA(3)  may  increase  or  decrease  somewhat  with  increasing 
moisture  content,  depending  on  the  nature  of  the  soil  and 
the  amount  of  moisture  present.   The  variation  due  to  mois- 
ture appears  to  depend  on  clay  content.  ALPHA(3)  for  the 
silty  clays  seemed  to  be  almost  independent  of  moisture 
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while  ALPHA(3)  for  the  "heavier"  clay  appears  to  vary  much 
more  with  moisture  content. 

3.  Both  ALPHA(3)  and  GNU (3)  increase  together  as  moisture  in- 
creases until  ALPHA(3)  reaches  some  level  (higher  maxima  for 
higher  density) ,  and  then  ALPHA(3)  continues  to  increase  while 
GNU (3)  sharply  decreases.   This  implies  that  at  some  point  of 
moisture  increase  the  fraction  of  total  strain  that  is  permanent 
decreases  (fraction  of  resilient  strain  increases)  although  the 
magnitude  of  permanent  strain  is  increasing.   This  is  probably 
due  to  limitations  of  a  triaxial  pressure  to  simulate  resistance 
to  large  strains  (consistent  with  higher  moisture  content)  im- 
posed by  subgrade  soils  in  the  field.   Stated  differently,  the 
passive  case  for  coefficient  of  lateral  earth  pressure  likely 
applies  after  many  wheel  loadings  in  the  field  so  that  cell 
pressure  needs  to  increase  as  the  test  progresses  in  order  to 
simulate  these  conditions  accurately. 

4.  ALPHA(3)  appears  to  be  relatively  independent  of  deviator 
stress  if  the  clay  content  is  low,  but  may  vary  considerably 
as  clay  content  increases. 

5.  GNU (3)  appeared  to  be  relatively  independent  of  deviator 
stress  at  lower  moisture  contents,  but  decreased  sharply  at 
higher  moisture  contents  as  deviator  stress  increases 
(possibly  due  to  pore  water  accepting  some  of  the  stress 
increase) . 

6.  Typical  values  of  ALPHA(3)  appear  to  range  from  0.65  to  0.95. 

7.  Typical  values  of  GNU(3)  appear  to  range  from  0.01  to  0.40, 
but  can  range  much  higher  for  unusual  cases  or  materials. 
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APPENDIX  D 

EVALUATION  OF  CORRECTED  HAJEK-HAAS  MODEL  FOR  PREDICTING 
LOW  TEMPERATURE  CRACKING  FOR  FLEXIBLE  PAVEMENTS 


Mr.  E.  C.  Novak,  Jr. ,  Research  Laboratory  Section,  Testing  and 
Research  Division,  Michigan  Highway  Commission,  published  in  Febru- 
ary 1976,  a  paper  entitled  "Evaluation  of  a  Model  for  Predicting 
Transverse  Cracking  of  Flexible  Pavements."  This  paper  was  aimed 
at  comparing  the  predicted  low-temperature  cracking  derived  from 
the  Hajek-Haas  model  (References  8  and  9  )  and  actual  cracking 
measured  on  various  sections  of  highways  in  Michigan.   This  involved 
a  very  significant  effort  including:   (1)  A  survey  of  some  32  sec- 
tions of  highway  around  Michigan  including  numerous  sub-observations 
within  each  section  to  arrive  at  actual  Cracking  Indices  I  as 
described  in  Reference  9.   These  surveys  were  conducted  using  the 
procedures  described  by  Fromm  and  Phang  in  Reference  52.*  (2)  Search- 
ing out  information  on  each  of  the  sections  to  arrive  at  the  data 
required  to  estimate  the  stiffness  of  the  asphalt  cement  at  20,000 
seconds.   (3)  Assignment  of  base  temperatures  and  winter  design 
temperatures  for  the  various  sections,  calculation  of  the  predicted 
values  of  cracking  index  for  the  32  observations  and  the  many  sub- 
observations.   (4)  Statistical  analyses  to  evaluate  the  applicability 
of  the  multiple  regression  model  to  observed  cracking  on  Michigan 
pavements.   (5)  Extensive  laboratory  studies  of  the  asphalt  proper- 
ties and  the  stiffnesses  of  cores  removed  from  various  pavements. 
These  cores  included  representation  of  areas  of  high,  medium  and 
low  cracking  indices.   The  purpose  of  this  testing  was  to  see  how 
well  the  estimated  asphalt  stiffness  values  correlated  with  the 
actual  stiffness  values  of  the  cracked  pavement  as  determined  by 
laboratory  testing.   The  test  method  used  was  a  tensile-creep  test. 

The  mean  actual  cracking  indices  and  those  predicted  are  shown 
along  with  the  data  used  in  the  Hajek-Haas  equation  in  Table  34. 
The  majority  of  Mr.  Novak's  summary  and  conclusions  is  quoted 
below  to  provide  the  detailed  results  of  this  study: 

Summary  and  Conclusions 

"In  this  study  the  Hajek-Haas  model,  designed  to  predict 
transverse  cracking  susceptibility  of  flexible  pavements,  was 
statistically  evaluated  by  testing  its  ability  to  predict  the 
mean  cracking  index  of  selected  Michigan  pavements.   A  supple- 
mental study  was  conducted  when  it  was  found  that  the  model  was 
not  suitable,  in  an  effort  to  determine  why  the  model  failed 
and  to  provide  direction  for  possible  future  work  in  this  area. 
These  studies  have  produced  the  following  observations  and  con- 
clusions: 


*52  Fromm,  H.J.,  Phang,  W.A.,  "A  Study  of  Transverse  Cracking  of 
Bituminous  Pavements,"  Department  of  Transportation  and  Communications. 
Ontario,  Report  No.  RR176,  January  1972. 
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1.  For  Michigan's  flexible  pavements,  the  Hajek-Haas 
model  lacks  the  ability  to  predict,  within  reason,  transverse 
cracking  performance. 

2.  Modification  of  the  model,  based  on  data  collected 

in  Michigan,  failed  to  improve  predictive  ability  sufficiently 
to  warrant  its  use. 

3.  Of  the  model's  five  independent  variables,  only 
stiffness  modulus  was  significantly  related  to  the  cracking 
frequency  of  Michigan  flexible  pavements. 

4.  Rational  studies  of  transverse  cracking  indicate  that 
it  is  the  bituminous  concrete  stiffness  modulus,  tensile 
strength,  and  coefficient  of  expansion  properties  acting  in 
combination  with  climatic  conditions,  that  governs  an  asphalt's 
susceptibility  to  transverse  or  thermal  cracking. 

5.  The  Hajek-Haas  model  demonstrated  poor  predictive 
ability  because  it  included  only  one  of  several  bituminous 
concrete  properties  which  affect  transverse  cracking. 

6.  No  direct  correlation  could  be  found  between  cracking 
index  and  the  following  stiffness  modulus  characteristics: 

a.  The  modulus  of  field  aged  bituminous  concrete 
determined  by  creep  test; 

b.  The  estimated  modulus  of  field  aged  asphalt 
cement  determined  using  McLeod's  method  of  estimating 
modulus  from  penetration  viscosity  data; 

c.  The  estimated  modulus  of  field-aged  bituminous 
concrete  determined  using  Heukelom  and  Klomp's  chart  for 
converting  estimated  cement  stiffness  to  mix  stiffness 
on  the  basis  of  the  volume  concentration  of  the  aggre- 
gate, Cy. 

7.  The  stiffness  modulus  of  asphalt  cement  can  be 
lowered  to  such  as  extent  that  other  properties  of  the  bitumi- 
nous mix  which  also  influence  transverse  cracking  are  dimin- 
ished to  a  non-effective  level. 

8.  The  ability  of  very  soft  asphalts  to  override  other 
mix  properties  apparently  enable  criteria  to  be  established 
that  will  permit  designers  to  essentially  eliminate  transverse 
cracking  based  on  the  estimated  stiffness  of  the  asphalt 
cement . 
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9  .  Based  on  the  transverse  cracking  performance  of 
Michigan's  flexible  pavements,  McLeod's  method  of  selecting 
penetration  grades  of  asphalt  that  will  produce  essentially 
crack-free  pavements  is  conditionally  recommended  depending 
on  the  cost  and  availability  of  soft  asphalt  grades  and  the 
ability  to  develop  mix  designs  of  suitable  high  temperature 
stability. 

10.  The  data  collected  for  this  study  indicate  that 
stiff  grades  of  asphalt,  60  to  70  penetration  grade,  can  also 
perform  in  a  crack-free  manner.   Because  of  this  and  the 
fact  that  harder  asphalt  cements  are  desirable  from  a 
stability  standpoint,  there  is  a  need  to  develop  a  method  of 
more  accurately  assessing  the  susceptibility  of  bituminous 
mixes  to  'transverse  or  thermal  cracking. 

11.  Michigan  flexible  pavements  consist  of  three 
separately  constructed  layers  of  bituminous  concrete.   For 
Convenience,  the  total  bituminous  concrete  layer  is  usually 
considered  homogeneous  and  isotropic  for  most  design  purposes. 
However,  this  study  indicates  this  assumption  to  be  incor- 
rect since  the  rheologic  and  fracture  properties  of  each 
layer  differ  significantly.   This  increases  the  difficulty 
with  which  the  properties  of  the  bituminous  layer  can  be 
characterized  and,  of  course,  increases  the  difficulty  of 
predicting  all  types  of  surface  cracking. 

12.  The  tensile  strength  of  bituminous  mixes,  at  low 
temperature,  is  dependent  on  the  tensile  strength  of  the  aggre- 
gate as  indicated  by  the  large  percentage  of  tensile  failure 
through  the  aggregate  itself." 

Unfortunately,  the  equation  published  in  Reference  9  was 
not  correct  and  this  is  the  equation  that  was  used  for  the  Novak 
study.   This  equation  is  shown  below: 

in1   o  /Q7   in30   c(6. 7966-0. 8740t  +  1.3388a) 
10  =  2.497  x  10   x  S 

x  (7.054  x  10_3)d  x  (3.193  x  10_13)m 

,0.6026s  ,1N 

x  d (1) 


where 


I  =  the  cracking  index  -  number  of  full  plus  half  of  the 
transverse  cracks  per  500  ft  (152  m) 

s  =  stiffness  of  the  asphalt,  kg/sq  cm 
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a  =  age  of  the  pavement ,  years 

t  =  thickness  of  asphalt,  in.  (1  cm  =  0.39  in.) 

d  =  type  of  subgrade,  dimensionless  code  (Sand-5,  Loam-3, 
Clay-2) 

m  =  winter  design  temperature,  degrees  centigrade 

As  the  equation  given  in  Reference  9  did  not  provide  the  same 
results  as  the  nomographs  in  the  reference,  Dr.  Haas  was  contacted 
and  questioned  in  this  regard.  He  advised  that  there  was  an  error 
in  the  equation  as  it  appeared  in  the  reference  and  that  the  cor- 
rected equation  could  be  obtained  by  multiplying  the  winter  design 
temperature  m  by  -0.1  and  multiplying  the  stiffness  s  by  0.1.  The 
resulting  corrected  equation  is  then: 

m1   o  /o-»   m30   n  i   (6. 7966-0. 8740t  +  1.3388a) 
10  =  2.497  x  10   x  0.1  s 

•  x  (7.054  x  10"3)d  x  (3.193  x  10"13)"°*lm 

,0.06026  s  /ON 

x  d         (2) 


Using  the  corrected  equation  (2)  with  the  data  appearing  in 
Table  34 ,  new  predictions  of  Cracking  Index  were  developed  and 
compared  to  the  actual  mean  cracking  indices.   These  values  appear 
in  Table  35 ,  along  with  new  residuals  and  the  ratio  of  actual  to 
predicted  cracking  indices.  As  can  be  readily  seen,  the  corrected 
equation  generally  improves  the  predictive  capability  of  the 
Hajek-Haas  model. 

Analysis  of  the  data  in  Table  35  for  the  predictions  using  the 
corrected  equation  indicate  the  following:   (1)  The  model  predicted 
correctly  that  the  10  sections  that  did  not  crack  would  not  crack. 
It  also  predicted  that  four  of  the  nine  sections  with  one  crack  or 
less  for  500  feet  (152  m)  would  not  crack.   This  appears  to  indicate 
that  the  model  does  a  quite  reasonable  job  of  determining  whether 
important  cracking  will  or  will  not  occur.   (2)  The  corrected 
equation  like  the  published  equation  generally  over  predicts  the 
cracking  where  moderate  cracking  actually  occurs,  but  it  tends  to 
over  predict  less  than  the  published  equation.   (3)  There  were  two 
actual  cases  of  cracking  index  above  20  for  observations  31  and  32. 
The  predicted  values  for  these  cases  were  very  close  to  the  actual, 
however,  observations  33  and  34  were  virtually  identical  pavements, 
ages  and  winter  design  temperatures,  but  they  displayed  much  less 
actual  cracking.  This  seems  to  imply  that  the  model  is  tuned  to 
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pavements  in  this  category  except  there  are  other  variables  involved 
that  are  not  considered  by  the  model. 

There  appears  to  be  two  direct  approaches  for  diciding 
whether  the  Hajek-Haas  equation  is  suitable  for  use.   One  is  the 
accuracy  of  the  predictions  in  terms  of  actual  number  of  cracks 
per  500  feet  (152  m)  and  how  many  cracks  are  predicted  with  rela- 
tion to  actual  cracks  experienced.  A  summary  analysis  of  the  re- 
siduals or  difference  between  actual  and  predicted  cracking  indices 
appears  in  Table  36.  A  summary  analysis  of  cracking  index  ratios, 
actual  divided  by  predicted,  appears  in  Table  37. 

It  can  be  seen  from  Table  36  that  some  47%  of  the  crack  pre- 
dictions were  within  one  crack  or  less,  56%  were  within  two  cracks 
or  less,  63%  within  4  cracks  or  less  and  that  some  31%  predicted 
seven  or  more  cracks  in  excess  of  those  actually  occurring.   Similar 
review  of  Table  37  indicates  that  38%  of  the  predictions  were 
within  10%  in  number  of  cracks,  14%  were  within  50%,  and  56%  of 
the  predictions  differed  from  the  actual  predictions  by  more  than 
50%.  As  neither  of  these  analyses  is  complete  within  itself,  it  is 
necessary  to  look  a  little  deeper.   For  instance,  a  prediction  of 
I  =  0.2  to  an  actual  cracking  index  of  1  crack  per  500  ft.  (150  m)  gives 
a  cracking  index  ratio  of  5.0,  but  really  is  a  quite  acceptable 
prediction.  A  combined  summary  analysis  appears  in  Table  38  to 
more  clearly  identify  the  success  levels  of  the  predictions.   Review 
of  Table  38  shows  that  some  53%  of  the  predictions  missed  the  num- 
ber of  cracks  either  by  1  or  less  or  were  within  10%  of  the  total 
number  or  less.   69%  predicted  within  5  cracks  or  within  a  differ- 
ence of  50%  of  the  actual  number  of  cracks  and  31%  predicted  in 
excess  of  these  latter  limits. 

It  is  clear  that  the  corrected  Hajek-Haas  equation  does  a 
considerably  better  job  of  predicting  than  the  erroneous  published 
one.   It  is  still  obvious  that  it  is  not  an  especially  accurate 
model.  However,  considering  the  wide  variations  in  actual  cracking 
indicated  in  Table  34  for  apparent  similar  pavement  sections  and 
the  limited  basis  for  its  development,  it  does  a  remarkably  good 
job.   After  review  of  several  other  models  that  have  been  developed 
in  recent  years  and  the  overall  nature  of  the  prediction  problem,  the 
writer  has  concluded  that  this  much  more  simplistic  model  will  give 
almost,  if  not  as  good,  predictions  as  some  more  complex  models 
such  as  that  of  Shahin  and  McCullough  (References  6  and  7 ) .   The 
real  problem  is  that  low  temperature  cracking  (and  indeed  any  frac- 
ture problem  involving  tensile  stresses)  is  subject  to  many  more 
parameters  than  can  reasonably  be  included  in  a  practical  model.   In 
addition,  the  relationship  of  these  parameters  to  the  formation  of 
the  crack  is  in  general  poorly  defined.   It  is  concluded  that  the 
corrected  Hajek-Haas  model  does  a  reasonable  job  of  predicting 
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whether  important  low  temperature  cracking  will  or  will  not  occur 
and  that  it  does  a  marginal  job  of  predicting  number  of  cracks  or 
crack  spacing  if  cracking  does  occur. 
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APPENDIX  E 
STRUCTURE  OF  THE  SAVED  STATIC  SOLUTION  FILE 


The  file  created  by  VESYS  A  when  a  SAVETAPE  keyword  is 
encountered  for  a  type  1  run  is  as  follows : 

ITITLE  (20  words,  alpha  numeric) 

NTEMP  (=  the  number  of  seasons) 

NDEL  (=  the  number  of  terms  in  each  Dirichlet  series) 

NRAD  (=  the  number  of  tire  footprint  radii  considered) 

NTXNR  (NTEMP  x  NRAD) 

r     thickness  of  the  asphalt  and  base  layers,  respectively. 

(TEMP(I),  1=1,  NTEMP)   NTEMP  values  of  season  temperature 
(RADIUS(I),  1=1,  NRAD)   NRAD  values  of  radius 
(GCVAR(I,4),  1=1,  NTXNR)  a  coefficient  of  variation  for  the 

Dirichlet  series  for  vertical 
deformation  for  each  static  solution 
(GCVAR(I,  7),  1=1,  NTXNR)  the  same  as  above,  but  for  radial 

strain 
(ALF(I),  I  ==  1,  NTXNR)  |  permanent  deformation  parameters  for 
(GMU(I),  I  =  1,  NTXNR) J  each  static  solution 


HI 
H2 


•the  NDEL  Dirichlet  co- 

.  efficients  for  the  radial 
GRAD(I,J)   I  =  1  NTXNR)   J  =  1  «DEL>  J  fltrain  and  vertical  defor. 
((GN0RM(I,J),  II,  NTXNR),  J  -  1,  NDEL)  ]  ^^  responses  for  each 

static  solution. 


( 


It  can  be  seen  that  the  file  will  be  of  variable  length, 
depending  on  the  number  of  static  solutions.   This  is  helpful  in 
reducing  the  storage  charges  for  desk-resident  permanent  files, 
since  no  unused  space  is  required. 


186 


APPENDIX  F 

VARIANCE  OF  CRACKING  AS  A  FUNCTION  OF 
EVALUATION  FREQUENCY 

In  the  process  of  verifying  changes  made  for  Task  A-l,  which 
involved  allowing  seasonal  variation  of  ALPHA  and  GNU  for  the  three 
layers,  and  of  the  creep  compliance  arrays  for  layers  2  and  3,  a 
comparison  was  made  between:   (1)  a  run  on  the  old  version  with  a 
particular  data  set  (using  NTEMP  =  1);  (2)  a  run  on  the  modified 
program  with  the  same  data  set;  and  (3)  a  run  on  the  modified  pro- 
gram with  the  same  data  except  that  NTEMP  was  set  to  4,  and  four 
identical  values  of  ALPHA  and  GNU  were  entered  for  each  layer  (the 
same  as  the  single  value  previously  used) ;  the  four  "season" 
multiplier  for  creep  compliance  for  layers  2  and  3  were  set  to 
1.0,  as  had  the  single  values  for  run  (2).   The  four  values  of 
TEMP  were  likewise  equal  to  the  single  value  used  before. 

The  results  of  runs  (1)  and  (2)  agreed  exactly.  The  results 
of  runs  (2)  and  (3)  differed,  but  only  in  the  predicted  variance  of 
cracking;  the  value  from  run  (3)  was  almost  exactly  1/4  the  value 
from  run  (2) .  Examination  of  the  equations  used  to  calculate  the 
variance  of  cracking  soon  showed  that,  at  least  under  some  condi- 
tions, this  was  a  result  to  be  expected.  The  relevant  equations 
from  subroutine  CRACKS  can  be  written  as  follows: 


(We  refer  to  a  "season"  as  the  period  of  time  into  which  the 
year  is  divided  by  NTEMP.) 

TDEL  =  7200.   *365.25  *  12/NTEMP  =  number  of  seconds  per 
season. 
TDDAY  =  TDEL/3600./24.  =  365.25/NTEMP  =  number  of  days  per 
season. 
=  average  daily  traffic  (axles/day)  during  the  jth 

period  of  time  (jtn  entry  in  TRANDOM). 
=  ADTj  *  TDDAY  =  number  of  axles  in  the  ith  (cumulative) 
season  (which  occurs  during  the  jtn  entry  in  TRANDOM). 


ADT, 


EN 


EMF, 
VEMF, 


) 


complicated  functions  of  the  results  from  GRESP  and 
of  the  fatigue  relation  for  the  kth  season,  where  k 
runs  from  1  to  NTEMP. 


then 
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I 

DMG  =  I     •  EMF,  •  EN.  where  k  =  1  +  [(i-1)  modulo  NTEMP] 
i=l  and  £  is  the 

(cumulative)  season 
number  at  which  the 
estimate  is  being 
calculated 


and 


I 
i=l 


VDMGn  =   T  EN.  I  (EMF,)   +  EN,  -  VEMF,  I 
x-    .**„    x  «-    k       i       k1 


In  our  problem,  assume  that  ADT  is  held  constant.   Then  EN 
(and  hence  each  term  in  the  sum  for  DMG)  will  be  inversely  propor- 
tional to  NTEMP  (see  the  definition  of  TDDAY)  and  at  any  particular 
time  DMG  will  be  the  same  for  any  value  of  NTEMP  because  the 
number  of  terms  in  the  sum  will  be  directly  proportional  to  NTEMP. 

The  situation  for  VDMG  depends  on  which  term  within  the 
parentheses  dominates.   If  the  first,  then 

I  2 

VDMG  -     I     EN.  (EMF  )    (k  depends  on  i) 

X/     ,  -     1       R 

1=1 

and,  since  EMF  does  not  depend  on  the  amount  of  traffic  or  the 
period  of  time  involved,  the  same  argument  applies  as  for  DMG. 

In  our  situation,  however,  the  second  term  dominates: 

(VEMF  •  EN) » (EMF)2 

and  hence  we  have 

1  2 

VDMG„  -  I      (EN  )   •  VEMF   (k  depends  on  i) 

Thus  the  number  of  terms  depends  directly  on  NTEMP  but  the  size  of 
each  is  inversely  proportional  to  the  square  of  NTEMP;  hence  the 
value  of  the  sum  is  inversely  proportional  to  NTEMP.   It  is  this 
result  which  does  not  seem  grounded  in  physical  reality,  namely 
that  the  variance  of  cracking  at  a  given  time  depends  on  how  often 
it  has  been  evaluated  prior  to  that  time. 

For  our  particular  case,  we  used 
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ADT  =  8000  (constant) 

AMPLITUD  =  75. 

DURATION  =  .0125 

TEMP  =  60.   (for  both  NTEMP  =  1  and  NTEMP  -  4) 

and,  with  this  particular  data  set, 

EMF  =  .289  x  10'10 

VEMF  =  1.989  x  10~22 

EN  =  2.922  x  106 

so  that  the  two  terms  in  the  equation  for  VDMG  are 

EN  (EMF)2  =  2.44  x  10_15 
(EN)2  •  VEMF  =  1.698  x  10"9 

and  the  dominance  of  the  second  term  is  quite  complete. 
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APPENDIX  G 

GUIDE  FOR  DEVELOPING  INPUT  FOR 
LOW  TEMPERATURE  CRACKING  PREDICTIONS 

The  values  required  to  evaluate  the  regression  equation  for  low 
temperature  cracking  (Hajek-Haas  model)  are: 

1.  age  of  the  pavement,  in  years. 

2.  thickness  of  the  bituminous  layer,  in  inches. 

3.  subgrade  soil  type,  in  a  numerical  code 
(5-sand;  3-loam;  2-clay) . 

4.  winter  design  temperature,   C. 

5.  stiffness  of  original  asphalt  cement,  according  to 
McLeod's  methods,  at  the  winter  design  temperature, 
in  Kg/cm  . 

Items  1  and  2  are  provided  automatically  by  input  already  required 
by  VESYS  A.   These  are  the  values  input  under  the  TRANDOM  and  THICK1 
keywords,  respectively.   Item  3  is  a  simple  input  based  on  observation 
of  the  soil  type.   Item  4  and  5,  however,  are  very  specifically  defined, 
using  tables  and/or  charts  which  may  not  be  readily  available.   There- 
fore, these  two  parameters  are  described  in  more  detail,  including  the 
specific  steps  needed  to  go  from  readily  available  data  to  the  desired 
input . 

The  winter  design  temperature  is  defined  to  be  "the  lowest  tempera- 
ture at  or  below  which  only  1%  of  the  hourly  air  temperatures  in  January 
occur  for  the  severest  winter  during  a  ten  year  period"  (Haas,  Ref.  9). 
Where  such  detailed  records  are  not  available,  Haas  (ibid)  has  given  a 
relationship  between  air  freezing  index  (40  year  average)  and  winter  de- 
sign temperature,  covering  a  range  from  500  to  more  than  3000  degree  days 
in  freezing  index,  as  well  as  a  suggested  extrapolation  to  lower  value  of 
freezing  index.   This  is  shown  in  Figure  65. 

Freezing  index  contour  maps  for  Canada  and  for  the  United  States  are 
given  in  Figures  66  and  67,  respectively,  and  may  be  used  with  Figure 
65  to  obtain  approximate  values  of  winter  design  temperature  if  other 
more  accurate  data  are  unavailable. 

The  asphalt  stiffness  value  to  be  used  is  to  be  derived  in  a  very 
specific  manner,  and  hence  considerable  detail  is  given  here  to  that 
derivation.  The  required  input  data  to  the  procedure  are  the  viscosity  at 
279°F  (137°C)  and  the  penetration  at  77°F  (25°C)  derived  from  laboratory 
tests.  From  Figure  68,  obtain  the  penetration  index  (PI)  from  the  above 
data  by  visual  interpolation,  or  by  the  following  interpolation  forumula 
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\ 


log  L  -  log  X   . 

log  L  -  log  M     K        J 

where : 

X  =  viscosity  in  centistokes  at  275°F  (135°C)  for  the 
asphalt  cement  for  which  P.I.  is  to  be  determined. 

L  =  viscosity  in  centistokes  at  275°F  (135°C)  for  a  P.I. 
of  0.0,  determined  from  Figure  68  for  the  penetration 
at  77°F  (25°C)  for  the  given  asphalt  cement. 

M  =  viscosity  in  centistokes  at  275°F  (135°C)  for  a  P.I. 
of  -1.5,  determined  from  Figure  68  for  the  penetration 
at  77°F  (25°C)  for  the  given  asphalt  cement. 

The  ring  and  ball  temperature  (Tr&b)  or  base  temperature  is  required 
for  the  next  step  and  may  be  measured  in  the  laboratory  or  determined 
from  P.I.  and  penetration  using  the  nomograph  of  Figure  69,  which  gives 
a  temperature  difference  AT  between  base  temperature  and  the  temperature 
employed  for  the  pentration  test  (Tjen  )  in  the  sense. 

Base  T  =  T    +  AT  (nomograph) 
pen  or/ 

Determine  the  stiffness  at  a  loading  time  of  20,000  seconds 
(5.6  hours)  from  the  nomograph  of  Figure  70,  entering  the  temperature 
scale  with  the  difference  between  the  base  temperature  and  the  winter 
design  temperature  (°C) .   Note  that  at  AT  from  Figure  69,  is  in  °C,  so 
the  value  of  T    must  be  converted  to  °C  before  the  base  temperature  is 
obtained. 

The  regression  equation  is  not  meaningful  for  values  of  asphalt 
cement  stiffness  S  equal  to  or  less  than  10  kg  /cm  ,  so  a  value  of 
zero  would  be  input  and  no  low-temperature  cracking  predicted.   This 
is  essentially  realistic  as  low-temperature  cracking  is  generally  not  a 
serious  problem  in  the  southern  United  States  where  such  low  winter 
stiffness  would  occur. 
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Figure  68.  McLeod's  chart  for  estimating  the  P.I.  of  an  asphalt  cement 
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Figure  69.   Suggested  modification  of  Heukelom's  version  of  Pfeiffer's  and 
Van  Doormal's  Nomograph  for  relationship  between  penetration, 
penetration  index  and  base  temperature.   (After  McLeod) 
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Figure  70.   Suggested  modification  of  Heukelom's  and  Klomp's  version 
of  Van  Der  Poel's  Nomograph  for  determining  modulus  of 
stiffness  of  asphalt  cements.   (After  McLeod) 
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FEDERALLY  COORDINATED   PROGRAM  OF  HIGHWAY 
RESEARCH  AND   DEVELOPMENT    (TCP) 


The  Offices  of  Research  and  Development  of  the 
Federal  Highway  Administration  are  responsible 
for  a  broad  program  of  research  with  resources 
including  its  own  staff,  contract  programs,  and  a 
Federal-Aid  program  which  is  conducted  by  or 
through  the  State  highway  departments  and  which 
also  finances  the  National  Cooperative  Highway 
Research  Program  managed  by  the  Transportation 
Research  Board.  The  Federally  Coordinated  Pro- 
gram of  Highway  Research  and  Development 
(FCP)  is  a  carefully  selected  group  of  projects 
aimed  at  urgent,  national  problems,  which  concen- 
trates these  resources  on  these  problems  to  obtain 
timely  solutions.  Virtually  all  of  the  available 
funds  and  staff  resources  are  a  part  of  the  FCP. 
together  with  as  much  of  the  Federal-aid  research 
funds  of  the  States  and  the  NCHRP  resources  as 
the  States  agree  to  devote  to  these  projects." 


FCP  Category  Descriptions 

1.  Improved    Highway   Design    and    Opera- 
tion for  Safety 

Safety  R&D  addresses  problems  connected  with 
the  responsibilities  of  the  Federal  Highway 
Administration  under  the  Highway  Safety  Act 
and  includes  investigation  of  appropriate  design 
standards,  roadside  hardware,  signing,  and 
physical  and  scientific  data  for  the  formulation 
of  improved  safety  regulations. 

2.  Reduction    of    Traffic     Congestion    and 
Improved  Operational  Efficiency 

Traffic  R&D  is  concerned  with  increasing  the 
operational  efficiency  of  existing  highways  by 
advancing  technology,  by  improving  designs  for 
existing  as  well  as  new  facilities,  and  by  keep- 
ing the  demand-capacity  relationship  in  better 
balance  through  traffic  management  techniques 
such  as  bus  and  carpool  preferential  treatment, 
motorist  information,  and  rerouting  of  traffic. 


3.  Environmental  Considerations  in  High- 
way Design,  Location,  Construction,  and 
Operation 

Environmental  R&D  is  directed  toward  identify- 
ing and  evaluating  highway  elements  which 
affect  the  quality*  of  the  human  environment. 
The  ultimate  goals  are  reduction  of  adverse  high- 
way and  traffic  impacts,  and  protection  and 
enhancement  of  the  environment. 

4.  Improved  Materials  Utilization  and  Dura- 
bility 

Materials  R&D  is  concerned  with  expanding  the 
knowledge  of  materials  properties  and  technology 
to  fully  utilize  available  naturally  occurring 
materials,  to  develop  extender  or  substitute  ma- 
terials for  materials  in  short  supply,  and  to 
devise  procedures  for  converting  industrial  and 
other  wastes  into  useful  highway  products. 
These  activities  are  all  directed  toward  the  com- 
mon goals  of  lowering  the  cost  of  highway 
construction  and  extending  the  period  of  main- 
tenance-free operation. 

5.  Improved  Design  to  Reduce  Costs,  Extend 
Life  Expectancy,  and  Insure  Structural 
Safety 


*  The  complete  7-volume  official  statement  of  the  FCP  is 
available  from  the  National  Technical  Information  Service 
(NTIS),  Springfield,  Virginia  22161  (Order  No.  PB  242057, 
price  .$45  postpaid).  Single  copies  of  the  introductory 
volume  are  obtainable  without  charge  from  Program 
Analysis  (HRD-2),  Offices  of  Research  and  Development, 
Federal    Highway    Administration,    Washington,    D.C.    20590. 


Structural  R&D  is  concerned  with  furthering  the 
latest  technological  advances  in  structural  de- 
signs, fabrication  processes,  and  construction 
techniques,  to  provide  safe,  efficient  highways 
at  reasonable  cost. 

6.  Prototype  Development  and  Implementa- 
tion of  Research 

This  category  is  concerned  with  developing  and 
transferring  research  and  technology  into  prac- 
tice, or.  as  it  has  been  commonly  identified, 
"technology  transfer." 

7.  Improved  Technology  for  Highway  Main- 
tenance 

Maintenance  R&D  objectives  include  the  develop- 
ment and  application  of  new  technology  to  im- 
prove management,  to  augment  the  utilization 
of  resources,  and  to  increase  operational  efficiency 
and  safety  in  the  maintenance  of  highway 
facilities. 
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